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1 .  INTRODUCTION 


The  phenomenon  of  second  breakdown  was  first  reported  in  diodes  by 
Tauc  and  Abraham1  and  in  transistors  by  Thornton  and  Simmons.2  The  name 
second  breakdown  has  been  chosen  because  the  first  appearance  of  micro¬ 
amperes  of  current  in  a  reverse-biased  diode  was  called  breakdown. 
Experimentally,  second  breakdown  is  generally  identified  by  an  abrupt 
decrease  of  voltage  and/or  a  sudden  increase  of  current  after  a  time  of 
relatively  constant  current  and  voltage.  If  not  prevented  by  current 
limiting  in  the  external  circuit,  second  breakdown  usually  leads  to 
device  damage.  Second  breakdown  is  of  particular  interest  to  individ¬ 
uals  in  the  power  semiconductor  community,  in  the  electromagnetic  pulse 
(EMP)  community,3  and,  increasingly,  to  those  in  the  integrated-circuit 
community. 

Since  the  first  reports,  there  has  been  a  controversy  over  whether 
electronic  or  thermal  processes  are  predominant  in  the  initiation  of 
second  breakdown.  A  comprehensive  review  and  bibliography  of  early 
papers  has  been  given  by  Schafft.1 S*  Two  more  recent  reviewers  are  Kalab* 
and  Dickhaut.6  A  few  other  recent  contributions  will  be  cited  later  in. 
this  report. 

In  1976,  Ward7  described  an  electrothermal  model  of  second  break¬ 
down.  The  computer  program  used  in  this  model  was  originally  written 
for  study  of  electrical  breakdown  in  gases,®  modified  for  study  of 


1J.  Tauc  and  A .  Abraham,  Thermal  Breakdown  in  Silicon  PS  Junctions , 
Phys.  Rev.,  108_  (1957),  936-937. 

2C.  G.  Thornton  and  C.  D.  Simmons,  A  New  High  Current  Node  of 
Transistor  Operation,  IRE  Trans.  Electron  Devices,  ED-5  (January  1958), 
6-10 . 

3L.  w.  Ricketts,  J.  E.  Bridges,  and  J.  Niletta,  EMP  Radiation  and 
Protective  Techniques,  John  Wiley  and  Sons,  New  York  (1976). 

4h .  A.  Schafft,  Second  Breakdown — A  Comprehensive  Review,  Proc. 
IEEE,  55_  (August  1967),  1272-1278. 

*B.  Kalab,  Analysis  of  Failure  of  Electronic  Circuits  for  BMP 
Induced  Signals— Review  and  Contribution,  Harry  Diamond  Laboratories, 
HDL-TR-1615  (August  1973). 

®R.  H .  Dickhaut,  Electromagnetic  Pulse  Damage  to  Bipolar  Devices, 
Circuits  and  Systems  (Magazine) ,  10_  (April  1976),  8-21. 

7 A.  L.  Ward,  An  Electro-Thermal  Model  of  Second  Breakdown,  IEEE 
Trans.  Nu cl.  Sci.,  NS-23  (December  1976),  1679-1684 . 

SW,  Borsch-Supan  and  H.  Oser,  Numerical  Computation  of  the  Temporal 

Development  of  Currents  in  a  Gas  Discharge  Tube,  J.  .Res.,  National 

Bureau  of  Standards,  67B  (January  1963),  41-60. 


electrical  processes  in  semiconductors, 9  and  finally  modified  to  include 
study  of  thermal  effects7 * 9 * ll 12  important  in  second  breakdown.  Other 
papers10-*2  on  second  breakdown  have  followed.  It  is  the  purpose  of 
this  report  to  unify  these  publications  with  additional  material  into  a 
single  comprehensive  study  of  the  basic  physical  processes  in  the  initi¬ 
ation  of  second  breakdown  in  silicon  diodes. 


2.  FORMULATION  OF  COMPUTER  PROGRAM 
2.1  Basic  Equations 

The  one-dimensional  continuity  equations  for  electrons  and 
holes,  respectively,  in  a  semiconductor  are 

e3n/3t  *  -3J_/3x  +  oJ_  +  8J+  -  R  ,  (1) 

e3p/3t  =  3j+/3x  +  oJ_  +  3J+  -  R  ,  (2) 

where  x  and  t  are  the  space  and  time  coordinates,  e  is  the  electron 
charge,  n  and  p  are  the  electron  and  hole  number  densities,  a  and  8  are 
the  ionization  coefficients  for  electrons  and  holes,  and  R  is  the  recom¬ 
bination  rate.  The  electron  and  hole  current  densities,  J_  and  J+,  are 
given  by 


J  =  ney  E  -  eD  (dn/dx)  ,  (3) 

J+  =  pep+E  +  eD+(dp/dx)  ,  (4) 

where  u_  and  w.  are  the  electron  and  hole  mobilities,  D_  and  D+  are  the 
diffusion  coefficients,  and  E  is  the  electric  field. 

Space-charge  effects  are  determined  by  the  one-dimensional 
Poisson  equation 


7A.  L.  Ward,  An  Electro-Thermal  Model  of  Second  Breakdown,  IEEE 

Trans .  Nucl .  Sci ,,  MS-23  ( December  1976),  1679-1684. 

9J .  L.  Scales  and  A.  L.  Ward,  Effects  of  Space  Charge  on  Mobility, 
Diffusion,  and  Recombination  of  Minority  Carriers,  J.  Appl.  Phys.,  39 
(February  1968),  1692-1700. 

1  °A  •  L.  Ward,  Studies  of  Second  Breakdown  in  Silicon  Diodes,  IEEE 
Trans.  Parts  and  Hybrids  Packag.,  PHP-13  (December  1977),  361-368. 

llA.  L.  Ward,  Calculations  of  Second  Breakdown,  IEEE  Trans .  Hu  cl, 
Sci.,  HS-24  (December  1977),  2357-2360. 

12A.  L.  Ward,  Doping  Profiles  and  Second  Breakdown,  Proceedings,  1979 
Electrical  Overstress/Electrostatic  Discharge  Conference,  Denver,  CO 
(1979),  25-77.  Reliability  Analysis  Center,  Rome  Air  Development 

Center,  EOS-1 ,  109-115. 
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(5) 


dE/dx  =  e(n  -  p  +  N)k£q  , 

where  <  is  the  dielectric  constant  for  the  semiconductor ,  tQ  is  the 
permittivity  of  free  space,  and  N{x)  is  the  distribution  of  net  fixed 
charge.  The  sign  convention  in  equations  (1)  through  (5)  is  chosen  so 
that  E,  J_,  and  J+  are  all  normally  positive  quantities. 

The  one-dimensional  thermal  diffusivity  equation  is 

3T  32T  3k_  /3t\2 

P®  TE  •  ™  +  kT  ^  +  TT  \Sx)  '  (6) 

where  T  is  the  temperature,  J  the  total  current  density,  p  the  semicon¬ 
ductor  density,  c  the  heat  capacity,  and  kT  the  thermal  conductivity. 

It  may  be  shown  from  equations  (1)  through  (5)  that  the  total 
current  density. 


J  =  J  +  J+  +  <eo(dE/dt)  ,  (7) 

is  a  constant  in  space.  This  merely  expresses  current  continuity  in  one 
dimension.  The  accuracy  of  the  calculations  may  be  monitored  by  the 
constancy  of  J  across  the  diode. 

2.2  Supplementary  Equations 

The  cathode  (x  =  0)  to  anode  (x  =  d)  distance  is  divided  into  M 
equal  intervals  of  width  Ax.  Initial  (t  *  0)  arrays  of  M  +  1  values 
must  be  given  for  n,  p,  and  T,  and  a  similar  array  given  for  the  fixed 
charges,  N  =  ND  -  Nft,  where  N~  is  the  net  number  density  of  donors  and 
is  the  net  number  density  of  acceptors. 

Two  boundary  conditions  are  required  for  equations  (1)  and  (2): 
the  electron  current  density  at  x  =  0  and  the  hole  current  density  at  x 
=  d.  Alternately,  the  number  densities  can  be  given  and  the  current 
densities  calculated  from  equations  (3)  and  (4).  At  present,  the  boun¬ 
dary  currents  are  constant  in  time. 

The  boundary  condition  for  equation  (5)  is  supplied  by  the 
total  voltage  across  the  diode.  The  initial  voltage  across  the  diode 
must  be  given,  and  that  for  later  times  is  determined  by  the  external 
circuit.  Optional  external  circuits  are  available,  but  for  this  report, 
the  diode,  shunted  by  a  capacitance,  C,  is  in  series  with  a  load  resis¬ 
tance,  Rg,  and  a  voltage  source,  V(t).  The  voltage  source  may  be  con¬ 
stant  in  time,  have  one  discrete  step,  or  have  an  incremental  sinusoidal 
variation.  The  last  option  allows  a  constant  dV/dt  value  to  be  closely 
approximated. 
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Equations  (1)  through  (5)  have  been  solved  by  difference  equa¬ 
tions  whose  stability  has  been  discussed  by  Borsch-Supan  and  Oser.8 
Stability  conditions  require  that  the  time  step  At  be  less  than  the 
dielectric  relaxation  time,  t  =  tceQ/o ,  where  o,  the  conductivity,  is 
given  by 


°  =  e(PUp  +  ni^)  -  J/E  . 

To  an  excellent  approximation, 9 * * *  any  uncompensated  charge  in  a  semicon¬ 
ductor  will  die  out  exponentially  with  the  time  constant  t.  A  second 
stability  condition  requires  that  the  time  step  be  less  than  twice  the 
R8C  time  constant  of  the  external  circuit.  Also,  accuracy  of  calcu¬ 
lations  (and  usually  stability)  require  that  At  be  chosen  less  than  the 
transit  time  of  the  fastest  particle  across  Ax.  A  constant  F  £  1  is 
chosen  and  we  set 


At  =  F(Ax)/v  ,  (8) 

max 

where  is  the  maximum  value  of  v_  or  v+  at  any  space  point.  A  final 

stability  condition  is  that  the  change  in  E  in  the  distance  Ax  must  be 
less  than  E  itself.  This  condition  is  usually  violated  when  E(x) 
changes  sign.  If  the  negative  field  is  of  Bhort  duration,  or  if  the 
boundary  between  negative  and  positive  field  is  moving,  the  calculations 
quickly  recover.  However,  if  negative  fields  persist,  the  calculations 
are  halted. 

It  has  been  shown13  that  the  difference  equation  represen¬ 
tations  of  equations  ( 1 )  and  ( 2 ) ,  with  only  the  conduction  current 
density  terms  included  (that  is,  neglecting  the  diffusion  components  of 
eq  (3)  and  (4))  introduces  a  diffusion  factor  that  is  a  realistic 
analogy  to  particle  diffusion.  The  magnitude  of  the  diffusion  is  given 
by 


D/y  =  (1  -  F)V/2M  ,  (9) 

where  V  is  the  diode  voltage,  F  is  defined  in  equation  (8),  and'  M  »  d/Ax 
is  defined  at  the  beginning  of  this  section.  The  derivation  assumes 
that  the  field  is  a  constant.  Most  calculations  have  been  made  Wxth  M  - 
50  or  100.  For  these  values,  D/y  is  usually  found  to  be  greater  than 

8W.  Borsch-Supan  and  H .  Oser,  Numerical  Computation  of  the  Temporal 
Development  of  Currents  in  a  Gas  Discharge  Tube ,  J.  Res,,  National 
Bureau  of  Standards,  67B  (January  1963),  41-60, 

9J.  L.  Scales  and  A.  L.  Ward,  Effects  of  Space  Charge  on  Nobility, 

Diffusion,  and  Recombination  of  Minority  Carriers,  J.  A ppl.  Phys,,  39 

(February  1968),  1692-1700 , 

13A.  L,  Ward,  Effect  of  Space  Charge  Upon  the  Transport  of  Charge 

Carriers,  J.  Appl,  Phys,,  2L  (N**ch  1964),  469-474, 
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e/kT,  the  value  required  by  the  Einstein  relation,  where  k  is  the 
Boltzmann  constant.  Choosing  larger  values  of  M  greatly  increases  the 
cost  of  calculations.  Since  diffusion  is  simulated  by  the  difference 
equations,  the  diffusion  terms  of  equations  (3)  and  (4)  are  omitted. 
The  parameter  P  is  chosen  to  be  close  to  unity  if  the  time  step 
satisfies  stability  conditions.  For  most  calculations,  0.1  <  F  <  0.9* 

Stability  conditions  require  that  At  be  less  than  the 
dielectric  relaxation  time,  t.  For  semiconductor  material  of  interest, 
t  is  of  the  order  of  picoseconds.  However,  the  thermal  time  constant 
for  most  diodes  is  microseconds,  depending  largely  upon  the  distance  to 
the  heat  sink.  Ibis  disparity  of  electrical  and  thermal  time  constants 
has  made  it  impractical  to  include  the  two  heat  conduction  terms  in 
equation  (6)  in  the  computer  program.  Therefore,  only  the  heat  source 
term,  JE,  on  the  right  side  of  equation  (6)  has  been  programmed.  The 
impact  of  this  approximation  is  negligible  for  the  time  scale  imposed  by 
computation  cost. 

To  solve  equations  (1)  through  (6),  the  cycle  is  as  follows. 
The  initial  values  of  n(x),  p(x),  and  N(x)  are  used  to  calculate  dE/dx 
from  equation  (5).  Setting  E(0)  *  0,  /Edx  is  evaluated.  Then  the  field 
at  each  grid  point  is  adjusted  by  a  constant  increment  in  order  that 
/Edx  »  V,  the  desired  voltage  across  the  diode.  Knowing  E(x),  all 
field-dependent  quantities,  such  as  J_(x),  J+(x),  a(x),  &(x),  u_(x),  and 
U+(x) ,  may  be  calculated.  These  quantities  are  then  used  to  calculate 
new  values  of  n(x),  p(x),  and  T(x)  at  the  time  At,  with  equations  (1), 
(2),  and  (6),  respectively.  The  average  current  density  is  calculated 
from  equation  (7)  and,  with  the  given  diode  area,  the  total  current 
flowing  in  the  diode  is  determined.  This  current  then  is  used  to  obtain 
an  adjusted  voltage  across  the  diode  due  to  the  action  of  the  external 
circuit,  as  for  example  the  IRS  drop  across  the  external  resistor,  Rg. 
This  completes  one  time  cycle. 

The  form  of  the  recombination  rate  in  equations  (1)  and  (2)  in 
the  computer  program  was  chosen  to  be 

R  «  eR'(np  -  n2)/nA  ,  (10) 

where  n^  is  the  carrier  density  of  intrinsic  material  and  R’  is  the 
volume  rate  of  electron  hole  recombination.  This  form  for  the  recombi¬ 
nation  rate  was  chosen  for  its  versatility.  The  usual  expressions 
using  minority  carrier  lifetimes,  tr  and  t  ,  are  much  more  complex  in 
that  they  are  functions  of  the  doping  density.  Equation  (10)  reduces  to 
the  usual  minority  carrier  lifetime  equation  when  the  excess  minority 
carrier  number  density  is  small  compared  to  the  doping  density.  The 
result9  is  that 

~~  VTTT Scales  and  A.  L.  Ward ,  Effects  of  Space  Charge  on  Mobility, 
Diffusion,  and  Recombination  of  Minority  Carriers,  J.  A ppl.  Phys,,  39 
(February  1968),  1692-1700 . 


(11) 


ni/R,po 


n./R'n 
i  o 


where 


nQ  and  pQ  are  the  majority  carrier  densities, 
carrier  lifetime  is  a  function  of  n,  p,  n_,  and  p  as  desired. 


Thus,  the  minority 

uat  lici  uicuiuic  19  ci  luiivuuu  ui  ii  t  ^  g  f  auu  ^  as  desired.  R'  needs 
only  to  be  adjusted  for  the  quality  of  the  material  or  to  simulate  gold 


doping,  for  example.  Note  from  equation  (10)  that  thermal  generation  of 
carriers  exceeds  the  recombination  when  the  np  product  is  less  than  n2. 


3.  MATERIAL  PARAMETERS  FOR  SILICON 

The  most  critical  material  parameters  for  semiconductor  breakdown 
studies  are  the  ionization  coefficients.  The  ionization  coefficients, 
defined  as  the  number  of  electron-hole  pairs  created  by  collisions  per 
centimeter  of  travel,  are  given  by 

a  =  Ai  exp(-A2/E)  ,  (12) 

&  -  B!  exp(-B2/E)  ,  (13) 

where  Ai  and  B]  are  empirical  constants.  Most  of  the  calculations 
reported  herein  use  those  measured  by  Van  Overstraeten  and  DeMan.14  The 
computer  program  allows  different  empirical  constants  to  be  used  above  a 
fixed  value  of  E.  The  constants  used  were  generally  chosen  as  follows: 


Ai  -  7.03  x  10s 
A2  =*  1.23  x  106 
B i  =  1.582  x  106 
B i  =  6.75  x  105 
B2  =  2.036  x  106 
B2  -  1.693  x  106 


for  all  E, 
for  all  E, 

E  <  4  x  10s  V/cm, 
E  >  4  x  105  V/cm, 
E  <  4  x  105  V/cm, 
E  >  4  x  10s  V/cm. 


The  temperature  variation  of  the  ionization  coefficients  was  based 
on  the  measurement  of  Crowell  and  Sze . 1 5  They  found  variation  for  both 

14r.  Van  Overstraeten  and  H.  DeMan,  Measurement  of  the  Ionization 
Rates  in  Diffused  Silicon  p-n  Junctions ,  Solid-State  Electron.,  13_ 
(1970),  583-608. 

1  r,  Crowell  and  S.  M.  Sze,  Temperature  Dependence  of  Avalanche 
Multiplication  in  Semiconductors,  Appl.  Phys.  Lett.,  9_  (September  1966), 
242-244. 
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A i  and  A 2  with  temperature,  but  the  variation  of  A2  was  small  and 
irregular,  while  that  of  Ai  was  larger  and  more  regular.  Therefore, 
only  the  variation  of  Ai  was  used.  Ai  is  reduced  by  0.2  percent  for 
each  degree  increase  above  300  K.  No  measurements  are  available  for  the 
temperature  dependence  of  the  hole  ionization  coefficients,  so  an  equal 
decrease  was  assumed. 

A  second  material  parameter  of  importance  is  the  variation  of 
carrier  mobility  with  field,  temperature,  and  doping  level.  The  vari¬ 
ation  of  the  carrier  mobilities  with  field  is  given  by  equations  of  the 
form 


V  =  Vi0(l  -  CiE)  ,  E  <  C4 


(14) 


U  =  C2E-1/2(l  -  C3E“3/2)  ,  E  >  C4  ,  (15) 


where  uQ  is  the  low  field  mobility  and  the  constants  may  be 
considered  empirical.  Each  maximum  velocity  is  also  bounded  by  maximum 
values,  which  are  included  in  the  input  data.  The  low  field  mobilities 
were  taken  from  Sze16  and  are  *  1500  cm/v-s  for  electrons  and  iiQ  * 
600  cm/V-s  for  holes.  There  is  a  further  option  in  the  computer  program 
to  interpose  a  transition  range  between  two  given  field  values  where 
equations  (14)  and  (15)  are  used.  The  expression  is  a  power  series  in  E 
and  the  parameters  are  calculated  internally  to  insure  continuity  of  v 
and  dv/dE  at  the  end  points.  Several  expressions  have  been  used  in  the 
calculations,  but  the  expressions  used  most  often  give  the  velocity 
versus  field  curves  given  in  figure  1.  The  low  field  mobility  decreases 
with  increasing  temperature.  The  temperature  variation  given  by  Sze18 
is  that  the  electron  mobility  varies  as  T-2 *5  and  the  hole  mobility  as 
Tf^*7.  The  variation  of  the  saturated  drift  velocities  with  temperature 
is  less  than  for  the  low  field  mobility.  The  variation  chosen  was  based 
on  the  measurements  of  Duh  and  Moll17  and  of  Canali  et  al.18  The  vari¬ 
ation  of  drift  velocity  with  field  at  elevated  temperatures  is  also 
shown  in  figure  1 . 


16S.  M.  Sze,  Physics  of  Semiconductor  Devices,  tfiley-Int ersc ience. 
New  York  (1969),  39-41,  57-59,  121-126. 

17C.  Y .  Duh  and  J.  L.  Noll,  Electron  Drift  Velocity  in  Avalanche 
Silicon  Diodes,  IEEE  Trans.  Electron  Devices,  ED-14  (January  1967),  46- 
49. 

18C.  Canali,  G.  Majni,  R.  Minder,  and  G.  Ottaviani,  Electron  and  Bole 
Drift  Velocity  Measurements  and  Their  Empirical  Relation  to  Electric 
Field  and  Temperature,  IEEE  Trans.  Electron  Devices,  ED-22  (November 
1975),  1045-1047. 
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Figure  1.  Temperature  variation  of  (a)  electron  and  (b)  hole 
velocities  as  function  of  electric  field. 


There  is  also  an  option  in  the  computer  program  to  vary  the  drift 
mobilities  as  a  function  of  doping  level.  The  variation  follows  that 
given  in  Sze.lfc  However,  measurements  of  saturation  velocities  as  a 
function  of  doping  level  and  of  temperature  have  not  been  made.  Just  as 
the  variation  of  the  saturation  velocities  with  temperature  is  less  than 
that  of  the  low  field  mobilities,  one  should  expect  that  the  variation 
with  doping  level  should  be  less  for  the  saturated  velocities.  For  this 
reason,  the  option  to  vary  the  mobilities  with  doping  level  has  not  been 
used  for  second  breakdown  studies  where  saturation  velocities  are 
predominant. 

The  variation  of  the  intrinsic  density  with  temperature  was  chosen 
as 


n2(T)  »  2  *  1Oz0(T/3OO)3  exp[ (Eg/k)(300-1  -  T-1)]  ,  (16) 

where  Eg  is  the  band  gap  energy  of  silicon.  The  variation  of  the 
injected  (thermally  generated)  current  density  with  temperature  was 
assumed  to  be  the  same  as  for  n|(T),  since  np  ■  n?  and  the  majority 
carrier  density  is  fixed  by  the 1 doping  level.  Since  the  effective 
thermally  injected  current  is  determined  by  both  thermal  generation  and 
particle  diffusion,  the  calculated  saturated  current  density  for  the 
ideal  diode  equation  is  in  poor  agreement  with  measurements.  For  this 
reason  a  value  of  room  temperature  thermal  injection  current,  JQ,  of  1  x 
10"7  A/cm 2  was  chosen  for  each  boundary.  The  minority  electron  current 
density  injected  at  the  cathode  boundary  is  designated  JRO  and  the 
minority  hole  current  density  injected  at  the  anode  boundary  is 

M.  Sze,  Physics  of  Semiconductor  Devices,  Wiley-Interscience, 
See  York  (1969),  39-41,  57-59,  121-126. 
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designated  J  .  As  stated,  Jno  =  J_Q  =  JQ  =  1  x  10“7  k/cm  was  chosen  at 
room  temperature  and  was  varied  with  temperature  as  n~(T),  as  given  in 
equation  (16),  for  higher  temperature  runs. 

The  other  silicon  material  properties  are  taken  from  Sze.16  They 

are 


dielectric  constant 

11.8 

density 

2.328  g/cm 

energy  gap  at  300  K 

1.12  eV 

minority  carrier  lifetime 

2.5  x  10~3 

specific  heat 

0.7  J/gC 

thermal  conductivity  at  300  K 

1.45  W/cmC 

4.  PROCEDURES  AND  DISCUSSION 

Before  we  present  the  results  of  the  computer  calculations  it  is 
useful  to  further  discuss  the  formulation,  the  procedures,  and  the 
limitations  of  the  computer  program. 

First,  consider  the  continuity  equations,  (1)  and  (2).  If  3n/3t  * 
3p/3t  *  0,  then  a  steady-state  condition  prevails.  Hie  steady-state 
equation  for  electrons,  neglecting  recombination,  is 

3J_/3x  -  +  BJ+  .  (17) 

Since,  for  silicon,  the  ionization  constant  for  electrons  is  of  an  order 
of  magnitude  greater  than  for  holes,  we  may,  as  an  approximation, 
neglect  6  compared  to  ou  Then,  for  a  constant  field,  one  obtains 

J_(x)  -  JQ  exp (ax)  . 

This  shows  that  for  a  PIN  diode,  where  the  field  is  constant,  we  should 
expect  the  electron  current  density  to  increase  exponentially  across  the 
diode  at  low  currents,  where  space  charge  can  be  neglected.  This  ap¬ 
proximation  will  be  shown  to  be  surprisingly  good  for  the  high  field 
region  of  other  diodes,  also.  Since  dE/dt  *  0  in  the  steady  state, 
J+(x)  ■  J  -  J_(x),  from  equation  (7). 

16S.  N.  Sze,  Physics  of  Semiconductor  Devices,  Wiley-Intersclenoe, 
Hew  York  (1969),  39-41,  S7-59,  121-126 . 
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The  steady-state  solution  of  the  continuity  equations  determines  the 
static  current  density  versus  voltage  characteristic  for  the  device 
The  static  characteristic  separates  the  high-voltage  region,  where 
currents  increase  with  time,  from  the  low-voltage  region,  where  currents 
decrease  with  time.  In  dynamic  calculations  with  a  constant  applied 
voltage  and  neglecting  heat  generation,  a  point  of  the  static  character¬ 
istic  is  attained  after  a  period  of  time. 

Since  steady-state  calculations  are  much  simpler  and  cheaper  than 
dynamic  calculations,  it  may  be  suggested  that  they  should  be  used  in 
conjunction  with  dynamic  calculations.  This  was  done  by  the  author  in 
studies  of  breakdown  in  gases.  However,  the  difference  in  accuracy 
between  the  two  programs  made  the  static  program  of  little  use  in  sup¬ 
porting  the  dynamic  program.  Therefore,  the  static  program  was  not 
modified  to  supplement  the  dynamic  program  for  semiconductor  calcu¬ 
lations. 

A  second  useful  approximation  to  the  continuity  equations  considers 
the  case  where  3j/3x  is  small  and  may  be  neglected.  We  then  have  the 
approximation  that 


e3n/dt  =  aj_  *  cmev 


(18) 


Again  we  consider  the  constant  field  case  so  that  a  and  v  are 
constant.  The  solution  to  this  equation  is  that 


J_(t)  =  J_(o)  exp(ovt)  . 

This  shows  that  the  current  in  a  diode  will  initially  grow  exponentially 
with  a  time  constant  of  (av)-^.  This  technique  has  been  used  with  much 
greater  accuracy  by  Holway.19 


It  is  also  useful  to  consider  some  approximate  solutions  to  the 
thermal  diffusivity  equation,  equation  (6).  For  times  short  compared  to 
the  thermal  time  constant  (t^,  defined  later),  the  terms  of  equation 
(6),  including  the  spatial  derivatives  of  T  with  x,  are  small,  and 


dT  a  JE 
dt  pc  ’ 

This  expression  defines  an  adiabatic  calculation. 


(19) 


H .  Holway,  Jr.,  Electron-Hole  Avalanches  with  Constant 
Ionization  Coefficients,  IEEE  Trans ,  Electron  Devices,  ED-26  (June 
1979),  991-993. 
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For  the  one-dimensional  approximation,  where  the  heat  source  is 
planar  and  all  heat  flows  to  a  heat  sink  at  a  distance,  L,  which  is  at  a 
temperature  TQ,  the  junction  temperature,  Tj ,  as  a  function  of  t  is 

T =  To  +  ep[l  -  exp(-t/Tth)]  ,  (20) 

where  0  is  the  thermal  resistance  and  P  is  the  power  dissipated. 
Thermal  resistance  is  defined  from  the  junction  temperature  in  the 
steady  state, 

T.  =  T  +  0P  .  (21) 

3  ° 

The  thermal  time  constant  is  found  to  be 

Tth  =  '  (22) 

and 

0  =  *  (23) 

where  A  is  the  diode  area.  Combining  equations  (22)  and  (23), 


8/ Tth  =  ( pcAL)~ 1 


(24.) 


The  rate  of  temperature  rise  is  found  from  equation  (20)  to  be 


dTj/dt  =  (Pe/T^)  exp(-t/xth)  .  (25) 

Since  P  =  IV,  we  find,  with  the  use  of  equation  (24), 


dTj/dt  *  (IV/ALpc)  exp(-t/Tth) 
=  ( JE/pc)  exp(-t/xth)  . 


(26) 


Comparing  equation  (26)  with  (19),  we  see  that  the  heat  conduction 
processes  reduce  the  rate  of  temperature  rise  by  the  exponential  time 
function  of  t^. 


One  further  useful  thermal  coefficient  is  the  thermal  diffusivity, 

Dth  *  analogous  to  the  electrical  diffusion  coefficient.  The 

distance  that  heat  diffuses  in  a  given  time ,  t,  depends  on  the  geometry, 

but  in  general  is  given  by  A)  .t. 

th 


It  was  mentioned  in  section  2  that  stability  conditions  require  that 
At  be  less  than  the  dielectric  relaxation  time,  r  =  <tQ/o.  Since  t  is 
inversely  proportional  to  the  number  densities,  computations  increase  in 
cost  for  higher  doping  levels.  Fbr  this  reason,  only  the  high  field 
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regions  are  simulated.  Generally,  the  maximum  doping  at  the  boundaries 
is  set  at  about  1  x  1017  cm”3  or  greater.  This  allows  the  field  at  the 
boundaries  to  be  well  below  that  for  which  the  velocity  starts  to  satu¬ 
rate,  which  in  turn  is  more  than  an  order  of  magnitude  below  the  ava¬ 
lanching  field.  Calculations  made  with  higher  doping  levels  at  the 
boundaries  show  little  difference  if  the  above  criteria  are  observed. 
With  the  use  of  a  practical  number  of  distance  steps  (M  =  50  to  100),  it 
is  often  found  that  negative  fields  are  calculated  near  the  boundaries 
of  high  doping  levels.  In  this  case,  the  negative  fields  are  eliminated 
by  reducing  the  doping  levels  at  the  boundaries.  Then,  at  higher 
currents,  with  their  resulting  space-charge  field  distortion,  the  doping 
levels  may  again  be  increased  to  observe  the  above  criteria. 

The  cost  of  running  the  program  on  an  IBM  370/168  computer  varies 
with  doping  density,  diode  widths,  and  other  factors.  A  typical  cost  is 
$10  per  nanosecond  of  problem  time.  This  precludes  problem  times  over  a 
few  tens  of  nanoseconds.  However,  it  will  be  shown  that  longer  problems 
may  often  be  treated  in  reasonable  approximation. 

Or,w  important  factor  in  junction  diode  physics  is  the  punch-through 
factor,  PTF.  A  diode  is  termed  punched  through  if  the  field  extends 
from  the  junction  completely  through  the  low  doping  region  into  the  high 
doping  boundary  region.  Some  investigators  term  this  effect  in  diodes 
reach-through  and  use  the  term  punch-through  for  transistors.  However, 
the  term  punch-through  is  used  unanimously  for  IMPATT  diodes,  and  is  so 
used  in  this  report.  The  PTF  is 
considered  equal  to  the  ratio  of  ( 1 ) 
the  distance  from  the  junction  to 
where  the  field  extrapolates  to  zero 
at  breakdown  to  (2)  the  maximum  de- 
pletable  width,  W,  i.e.,  the  width 
of  the  low  doping  region.  There¬ 
fore,  PTF  *  1  for  a  diode  whose  high 
field  region  just  reaches  the  high 
doping  boundary  region  at  break¬ 
down.  For  the  same  width  diode,  the 
PTF  >  1  for  a  diode  doped  more 

lightly  and  the  PTF  <  1  for  a  more 
heavily  doped  diode.  The  definition 
of  PTF' s  is  illustrated  in  figure 
2.  Punch-through  factors  may  be 
extended  to  two-sided  diodes  and 
linearly  graded  diodes,  using  the 
appropriate  nonlinear  extrapolation 
in  the  latter. 


Figure  2.  Punch-through  factor 
defined.  Curved  dashed  lines 
show  field  distribution  of  non¬ 
ideal  device. 
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Most  junction  diodes  are  one-sided;  that  is,  one  side  has  a  much 
lower  doping  level  than  the  other.  One-sided  diodes  either  may  be 
abrupt,  with  steps  in  the  doping  profile,  or  diffused,  where  the  doping 
changes  gradually.  Diffusion  profiles  may  be  given20  either  by  the 
complementary  error  function  or  by  the  Gaussian  distribution  according 
to  the  experimental  procedure  used.  Either  distribution  gives  approx¬ 
imately  an  exponential  falloff  at  a  distance  from  the  surface.  Most 
calculations  have  been  made  for  exponential  decreases  in  doping  from 
each  boundary  added  to  a  constant  N-  or  P-type  doping.  Hie  diodes  are 
designated  N  type  for  a  P+NN+  doping  profile  and  P  type  for  an  N+PP+ 
doping  profile.  Two-sided  diodes,  either  abrupt  or  linearly  graded, 
will  also  have  high  doping  regions,  decreasing  exponentially  from  the 
boundaries . 


Many  thermal  models  of  second 
breakdown  use  the  variation  of  re¬ 
sistivity  with  temperature.  This 
variation  is  shown  in  figure  3. 
However,  resistivity  no  longer  is 
meaningful  when  the  mobility  becomes 
a  function  of  field.  As  shown  in 
figure  1,  this  occurs  for  fields 
above  about  1  x  103  V/cra.  In  sili¬ 
con,  the  velocity  of  electrons 
becomes  saturated  at  fields  only 
about  one-tenth  of  the  avalanching 
field.  Therefore,  the  use  of  re¬ 
sistivity  variation  with  temperature 
to  study  second  breakdown  is  inaccu¬ 
rate.  Further,  many  authors  term 
the  decrease  of  resistivity  with 
temperature  as  a  negative  resis¬ 
tivity.  This  is  a  misnomer.  Cur¬ 
rent  contraction  follows  directly 
from  a  true  negative  resistance,  but 
not  from  a  negative  temperature 
coefficient.  This  will  be  discussed 
further  in  later  sections. 


TEMPERATURE  (C) 

Figure  3.  Typical  variation  of 
resistivity  with  temperature. 
Parameter  is  doping  level  in 


5.  TYPICAL  REVERSE- 3IAS  STATIC  CHARACTERISTICS 

Quasi-static  current  density  versus  voltage  characteristics  have 
been  calculated  point  by  point  for  typical  N-  and  P-type  diodes.  Hie 
prefix  quasi  has  been  used  to  indicate  some  slight  current  decrease  as 
the  temperature  rises.  Initial  conditions  for  the  calculations  are 

20A.  S.  Grove,  Physics  and  Technology  of  Semiconductor  Devices ,  John 
Wiley  and  Sons,  New  York  (1967),  44-52,  194-201, 
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either  space-charge  free,  i.e,,  the  mobile  charges  equal  the  fixed 
charge  at  all  points,  or  the  equilibrium  condition  of  a  previous  calcu¬ 
lation.  Generally,  quasi -equilibrium  is  obtained  in  a  few  tenths  of  a 
nanosecond  of  problem  time. 


The  results  are  presented  more  clearly  by  dividing  the  current 
density  (J)  range  into  four  regimes  as  defined  below.  These  regimes 
will  be  exemplified  for  an  N-type  diode  at  room  temperature  (300  K). 
The  total  diode  width  is  4  un.  Each  of  the  P+  and  N+  regions  is  0.5  pm, 
and  the  3-pm  N-region  is  doped  to  5  *  1015  cm-3.  The  boundary  between 
the  N  and  N+  region  is  arbitrarily  set  at  where  the  dopin  •  level,  N,  is 
doubled.  The  device  is  punched  through  at  breakdown.  A  doping  level  of 
7  x  i015  is  required  for  PTF  =1. 


The  regime  below  J  =  1  x  io3  A/cm2  for  this  diode  is  designated 
the  space-charge-free  (SCF)  or  the  single-avalanche  regime.  Since  JQ  = 
1  x  io'7  A/cm  for  this  calculation,  the  current  multiplication  factor 
for  this  region  reaches  a  maximum  of  1  x  101®.  In  the  SCF  regime,  the 
current  increases  exponentially  with  voltage  as  shown  in  figure  4.  The 
data  are  fitted  to  the  equation 

JSCF  =  J10  exP[(V  *  VBl)/Vl°]  '  {27) 


where  Jjo  =  0.4  A/cm2,  VB1  =  79  V,  and  Vio  =  0.41  V. 

The  second  regime  is  the  space-rharge-limi ted  (SCL)  regime  and 
extends  from  1  x  103  to  8  x  103  A/cm2.  For  a  punched -through  diode,  the 
best  fit  to  this  regime  has  been  found  to  be  that  J  varies  as  the  square 
of  the  voltage  above  Vg1.  This  is  shown  for  this  diode  in  figure  5  for 
three  values  of  V0, .  TTie  PTF  of  this  diode  is  about  1.4  and  the  best 
fit  is  between  slopes  M  =  2  and  M  =  1.6.  For  higher  PTF's,  the  best  fit 
is  usually  M  =  2.  The  dashed  curve  labeled  Ji  in  figure  4  is 

Jl  =  50( V  -  79 ) 2  A/cm2  . 

The  current  density  falls  below  the  exponential  increase  with  voltage 
because  of  the  space  charge  of  the  mobile  electrons.  The  space  charge 
decreases  the  magnitude  of  dE/dx,  increasing  the  field  away  from  the 
junction.  When  the  field  near  the  NN+  junction  reaches  the  avalanche 
field,  the  current  rises  above  the  SCL  value.  This  third  J  regime  is 
designated  the  double-avalanche  regime  and  extends  from  J  =  6  x  10^  to 
2.3  x  io4  A/cm2.  In  the  double -avalanche  regime  the  current  exceeds  the 
SCL  current  by  another  exponential  increase  with  voltage.  This  extra 
current  is  labeled  J2  in  figure  4  and  is  given  by 

j2  =.  4.6  x  10-1  exp[  (V  -  79)/1 .69]  . 


VOLTAGE  (V) 


Figure  4.  Static  current  density 
versus  voltage  characteristic  for 
4-ym  N-type  diode  doped  to  5  x 
1015  cm  ^ .  Points  are  computed 
data  and  curves  are  analytic  fits 
as  described  in  text. 


VOLTAGE  (V) 


Figure  5.  Current  density  as 
function  of  overvoltage  for 
various  assumed  breakdown  volt¬ 
ages,  vBi*  Slope  of  each 
straight  line  js  indicated. 


Hie  current  in  the  SCL  current  regime  can  be  rewritten  as 

JSCL  "  (V  "  VBl)2/RscVB1A  '  <28) 

where  A  is  the  diode  area  and  R  is  the  space-charge  resistance.  V_1 
is  included  in  the  denominator  of  equation  (28)  in  order  to  retain  the 
correct  dimensions  for  Rsc»  It  will  be  shown  later  that,  for  diodes  not 
punched  through  at  breakdown,  the  current  is  proportional  to  the  excess 
voltage  above  V01 .  Both  cases  may  be  included  in  a  single  expression: 

JSCL  =  (v  -  Vbi )1+m/R8cvB1A  »  <29) 


where  m  *  0  for  PTF  <  1  and  m  *  1  for  PTF  >  1 .  For  m  »  0,  Rsc  is  often 
called  the  surge  resistance,  but  space-charge  resistance  is  more  inform- 
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ative.  For  m  *  1,  the  justification  of  calling  Rsc  the  space-charge 
resistance  follows  from  the  near  continuity  of  Rsc  as  the  RTF  changes 
from  below  to  above  unity. 

Hie  J-V  characteristic  reaches  infinite  slope  (zero  resistance) 
at  J  =  2.3  x  101*  A/cm  and  V  =  94.5  V.  The  regime  above  this  current 
density  is  the  negative-differential-resistance  (NDR)  regime,  and  it  is 
usually  impossible  to  calculate  steady-state  conditions  in  this 
regime.  The  points  shown  in  figure  4  for  the  NDR  regime  are  average 
values  for  dynamic  characteristics  with  increasing  and  with  decreasing 
currents. 


5.2  Diode  Profiles 


Figure  6  shows  the  electron-  and  hole-density  distributions  for 
currents  in  the  SCL  and  double-avalanche  regimes.  Figure  7  shows  the 
corresponding  field  distributions.  Note  that  in  the  double-avalanche 
regime  the  mobile  electron  density  exceeds  the  sum  of  the  doping  level 
and  the  hole  density.  This  results  in  a  change  in  the  sign  of  dE/dx  and 
a  field  minimum.  The  existence  of  double  avalanching  is  shown  more 
dramatically  in  figure  8,  where  aJ_  and  0J+  are  plotted  across  the 
diode.  aJ_  is  the  rate  of  increase  in  charge  density  due  to  ionization 
by  electrons,  as  may  be  seen  from  equation  (1).  Similarly,  0J+  is  the 
increase  due  to  ionization  by  holes.  Although  a  exceeds  0  by  a  factor 
of  ~10,  0J+  usually  exceeds  ctJ_  at  some  position,  since  generally  J+ 
exceeds  J_  sufficiently  in  the  p  doping  region.  Plots  of  J_(x)  and 
J+(x)  are  essentially  ^proportional  to  n(x)  and  p(x),  respectively, 
except  in  the  N+  and  P+  doping  regions.  The  proportionality  follows 
since  the  particle  velocities  are  saturated  in  the  N  region.  Figure  9 
shows  plots  of  aJ_  at  the  PN  junction  as  a  function  of  current 
density.  Also  shown  is  the  value  at  the  NN+  junction,  defined  to  be 
where  N(x)  is  twice  the  N  region  doping.  txJ_  increases  approximately 
linearly  with  J  at  the  PN  junction,  since  the  junction  field  is  almost 
constant  as  J  increases.  However,  aJ„  increases  with  the  3.3  power  of  J 
at  the  NN+  junction  as  a  result  of  the  increasing  field.  Ebr  J  >  1.5  x 
101*  A/ cm,  aJ_  is  greater  at  the  NN+  junction  than  at  the  P+N  junction. 


When  a  =  0,  the  criterion  for  breakdown  is 

Jq  a  dx  *  a^ygd  *  1  .  (30) 

This  criterion  is  still  often  used  as  a  breakdown  criterion  in  silicon, 
although  a  ranges  approximately  10  times  greater  than  0  in  silicon. 
Plots  of  a^yjd  and  0AVEd  as  a  function  of  J  are  shown  in  figure  10.  It 
is  seen  that  is  much  greater  than  unity.  A  good  approximation  to 

breakdown  is  that  the  product  of  ot^ygd  and  0AV£d  is  equal  to  one. 
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Figure  6.  Distribution  of  number 
density  for  diode  of  figure  4. 
Doping  profile  is  given  by  solid 
lines,  electrons  by  dashed  lines, 
and  holes  by  dotted  lines.  Pa¬ 
rameter  is  current  density. 


Figure  7.  Field  distributions  cor¬ 
responding  to  current  densities  in 
figure  6.  Diode  voltage  is  param¬ 
eter. 


Figure  8.  Distribution  of  ionization  rates  across  ti-type  diodes 
(a)  ionizations  by  electron  collisions,  and  (b)  ionizations  by 
hole  collisions.  Parameter  is  current  density  in  A/cm2. 
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rates  as  function  of  current 
density.  Solid  lines  are 
drawn  with  a  slope  of  1.0. 
Rates  for  NN+  junction  (dashed 
curves)  are  not  maxima,  except 
at  highest  current  densities. 


Figure  10.  (a)  Integrals  of  a, 

8,  and  a8  across  diode  as  func¬ 
tion  of  current  density  for  P- 
type  diode.  (b)  Dynamic  char¬ 
acteristic  of  diode,  calculated 
with  dV/dt  «  2  x  io 10  V/s. 


Complementary  N-type  and  P-type  diodes  have  the  same  doping 
levels,  diode  width,  and  area.  Calculations  of  quasi -static  character¬ 
istics  for  the  complementary  p-type  diode  with  a  doping  level  of  5  x 
1015'cra“3  and  a  4-ym  width  yield  a  first-breakdown  voltage  VB1  of  80  V 
compared  to  79  V  for  the  N-type  diode.  The  maximum  voltage  was  101.5  V 
at  1.6  x  io1*  A/cm  compared  to  94.5  V  at  2.3  x  io4  A/cm  for  the  N-type 
diode.  The  maximum  voltage  will  be  designated  the  second  breakdown 
voltage,  VB2"  The  difference  in  Vfi2  f°r  the  two  diodes  is  7  V  compared 
to  a  difference  of  1  V  in  Vg1 .  This  indicates  that  the  p-type  diode  has 
a  greater  space-charge  resistance  than  does  the  N-type.  From  equation 
(28),  with  an  assumed  area  of  1  x  IO-4  cm,  the  calculated  value  of  Rsc 

for  the  N  diode  is  2.5  0  and  for  the  P  diode,  6.6  ft. 

For  the  same  4-itn  width  diode,  a  doping  level  of  1  x  io3^  cm”3 
gives  a  PTF  of  about  0.7;  that  is,  the  field  at  breakdown  extrapolates 
to  zero  at  about  0.7  of  the  depletable  width.  For  the  p  diode,  V01  ■  62 
V  and  VB2  *  101  V  at  1.5  x  io4  A/cm2.  Note  that  whereas  VB1  is  17  V 
lower  than  that  for  the  same  width  diode  doped  to  5  x  io3^  cm”3,  the 
values  of  VB2  are  essentially  the  same.  The  static  characteristic  of 
this  diode  is  shown  in  figure  11.  This  time,  a  linear  scale  is  used  to 

show  that  the  current  rises  linearly  with  voltage  above  VB1 .  A  surge 

resistance  can  be  defined  in  the  conventional  manner  and  is  determined 
to  be  23.8  0  for  a  diode  of  1  x  io”4  cm2  area. 
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It  has  been  shown  that  good  analytical  fits  to  the  computer 
data  can  be  made  with  J  «  (v  -  VBj  )  for  diodes  with  PTF  <  1  and  with 
j  «  (v  -  V  )2  for  PTF  >  1.  How  good  are  the  fits  for  PTF  z  1?  For  the 
4-pm  width8 diode ,  a  doping  of  7  x  1015  cm"3  gives  PTF  z  1.  Figure  12 
shows  a  dynamic  J-V  characteristic  calculated  for  the  P  diode.  Hie 
calculation  was  made  with  an  applied  voltage  rising  at  the  rate  of  1  x 
10*®  V/s.  Hiis  technique  greatly  reduces  the  cost  of  calculation.  It 
will  be  shown  in  a  following  section  that  this  rate  of  rise  gives  a 
maximum  error  of  about  1  V  at  any  current  density  up  to  VB2.  Also  shown 
in  figure  12  is  the  fit  of  three  analytic  expressions.  Dotted  curve  "O" 
is  the  exponential  variation  of  current  in  the  SCF  region  of  J  <  5  x  102 
A/cm.  The  fit  to  equation  (27)  gives  JlO  =  1.0  A/cm,  *  73.4,  and 

Vio  =  0.36  V.  Solid  curve  "1"  is  the  best  fit  to  J  «  V  -  V^,  and  the 
dashed  curve  ”2"  is  the  best  fit  to  J  <*  (v  -  VB1)2*  The  analytic  ex¬ 
pressions  are 

J  =  410(V  -  77)  and 

j  w  10(V  -  67) 2  . 


Figure  11.  Static  current  den¬ 
sity  versus  voltage  character¬ 
istic  for  P-type  diode.  Diode 
width  is  4  ym  and  doping  level 
is  1  x  1016  cm-3. 


Figure  12.  Characteristic 
for  diode  with  PTF  “  1 . 

Individual  points  were  cal¬ 
culated  with  dV/dt  *  1  x 

1010  V/s.  Equations  for 
analytic  curves  are  given  in 
text. 
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It  may  be  noted  that  VB1  from  equation  (27)  is  intermediate  between 
V'  and  V£r  Again  assuming  A  *  1  x  10"4  cm2,  Rsc  calculated  from 
elation  (29)  is  found  to  be  24.4  fi  for  m  *  0  and  14.9  0  for  m  *  1. 
Considering  that  PTF  ■  1  is  a  worst  case  for  each  expression,  there  is  a 
reasonable  fit  by  both  expressions  in  the  SCL  region,  which  extends  to 
about  8  x  103  A/cm2.  An  additional  exponential  terra  extends  the  fit 
into  the  double  avalanche  as  shown  previously.  As  expected,  curve  "1" 
is  the  better  fit  at  lower  J  values,  since  the  depleted  region  widens  a 
little,  and  curve  "2"  is  a  better  fit  at  higher  J  values,  since  no 
further  widening  of  the  depletion  region  is  possible.  Fbr  engineering 
analysis,  it  seems  satisfactory  to  use  J  «  V  -  VB1  for  PTF  £  1  and 
j  «  (v  -  VJ2)2  for  PTF  >1*  If  a  greater  accuracy  is  desired,  an 
intermediate  value  of  m  between  0  and  1  may  be  used  for  PTF  =  1 . 

6.  DYNAMIC  CHARACTERISTICS 


Constant  dV/dt 


A  series  of  calculations  has  been  made  with  varying  constant 
rates  of  rise  of  the  applied  voltage,  dV/dt.  The  initial  charge  distri¬ 
butions  were  chosen  to  give  space-charge  neutrality  and  the  initial 
voltage  somewhat  below  VB1 .  Calculations  were  made  for  the  complemen¬ 
tary  diodes  discussed  in  section  5:  d  =  4  ym  and  doping  of  5  x  1015 
cm-3.  the  dynamic  characteristics  are  compared  with  the  static 
characteristic  in  figure  13.  The  initial  transients  from  the  rather 
high  initial  currents  are  not  shown,  since  they  would  tend  to  confuse 
the  figure.  The  variation  of  voltage  with  time  for  four  values  of  dV/dt 
is  shown  in  figure  14.  In  order  to  allow  the  useful  information  to  be 
displayed  on  a  shorter  time  scale,  zero  time  in  this  plot  was  chosen  to 
be  the  first  time  that  Vg1  (80  V)  was  attained.  Both  figures  13  and  14 


istics  calculated  with  labeled  Figure  14.  Variation  of  voltage 

rates  of  voltage  rise  compared  with  time  for  diode  of  figure 

with  calculated  static  char-  13.  Time  is  measured  from  attain- 

acteristic.  The  4-ym  P-type  raent  of  VB1  *  80  V.  Static  value 

diode  is  doped  to  5  x  10* 5  of  VB1  is  indicated  by  dashed 

cm”3.  curve.  Parameter  is  dV/dt,  V/ns. 


Mn(V) 


show  that  the  maximum  voltage  attained  increases  with  the  rate  of  volt¬ 
age  rise.  The  difference  between  this  maximum  voltage  (vB2)  and  that 
obtained  in  static  calculations  (vB20)  is  plotted  in  figure  1%  on  a  log- 
log  scale  against  dV/dt.  The  slope  of  this  curve  is  2/3.  Two  of  the 
points,  those  for  dV/dt  -  2  and  4  x  10n  V/s,  are  slightly  above  the 
straight  line  because  the  temperature  of  the  diode  was  allowed  to  in¬ 
crease  in  time.  All  other  points  were  calculated  for  a  constant  temper¬ 
ature  of  300  K.  The  two  points  for  the  highest  rates  of  dV/dt  fall 
below  the  line  because  the  initial  voltage,  required  to  prevent  negative 
fields  from  developing  for  early  times  (small  t) ,  maintained  an  appre¬ 
ciable  current  minimum.  It  will  be  shown  later  that  the  presence  of  a 
residual  current  lowers  the  breakdown  voltage.  Also  plotted  in  figure 
15  as  a  function  of  dV/dt  is  the  difference  in  current  densities  at  VB2 

.  The  straight  line  was  obtained 
that  obtained  from  static  calcula¬ 
tions.  The  adjustment  was  less  than  the  probable  error  in  determining 
JB2Q.  The  current  differc-nce  points  show  more  scatter  than  the  voltage 
difference  points.  This  should  be  expected  from  the  greater  difficulty 
in  determining  the  JB2  values  compared  to  the  values  of  VB2.  The  slope 
of  the  AJ  curve  is  less  than  2/3,  and  other  similar  calculations  show  a 
scatter  in  slopes.  Data  foi  the  complementary  N  diode  fell  on  the  same 
line  as  for  the  P  diode  when  VB2Q  was  increased  by  0.5  V  from  that 
determined  by  the  static  calculations.  Again,  this  is  less  than  the 
error  range  expected  for  determination  of  VB2g .  Figure  16  shows  the 


Figure  15.  Log-log  plots  of  ex-  ev (V/») 

cess  voltage  and  current  as  * 

function  of  dV/dt  for  diode  of 

figure  13.  Excess  voltage  is  Figure  16.  Excess  voltage  plotted 

defined  as  maximum  voltage  aB  function  of  dV/dt  for  N-type 

attained  with  dV/dt  >  0,  less  the  diode.  Assumed  value  for  VB2  is 

static  value  with  dV/dt  =*  0.  labeled  for  each  curve. 

Excess  current  is  similarly 
defined. 


and  at 


—  VB20*  i,e*'  AJ»2  “JB2  ”B20 
with  a  small  adjustment  in  Jsj.  from 
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sensitivity  of  this  type  of  plot  to  the  choice  of  VB2£).  It  follows  that 
one  can  (usually)  obtain  VB20  and  Jjjjo  acre  accurately  (and  certainly 
much  more  cheaply)  by  a  series  of  dynamic  rather  than  static  computa¬ 
tions  • 

6.2  Constant  Overvoltage 

Calculations  have  also  been  made  for  the  same  P  dioae  at  300  K 
with  various  constant  overvoltages.  The  external  circuit  used  Rg  =  1  fi 
and  C  =  1  x  10-H  F.  The  initial  voltage  was  20  V.  A  log-log  plot  of 
breakdown  time  against  applied  overvoltage  is  shown  in  figure  17.  The 


time  of  breakdown  was  arbitrarily 
chosen  to  be  that  which  was  required 
to  attain  a  current  density  of  1  x 
10s  A/cm2.  Calculations  become  un¬ 
stable  and  are  halted  at  somewhat 
higher  (but  variable)  currents.  The 
overvoltage  as  plotted  is  adjusted 
for  the  voltage  drop  across  Rs  due 
to  the  breakdown  current.  The 
straight  line  in  figure  17  is  drawn 
with  a  slope  of  -2/3.  For  the  over¬ 
voltage  of  0.2  V,  the  breakdown  time 
was  also  calculated  with  the  diode 
temperature  increasing  due  to  the 
power  dissipated.  The  breakdown 
time  increased  about  25  percent. 
The  temperature  increase  at  the 
junction  is  only  5  K.  Longer  break¬ 
down  times  become  costly  to  calcu¬ 
late,  but  it  is  evident  that,  since 
V  increases  with  temperature,  the 
time  to  breakdown  must  increase  more 
and  more  rapidly  as  the  overvoltage 
decreases • 


Figure  17.  Log-log  plot  of 
time  to  breakdown  versus  over¬ 
voltage  for  P-type  diode.  All 
calculations  were  for  a  con¬ 
stant  300-K  temperature,  ex¬ 
cept  for  indicated  point, 
where  diode  heating  was  sim¬ 
ulated  . 


7.  TEMPERATURE  DEPENDENCE 

7.1  Isothermal  Calculations 

Isothermal  or  constant  temperature  calculations  have  also  been 
made  at  various  higher  temperatures  for  the  P-type  diode  considered  in 
section  6.  The  dynamic  characteristics  calculated  with  dV/dt  »  1  x  101 
V/s  at  various  temperatures  are  shown  in  figure  18.  It  is  seen  that 
both  VB1  and  VB2  increase  with  temperature.  These  calculations  were 
made  at  enough  different  rates  of  dV/dt  to  obtain  extrapolations  as 
shown  in  figure  15.  The  results  of  the  extrapolation  are  shown  in  table 
1.  Note  that  both  VB.  and  VB2«  increase  with  temperatures  for  400  and 


Figure  18.  Dynamic  characteris¬ 
tics  of  P-type  diode  of  figure  13 
at  different  temperatures,  as 
labeled.  Calculations  were  made 
with  dV/dt  -  1  x  10 10  V/s. 


TABLE  1.  BREAKDOWN  DATAa  FOR  4- Mm  DIODES  DOPED  TO  5  x  1015  ca~3, 

DIODE  AREA  -  1  x  1 0“4  cm2 


Type 

T 

(K) 

ft 

K* 

(ohms) 

(kV^.2) 

N 

300 

79 

2.5 

94.5 

23 

P 

300 

80 

6.6 

101.5 

16 

P 

400 

87 

9.4 

105.6 

9 

P 

500 

92 

18.4 

112 

6 

P 

600 

88 

17.6 

106 

12 

aData  from  fit  to  equation  (29),  except  600  K, 


500  K,  but  decrease  for  600  K.  The  increase  in  breakdown  voltages  with 
temperature  is  due  to  the  decrease  in  ionization  coefficients  with 
temperature  and,  to  a  very  small  degree,  the  decrease  in  saturation 
velocities  with  temperature.  The  decrease  in  VB2  at  600  K  is  due  to 
thermal  injection  of  minority  carriers  into  the  high  field  region  from 
the  high  doping  boundary  regions.  This  will  be  discussed  further  in  the 
next  section. 


A  similar  series  of  calculations  was  made  for  a  10-pm  N-type 
diode  doped  to  4  x  1015  cm-3.  An  area  of  5  x  10~5  cm2  was  assumed,  and 
the  calculations  were  made  with  Rg  -  50,  C  *  1  x  10“*2  F,  and  dV/dt  «  1 
x  1011  V/s.  Dynamic  characteristics  for  various  temperatures  are  shown 
in  figure  19.  The  PTF  of  these  diodes  is  about  0.65  and  VB2  is  nearly 


double  VB1*  Transit  time  oscillations  may  be  noted  in  figure  19.  At  a 
given  temperature,  the  frequency  of  oscillations  increases  with  current, 
and  at  a  given  current  the  frequency  decreases  with  temperature.  Both 
of  these  effects  are  observed  experimentally.  The  range  of  frequencies 
for  this  series  of  calculations  is  from  4  to  15  GHz.  Avalanche 
oscillations  will  be  discussed  further  in  section  14. 

The  rate  of  temperature  rise  in  a  semiconductor  is  proportional 
to  the  field  and  the  current  density  (see  eq  (6)).  Figure  20  shows  a 
plot  of  the  rate  of  temperature  rise  for  silicon  as  a  function  of 
current  density,  J,  and  with  the  electric  field,  E,  as  the  parameter. 
To  use  these  curves  one  needs  to  know  the  diode  area  in  order  to 
calculate  J  from  the  total  current,  and  the  diode  depletion  width  in 
order  to  calculate  E  from  the  voltage.  If  the  diode  area  and  width  are 
unknown,  they  may  be  estimated  from  measuring  the  breakdown  voltage  and 
the  capacitance  versus  voltage  curves. 

The  excess  second  breakdown  voltage,  AVQ2,  was  plotted  against 
dV/dt  for  a  P-type  diode  at  300  K  in  figure  15.  Similar  data  at  higher 
temperatures  are  plotted  in  figure  21 .  Also  shown  in  the  same  figure 
are  experimental  measurements  of  Yanai  et  al.21  Their  data  were 
obtained  in  a  study  of  the  initiation  of  avalanche  oscillations.  It  is 
seen  that  their  data  fit  just  below  the  400  K  calculated  curve.  Their 


characteristics  with  temperature  Figure  20.  Bate  of  temperature 
for  10-ym  diode  doped  to  N  *  4  x  rise  for  silicon  diodes  as  func- 
1015  cm"3.  Transit  time  oscil-  tion  of  current  density.  Param- 
lations  are  observed.  eter  is  electric  field. 

21W.  Yanai,  H.  Torizuka,  N.  Yamada,  and  K .  Ohkubo,  Experimental 
Analysis  for  the  Large- Amplitude  High-Efficiency  Mode  of  Oscillation 
with  Si  Avalanche  Diode,  IEEE  Trans.  Electron  Devices,  ED- 17  (December 
1970),  1067-1076. 


Figure  21.  Plot  of  overvoltage  as 
function  of  rate  of  rise  of  applied 
voltage.  Experimental  data  of  Yanai 
et  al  (ref  21)  are  shown  as  crosses 
(dashed  line)  and  computed  data  as 
circles.  Solid  lines  are  drawn  with  5 
slope  of  2/3.  Parameter  is  temper-  j ■ 
ature  (K). 


measured  slope  is  slightly  less,  as  would  be  expected  if  their  longer 
pulses  (low  dV/dt)  resulted  in  a  higher  temperature.  It  was  stated 
earlier  that  the  curves  for  complementary  (N  and  P)  diodes  fell  on  the 
same  curves.  Further  calculations  have  shown  that  these  curves  are 
essentially  independent  of  doping  density.  The  variation  with  temper¬ 
ature  is  predominantly  a  result  of  the  decrease  in  the  saturation  veloc¬ 
ities  of  the  holes  and  electrons  as  the  temperature  increases. 


7.2  Accuracy  of  Isothermal  Calculation 


The  question  arises  as  to  how  useful  are  the  isothermal  calcu¬ 
lations  when  the  temperature  across  the  diode  is  not  uniform.  Since  the 
total  current  across  a  diode  is  constant  due  to  current  continuity,  the 
temperature  rise  distribution  across  a  steady-state  diode  will  be  pro¬ 
portional  to  the  field  at  each  point.  Therefore,  neglecting  heat  flow, 
the  temperature  increase  at  any  time  will  also  be  proportional  to  the 
field  at  that  point:  the  temperature  distribution  has  the  same  shape  as 
the  field  distribution.  Calculations  were  made  at  300,  400,  and  500  K 
for  a  12-um  N-type  diode  doped  to  3  x  10 * 5  cm  Portions  of  the  static 
characteristics  at  400  and  500  K  are  shown  in  the  inset  of  figure  22. 
The  field-distribution  calculation  for  400  K  at  an  applied  voltage  of 
200  V  was  used  to  determine  a  temperature  profile  with  an  arbitrary 
maximum  of  600  K,  an  increase  of  300  K.  The  temperature  increase  at 
each  of  the  other  grid  points  was  the  same  fraction  of  300  K  as  the 
field  at  that  grid  point  was  of  the  maximum  field.  The  temperature 
profile  is  shown  in  figure  22  as  the  short  dashed  curve.  It  has  an 
average  temperature  of  440  K.  The  inset  in  figure  22  also  shows  a 
comparison  of  static  characteristics  for  the  given  temperature  profile 
(dashed  line  2)  and  the  isothermal  440-K  profile  (solid  curve  1).  The 
difference  in  current  for  the  two  profiles  is  quite  small  for  200  V 
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Figure  22.  Variation  of  field 
profiles  with  temperature  dis¬ 
tribution.  Inset  shows  por¬ 
tions  of  static  characteristics 
for  temperatures  of  400,  440, 
and  500  K  (solid  lines)  and  for 
temperature  distribution  shown 
in  main  figure  (dotted  line, 
right  scale).  Average  temper¬ 
ature  of  distribution  is  440 
K.  Field  distributions  1  and  2 
correspond  to  numbered  points 
in  inset. 

applied.  Fbr  220  V  applied,  the  difference  is  greater.  Also  shown  in 
figure  22  are  the  calculated  field  distributions  for  the  given  temper¬ 
ature  profile  (dashed  curve  2)  and  for  the  same  average  temperature 
(solid  curve  1).  The  higher  maximum  field  for  the  junction  at  600  K  is 
due  mainly  to  the  temperature  variation  of  the  ionization  coefficients, 
whereas  the  lower  field  in  the  low-temperature  region  results  from  the 
higher  mobility  at  the  lower  temperature.  Despite  these  rather  large 
field  variations,  it  is  seen  that  there  is  no  great  error  in  using 
constant  temperatures  instead  of  appropriate  temperature  profiles  in 
calculating  static  curves.  Farther  results  for  the  temperature  vari¬ 
ation  of  static  and  dynamic  characteristics  will  be  presented  in  later 
sections. 

7.3  Si li con -on-Sapphire  Diode 

The  silicon-on-sapphire  diode  used  by  Sunshine  and  Lampert22 
for  measurement  of  the  temperature  dependence  of  the  current-voltage 
static  characteristics  was  simulated  through  calculations.  Their  exper¬ 
imental  measurements  are  shown  in  figure  23.  Figure  23  also  shows  I-V 
characteristics  calculated  for  8  *  10* 6  cm”3 .  These  calculations  were 
made  with  dV/dt  =  1  x  1010  V/s.  Attempts  to  reach  the  steady  state 
below  the  current  "shoulders"  were  unsuccessful.  When  the  breakdown 
voltage  is  extrapolated  to  dV/dt  =  0,  N  =  8  x  io16  cm-3  gives  a  good  fit 
to  the  experimental  data.  An  area  of  5  x  10~8  cm2  was  required  to  fit 
the  current  shoulder  at  room  temperature.  A  diode  width  of  1  ym  is 
sufficient  to  fit  the  experimental  data  over  the  range  of  measured 
data.  A  wider  width  would  be  required  to  fit  the  measured  second- 
breakdown  voltages  (not  shown),  but  the  calculations  would  be  much  more 
costly  and  could  not  be  made  to  as  low  currents,  due  to  calculation 
instabilities. 

22R.  A.  Sunshine  and  M .  A .  Lampert,  Second  Breakdown  Phenomena  in 
Avalanching  Silicon-on-Sapphire  Diodes,  IEEE  Trans .  Electron  Devices, 
ED-19  (July  1972),  873-885. 
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Figure  23.  Measured  and  calcu¬ 
lated  current-voltage  character¬ 
istics  for  silicon-on-sapphire 
diode  with  temperature  (K)  as 
parameter.  Solid  curves  are 
experimental  measurements  of 
Sunshine  and  Lampert  (ref  22). 
Calculation  points  are  for  1-pm 
diode  doped  to  8  x  1016  cm-3; 

dV/dt  =  1  x  10 1°  V/s.  Dashed- 
line  extensions  are  guided  by 
other  calculations . 


Sunshine  and  Lampertz/  state  that  their  thin  film  material  had 
a  carrier  mobility  of  from  50  to  75  percent  of  bulk  material.  The  low- 
field  mobility  of  silicon  doped  to  about  1017  cm-3  is  only  about  50 
percent  of  undoped  material.  To  our  knowledge,  no  one  has  measured 
saturation  velocities  for  heavily  doped  material  at  any  temperature. 
Just  as  the  temperature  dependence  of  the  saturation  velocity  is  less 
than  that  of  the  low- field  mobility,  one  would  expect  less  dependence 
upon  the  doping  level  for  the  saturation  velocity  than  for  the  low-field 
mobility.  Moreover,  one  might  expect  less  variation  with  temperature  of 
the  saturation  velocity  of  highly  doped  material  than  is  measured  for 
intrinsic  material.  For  these  reasons,  an  option  to  use  mobilities 
dependent  upon  doping  level  was  not  used.  The  particle  velocities  are 
saturated  across  virtually  the  entire  depletion  region,  so  that  the  use 
of  the  variation  of  the  low-field  mobility  with  doping  level  would  be 
questionable.  However,  the  use  of  the  intrinsic  material  mobility 
overstates  the  calculated  current  density.  If  then  an  a  posteriori 
correction  is  made  to  account  for  the  lower  mobility  of  the  highly  doped 
thin  film  material,  the  area  used  to  fit  the  experimental  data  is  in  the 
range  of  those  listed  by  Sunshine  and  Lampert. 

lhe  measured  current  at  the  maximum  voltage  was  reported  to  be 
11.8  mA.  When  the  calculated  currents  for  were  extrapolated  to 

dV/dt  =  0,  a  current  range  of  from  10  to  15  mA  was  the  best  estimate. 
When  constant  applied  voltages  were  used  to  better  pinpoint  this 
current,  IMPATT  oscillations  of  approximately  60  GHz  were  calculated. 


22«.  A.  Sunshine  and  M.  A.  Lampert,  Second  Breakdown  Phenomena  in 
Avalanching  Silicon-ort-Sapphire  Diodes,  IEEE  Trans.  Electron  Devices, 
ED- 19  (July  1972),  873-88 5. 
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There  is  a  small  difference  between  the  measured  and  calculated 
variation  of  the  first  breakdown  voltage  with  temperature.  This  could 
be  due  to  the  omission  of  field  emission  in  the  calculations.  The 
maximum  fields  calculated  vary  from  6  x  10^  V/cm  at  300  K  to  8  x  10^ 
V/cm  at  650  K.  At  these  fields  some  field  emission  can  be  expected, 
since  field  emission  dominates  at  1  x  106  V/cm.  Field  emission 
increases  with  temperature,  whereas  ionization  coefficients  decrease 
with  temperature.  Experimentally,  it  has  been  observed  that  breakdown 
voltages  decrease  with  increasing  temperatures  for  diodes  whose  break¬ 
down  voltages  are  less  than  about  8  V.  There  is  qualitative  but  poor 
quantitative  agreement  in  the  variation  of  the  thermal  emission 
current.  In  this  case,  the  experimental  error  is  probably  greater.  The 
experimental  data  show  a  large  drop  in  the  maximum  voltage  obtainable  as 
the  temperature  exceeds  600  K.  No  decrease  is  calculated,  even  at  650 
K,  but  based  on  other  calculations,  one  should  expect  a  decrease  at  a 
slightly  higher  temperature^  If  one  examines  the  measured  curve  for  550 
K  carefully,  one  notes  an  "s"  shaped  curve  segment  in  the  current  at 
about  0.2  to  0.3  mA.  In  the  original  figure,  many  data  points  are 
involved.  One  could  easily  dismiss  this  as  experimental  error,  except 
one  notices  the  same  shape  in  the  600-  and  650-K  calculated  character¬ 
istics. 


Considering  the  above  discussion,  the  overall  fit  between  the 
measured  and  calculated  characteristics  is  considered  to  be  more  than 
satisfactory.  It  would  be  useful  to  similarly  compare  diodes  doped  at 
lower  levels,  but  no  similar  experimental  data  are  known  to  be 
avai lable . 


8.  THERMAL  INJECTION  CURRENTS 
8.1  Initial  Discussion 


In  section  4  it  was  pointed  out  that  the  use  of  low-field 
resistivity  versus  temperature  curves  at  high  fields  was  erroneous. 
Also,  in  the  literature  of  second  breakdown,  the  vague  term  "leakage 
current"  is  often  used  (or  misused).  What  is  probably  meant  is  the 
reverse  saturation  current  of  the  ideal  diode  theory.  The  term  leakage 
current  is  misused  when  it  is  considered  to  be  a  function  of  the 
junction  temperature. 

In  section  7  it  was  shown  that  breakdown  voltages  increased 
with  temperature  initially — primarily  because  of  the  decrease  of  ioni¬ 
zation  coefficients  with  temperature,  but  secondarily  because  of  the 
decrease  of  saturation  velocity  with  temperature — then  decreased  with 
temperature  because  of  thermal  injection.  The  process  of  thermal  injec¬ 
tion  will  now  be  discussed. 
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The  reverse  saturation  current  of  the  ideal  diode  equation  has 
two  components:  a  diffusion  component  and  a  thermal  generation  compo¬ 
nent.  The  ratio  of  these  components  is  a  function  of  temperature  and  of 
voltage.  Practically,  the  diffusion  component  is  most  important  in 
forward-biased  diodes  and  the  thermal  generation  in  reverse -biased 
diodes.  The  current  injected  from  the  low-field  boundary  region  into 
the  high-field  depletion  region  is,  in  this  paper,  termed  the  thermal 
injection  current.  The  term  leakage  current  is  rejected  because  it' can 
also  pertain  to  surface  leakage  current  and  because  it  has  been  so  often 
misused.  The  term  reverse  saturation  current  is  not  as  descriptive  in 
picturing  the  actual  physical  mechanism. 


The  thermal  injection  cur¬ 
rent  consists  of  two  components : 
electron  current,  AJnQ,  injected 
from  the  cathode  low- field  region 
and  hole  current,  AJpo,  injected 
from  the  anode  low-field  region; 
depending  upon  the  types  of  contact, 
injecting  or  ohmic,  one  may  dominate 
over  the  other.  For  most  of  our 
calculations  we  have  assumed  that 
Jno  =  Jpo  =  Jo-  The  dependence  of 
JQ  on  temperature  is  the  same  as 
that  of  n^,  given  in  equation 
(16).  The  experimentally  measured 
reverse  saturation  currents  are  in 
poor  agreement  with  theoretical  cal¬ 
culations,  so  an  arbitrary  choice 

has  been  made  for  300  K  that  J_  =  1 
—  n  o 

x  10  '  A/cm.  The  variation  of  JQ 

with  temperature  is  given  in  figure 

24. 


Figure  24.  Assumed  variation  of 
JQ  with  temperature.  Variation 
is  based  on  equation  (16). 


While  the  thermal  generation  of  carriers  is  important  in  the 
low-field  boundary  regions,  it  is  not  important  in  the  high-field  deple¬ 
tion  region,  because  of  two  phenomena.  First,  in  practically  all  the 
calculations  made  with  this  program,  the  np  product  of  the  carrier 
densities  has  greatly  exceeded  n£  once  V01  has  been  exceeded.  Then, 
according  to  equation  (10),  recombination  rather  than  generation  will 
occur.  However,  since  n?  depends  exponentially  on  temperature,  it  will 
exceed  the  np  product  at  high  enough  temperature.  Even  at  600  K,  where 
n^  =  4  x  io^°  cm-6,  typical  calculated  np  products  exceed  n?  well  below 
VB2.  Second,  even  when  n?  exceeds  np,  thermal  generation  will  be 
minimal,  since  the  transit  time  across  the  high-field  region  is  so 
short.  At  the  saturation  velocity,  this  time  is  about  1  x  10-11  s  per 
micrometer  of  distance  at  300  K.  This  time  is  orders  of  magnitude 
shorter  than  the  re combi nation -gene rati on  lifetime  of  silicon  (even  if 
gold  doped).  Further,  the  current  could  only  increase  by  a  factor  of  e 
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in  one  lifetime.  Therefore,  even  at  a  temperature  where  n|  >  pn,  little 
thermal  generation  could  be  expected  in  the  high-field  (temperature) 
regions.  It  is  the  temperature  of  the  low-field  injecting  regions  that 
is  important  in  reducing  the  voltage  required  for  current  multiplica¬ 
tion. 

Returning  to  the  discussion  of  the  (low-field)  resistivity 
versus  temperature  curves  (fig.  3),  it  is  the  reduction  of  mobility  with 
temperature  that  causes  the  increase  of  resistivity  with  T  at  the  lower 
temperatures,  whereas  the  decrease  of  resistivity  with  T  at  higher 
temperatures  is  caused  by  the  increase  of  n^  with  T.  Thus,  we  can 
expect  some  qualitative  agreement  from  thermal  models  that  use  the 
resistivity  versus  temperature  curves,  but  not  quantitative  agreement. 
The  main  reason  that  VB2  increases  initially  with  temperature  is  the 
decrease  of  the  ionization  coefficients  with  temperature,  not  that  of 
the  mobility  dependence.  The  decrease  of  VB_  at  still  higher  temper¬ 
ature  results  from  thermal  injection,  not  thermal  generation  at  the 
juncti on. 


Before  presenting  results  of  the  effect  of  thermal  injection 
currents,  it  should  be  pointed  out  that  JQ  may  also  be  used  to  simulate 
photoinjection  of  carriers,  to  simulate  the  injection  of  emitter  current 
into  the  base-collector  region  of  a  transistor,  and/or  to  simulate  the 
effect  of  ionizing  radiation  upon  second  breakdown.  The  last  effect 
will  be  discussed  in  section  12. 

8.2  Calculations  of  Thermal  Injection  Currents 

In  order  to  better  understand  the  influence  of  temperature  on 
diode  operation,  the  bulk  temperature  effects  (ionization  coefficients, 
saturation  velocities,  and  recombination  or  generation)  have  been 
separated  from  the  boundary  conditions  (thermal  injection).  Again,  in 
the  study  of  the  thermal  injection  currents,  it  is  useful  to  consider 
separately  diodes  which  punch  through  before  breakdown  and  those  not 
punched  through  at  breakdown.  Figure  25(a)  shows  characteristics  calcu¬ 
lated  for  a  4-un  P-type  diode  doped  to  5  x  1014  cm-3.  fliis  diode  is 
punched  through  well  below  breakdown  and  has  no  space-charge-limited 
(SCL)  current  region.  The  calculations  are  for  300  K  and  are  made  with 
dV/dt  =  2  x  io 10  V/s.  VB2  is  seen  to  be  reduced  as  JQ  is  increased,  but 
the  current  at  second  breakdown  increases.  Similar  calculations  were 
made  for  the  complementary  N-type  diodes.  The  results  are  shown  in 
figure  25(b).  Also  shown  in  the  figure  are  calculations  made  with  Jno 
unequal  to  J^.  It  is  seen  that  Jno*  the  electron  injection  current 
density,  is  more  effective  than  the  hole  injection  current.  This 
results  from  the  higher  ionization  coefficient  for  electrons  than  for 
holes  in  silicon.  When  the  multiplication  factor,  that  is,  M  «  J/J Q,  is 
plotted  against  voltage  for  these  punched -through  diodes,  the  curves  all 
merge  into  a  single  curve,  with  deviations  occurring  at  low  M  values  due 
to  initial  conditions,  and  at  high  M  values  due  to  space-charge  distor¬ 
tion  and  double  avalanching. 
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Figure  25.  Calculated  dynamic  characteristics  for  (a)  4-pm  P-type 
diodes,  doped  to  5  x  10 14  cm.  Parameter  is  JQ  in  A/cm2.  Arrows 
indicate  voltage  maxima.  Calculated  dynamic  characteristics  of  (b)  N- 
type  diodes,  doped  to  indicate  voltage  maxima.  Dashed  curve  A  is 
for  Jno  =  102  and  Jp0  =  1 0 1 ;  curve  B  is  for  Jno  =  10 1  and  Jpo  =  102. 
Calculations  were  made  with  dV/dt  =  2  x  10 10  V/s. 


Calculations  were  also  made  with  jQ  =  1  x  io~7  A/cm2,  corre¬ 
sponding  to  room  temperature,  but  with  constant  elevated  temperatures 
for  the  bulk  of  the  P-diode.  The  results  are  shown  in  figure  26  as  the 
dashed  curves.  Since  this  condition  assumes  no  heat  flow  to  the  low- 
field  regions,  this  will  be  termed  the  adiabatic  case.  Also  shown  in 
the  same  figure  are  characteristics  where  both  the  bulk  and  thermal 
injection  region  are  at  the  same  temperature — the  isothermal  case.  Hie 
maximum  voltages,  VB2«  for  each  case  are  plotted  as  a  function  of  tem¬ 
perature  in  figure  27.  When  the  thermal  injection  region  is  kept  at 
room  temperature,  VB2  1  ncreases  more  and  more  rapidly  at  higher  temper¬ 
atures,  but, with  increasing  thermal  injection,  a  maximum  is  reached  at 
about  500  K.  The  rate  of  decrease  in  VB2  with  temperature  becomes 
continuously  greater  above  the  temperature  for  maximum  voltage. 
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Figure  26.  Characteristics  for 
diode  of  figure  25(a)  with  bulk 
temperature  (K)  as  parameter. 
Fbr  300  K  and  dashed  curves ,  JQ  * 
1  x  10“7  A/cm2.  Fbr  solid  curves 
JQ  is  determined  by  temperature 
according  to  figure  24.  Arrows 
indicate  current  density  for  max¬ 
imum  voltage. 


TEMPERATURE  (K) 


Figure  27.  Second  breakdown 
voltage  as  function  of  temper¬ 
ature  for  diode  of  figure  26. 
Thermal  injection  current  is 
determined  by  bulk  temperature 
(isothermal  calculation)  for 
solid  curve,  whereas  thermal 
current  corresponds  to  300  K 
(adiabatic  calculation)  for 
dashed  curve. 


While  there  is  some  lowering  of  VB2  with  increasing  JQ,  even  at 
the  lo\<est  value  of  JQ  used  for  the  punched-through  diodes  discussed 
thus  far,  this  is  not  true  for  diodes  not  punched  through  at 
breakdown.  Calculated  characteristics  as  a  function  of  JQ  for  a  p-type 
diode  of  the  same  width,  but  doped  to  5  *  1015  cm-3  (a  factor  of  10 
larger),  are  shown  in  figure  28.  This  diode  is  not  punched  through  at 
breakdown  and  a  regime  of  SCL  current  is  noted.  The  SCL  regime,  where 
the  I-v  characteristic  is  linear,  shows  up  on  the  semilog  plot  as  a 
curved  region.  For  1  x  10"7  <  JQ  <  1  x  to"2,  no  difference  is  noted  in 
the  calculated  characteristics.  At  higher  JQ  values,  the  decrease  in 
Vgj  i»  again  observed.  A  comparison  of  the  decrease  in  VB2  a®  * 
function  of  JQ  for  each  of  the  diodes  is  shown  in  figure  29.  Fbr  high 
values  of  JQ,  the  lowering  of  V02  is  comparable. 
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Figure  28.  Calculated  dynamic 
characteristics  for  4-inn  P- 
type  diode,  doped  to  5  x  10 15 
cm"^,  as  a  function  of  JQ.  No 
lowering  of  VB2  is  noted  for 
this  punched-through  diode  as 
J  is  increased  from  10-7  to 
10~2  A/cm2. 


VOLTAGE  (V) 


Figure  29.  Plots  of 
VB2  for  4-ym  diodes  as 
function  of  JQ.  Dop¬ 
ing  density  (cm"3)  is 
indicated  for  each 
curve.  Points  A  and  B 
are  defined  in  figure 
25. 


10-*  10-1 
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Characteristics  for  the 
higher  doped  diode  were  also  cal¬ 
culated  for  both  the  adiabatic  and 
isothermal  cases  at  higher  temper¬ 
atures*  The  calculations  were 
made  at  dV/dt  =  1  x  10 10  V/s  and 
also  at  higher  rates  to  allow 
extrapolation  to  dV/dt  —  0.  The 
results  are  shown  in  figure  30. 
The  isothermal  and  adiabatic 
curves  are  the  same  up  to  nearly 
500  K;  500  K  corresponds  to  a  JQ 
value  of  20  A/cm2.  Only  the  iso¬ 
thermal  data  were  extrapolated  to 
dV/dt  =  0. 

The  calculations  of  this 
section  show  the  importance  of  the 
thermally  injected  current,  JQ. 
The  implications  of  this  effect  in 
comparing  differing  doping  pro¬ 
files  will  be  discussed  in  section 
10. 


9.  VARIATION  OF  DIODE  WIDTHS  AND 
DOPING  DENSITIES 


9.1  Varying  Doping  Densities 


Figure  30.  Second  breakdown  volt¬ 
ages  for  diode  of  figure  28  as 
function  of  temperature.  Dashed 
curves  are  for  dV/dt  =1  x  1010 
V/s.  Crosses  indicate  isothermal 
calculation  and  open  points 
indicate  adiabatic  calculations. 
Solid  curve  (pluses)  indicates  VB2 
extrapolated  to  dV/dt  =  0.  ( See 

fig.  16.) 


In  section  5  it  was  shown  that  for  a  constant  diode  width,  the 
V01  and  Rgc  varied  as  the  doping  level  was  changed  by  a  factor  of  2  but 
V02  was  essentially  unchanged.  These  observations  are  valid  over  a  much 
wider  range  of  doping  levels.  Figure  31  shows  the  family  of  character¬ 
istics  for  a  4-um  P-type  diode  at  300  K,  with  doping  ranging  from  1  x 
1015  to  2  x  1016  cm.  Again,  these  characteristics  were  calculated  with 
dv/dt  =1  x  10 10  V/s  to  reduce  costs  and  printout  paper.  As  may  be  seen 
from  the  values  of  VB1  and  VB2,  the  dynamic  characteristic  agrees  quite 
well  with  the  data  given  in  the  previous  section  for  static  calcu¬ 
lations.  Note  in  figure  31  that  while  VB1  varies  from  about  45  to  105 
V,  VB2  varies  only  by  about  3  V.  The  arrows  in  figure  31  indicate  the 
current  density,  JB2,  for  the  maximum  voltage.  Note  also  that  JB2 
increases  with  the  doping  level.  This  is  expected,  since  it  was  shown 
in  the  previous  section  that  double  avalanching  occurs  when  the  mobile 
carriers  exceed  the  doping  level.  The  punch-through  factor  (PTF)  varies 
from  about  0.4  for  the  doping  level  of  2  x  10 16  cm”3  to  about  7  for  the 
doping  of  1  x  1015  cm”3.  The  use  of  approximate  numbers  may  be 
explained  by  the  consideration  of  figure  32.  The  space-charge  free 
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Figure  31.  Characteristics  of  4- 
lim  P-type  diode  with  doping  level 
in  cm-3  as  parameter.  Calcu¬ 
lations  were  made  with  dV/dt  =  1  x 
1010  V/s.  Alternating  solid  and 
dashed  lines  are  used  only  for 
clarity. 


Figure  32.  Distributions  of  elec¬ 
tric  field  across  diode  of  figure 
31.  Currents  are  all  in  SCF  cur¬ 
rent  regime. 


(SCF)  field  distribution  across  each  diode  of  figure  31  is  plotted  in 
figure  32.  These  fields  are  determined  by  the  doping  profiles 
simulated,  since  space  charge  from  mobile  carriers  is  negligible.  While 
the  N+P  junction  is  exactly  defined,  the  gradual  increase  in  doping  at 
the  PP+  junction  leads  to  an  uncertainty  in  the  width  of  the  depletable 
region  and  thus  to  an  ambiguity  for  the  PTF. 


The  variation  of  the  field  distribution  with  current  density  is 
shown  for  the  diode  with  a  doping  level  of  2  x  101®  cm”3  in  figure  33. 
Note  the  apparent  similarity  of  the  curves  in  figures  32  and  33.  This 
results  from  the  compensation  by  the  mobile  holes  for  more  and  more  of 
the  fixed  doping  charge  as  the  current  increases.  As  noted  earlier,  the 
number  density  of  holes  rises  exponentially  in  the  high-field  junction 
region  and  is  relatively  constant  in  the  low-field  region.  At  JB_, 
field  distribution  is  quite  similar  for  all  doping  levels  as  might  be 
expected  from  the  similar  VB~  values. 


Figure  33.  Field  distributions 
across  P  =  2  x  1016  cm"3  diode 
for  various  current  densities 
(A/cm2).  Calculations  were  made 
with  dV/dt  =  1  x  1010  V/s. 


9.2  Varying  Widths 

.  A  second  family  of  diodes  of 
interest  is  that  with  the  same  doping 
density,  but  with  varying  widths. 
Figure  34  shows  300  K  characteristics 
for  N-type  diodes  doped  to  3  x  101S 
cm"3  and  with  width  varying  from  4  to 
12  pm.  Again  these  are  dynamic 
characteristics,  calculated  with 


dV/dt  =  1  x  10l1  V/sj  horizontal  arrows  indicate  current  for  voltage 
maxima.  Also  shown  by  crosses  are  the  values  of  V02  and  J02  obtained  by 
making  calculations  at  higher  dV/dt  values  and  extrapolating  to  dV/dt  * 
0.  For  the  two  diodes  of  extreme  widths,  meaningful  extrapolations 
could  not  (easily)  be  made.  The  current  extrapolation  was  uncertain  for 
the  10-pm  diode.  The  dashed  curves  at  low  J  for  the  10-  and  12-pm 
diodes  are  extrapolations  of  the  computer  data.  Since  the  PTF  is  less 
than  1  for  these  diodes,  VB1  is  the  same  as  for  the  8-pm  diode.  The 
field  distribution  at  low  J  values  for  each  diode  is  shown  in  figure 
35.  The  diode  widths  are  not  labeled,  but  can  be  identified  by  the 
point  where  E  approaches  zero.  Also  shown  is  the  field  distribution  for 
the  12-pm  diode  at  a  current  in  the  double  avalanche  region. 


Figure  34.  Characteristics  of 
N-type  diodes  doped  to  3  x  10* 3 
cm"3.  Diode  width  is  given 
in  pm  for  each  curve.  Calcu¬ 
lations  are  for  dV/dt  =»  1  x 
1010  V/s.  Horizontal  arrows 
indicate  current  density  for 
maximum  voltage.  Crosses  indi¬ 
cate  extrapolated  values  of  VBj 
and  JB2*  Vertical  double  arrow 
indicates  range  of  uncertainty 
in  current  density  for  10-pm 
diode. 
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9.3  Whole  Families 


Figure  35.  Field  distributions 
across  3  x  to1 5  cm'3  diodes  of 
various  widths  of  figure  34  at 
low  current  densities.  Also 
shown  is  field  across  12-um 
diode  in  double-avalanche  cur¬ 
rent  regime. 


One  needs  a  three-dimensional  plot  to  present  breakdown 
voltages  as  a  function  of  both  doping  level  and  diode  width  for  the 
entire  family  of  diodes  with  constant  doping  levels.  However,  the  same 
data  may  be  presented  on  two-dimensional  plots,  with  the  other  variable 
as  a  parameter.  Sze16  shows  a  log-log  plot  of  VB1  versus  doping  density 
for  one-sided  abrupt  P-N  junctions.  His  curves  were  calculated  with  the 
use  of  the  ionization  coefficients  of  Lee  et  al23  and  are  reproduced  in 
figure  36.  Also  shown  in  figure  36  are  the  VB,  and  VB2  data  of  this 
work.  The  data  are  for  either  P-  or  N-type  diodes;  the  small  difference 
does  not  show  on  this  voltage  scale.  Dashed  lines  connect  the  family  of 
diodes  with  the  same  width.  The  agreement  of  the  present  data  with 
Sze's  data  seems  satisfactory.  The  difference  is  that  which  is  expected 
from  the  different  ionization  coefficients  used.  In  figure  36,  all  non 
punched-th  ough  diodes  fall  along  a  common  diagonal  line  independent  of 
width.  As  the  diode  width  decreases,  punch-through  occurs  at  larger 
doping  densities.  V01  is  essentially  constant,  independent  of  doping 
level,  but  a  function  of  diode  width,  for  diodes  that  are  well  punched 
through.  The  short  curved  segments  between  the  straight  lines  corre¬ 
spond  to  PTF  21. 

Sze16  states  that  for  non-punched-through  diodes  of  Si,  Ge, 
GaAs ,  and  GaP,  V01  can  be  approximated  by 

VB1  -  60(Eg/1 .1  (n/1016)“3/4  .  (31) 

For  silicon,  Eg  *  1.12  eV,  so  that  for  silicon 

VB1  -  62(l016/N)3/4  .  (32) 

1  %.  M.  , Sze,  Physics  of  Semiconductor  Devices,  Wiley-Interscience, 
New  York  (1969),  39-41,  57-59,  121-126. 

2fc.  A .  Lee,  R.  A.  Logan,  R.  L.  Batdorf ,  J.  J.  Kleimack,  and  W. 
Wiegmann,  Ionization  Rates  of  Holes  and  Electrons  in  Silicon,  Phys. 
Rev.,  A134  (May  1964),  A761-773. 
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Our  data  fit  closely  to 


VB1  =  65(1016/n)2/3  . 


Grove20  shows  that  the  measured  VB1  data  for  silicon  have  a  slope  less 
than  three -fourths .  The  fit  of  that  experimental  data24  to  equation 
(33)  is  shown  in  figure  37. 

The  merged  curve  of  figure  36  may  be  separated  by  plotting  V01 
against  diode  width,  as  shown  in  figure  38.  Now  the  data  for  the 
punched-through  diodes  are  merged  into  the  single  curve  applicable  for 
intrinsic  material.  The  equation  for  intrinsic  material  is 


as  compared  to 


VB1  =  36d0*81 


VB1  =  33d0*92  , 


for  Sze's  data.  One  cautionary  word  is  necessary  here.  Hie  width  of 
the  entire  high-field  region  has  been  used  for  d  in  equation  (34), 
whereas  the  usual  practice  is  to  use  the  width  of  the  low  doping  region. 


ion  io«  10’*  10V 

DENSITY  (cm-5) 

Figure  36.  First  (VB1)  and 
second  (V^)  breakdown  volt¬ 
ages  as  function  of  doping 
density.  Diode  width  in  ym  is 
given  as  parameter.  Solid 
curves  are  VBj  curves  taken 
from  Sze  (ref  16).  Dashed 
curves  connect  calculated 
points  of  same  diode  width. 
Arrows  at  left  indicate  calcu¬ 
lations  for  intrinsic  diodes. 


Figure  37.  Comparison  of  measured 
and  calculated  VB1  variation  with 
doping  density.  Points  are  exper¬ 
imental  measurements  of  Miller  (ref 
24).  Dashed  line  is  equation  (33) 
from  Sze  (ref  16)  and  solid  line  is 
equation  (33)  from  our  calculations. 


20A.  S.  Grove,  Physics  and  Technology  of  Semiconductor  Devices,  John 
Wiley  and  Sons,  New  York  (1967),  44-52,  194-201, 

2  S.  L,  Miller,  Ionization  Rates  for  Holes  and  Electrons  in  Silicon, 
Phys.  Rev,,  105  (February  1957),  1246-1249 . 
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A  third  way  to  plot  breakdown  voltages  is  to  use  vB1  as  a 
parameter  on  a  log-log  plot  of  d  as  ordinate  and  N  as  abscissa.  Hiis 
plot  eliminates  the  merged  curves,  as  shown  in  figure  39.  It  gives  an 
overall  picture  of  the  breakdown  voltages  of  all  possible  combinations 
of  N  and  d.  A  great  majority  of  diodes  fall  in  or  near  the  curved 
regions  where  the  PTF  is  near  unity. 

For  the  abrupt  one-sided  junction  with  constant  doping  dE/dx  = 
eN/tce0,  equation  (5),  when  mobile  carriers  can  be  neglected.  Thus, 
below  breakdown, 


V  =  eNd2/2  K£0 


(36) 


This  is  the  basic  equation  to  determine  the  doping  profile  from  capa¬ 
citance  measurements  as  a  function  of  reverse  bias.  Knowing  the  area, 
the  capacitance  determines  d,  allowing  N  to  be  determined.  Plots  of  V 
for  silicon  as  a  function  of  N  with  d  as  parameter  are  shown  in  figure 
40.  The  intersection  of  these  curves  with  the  V01  curve  gives  the  PTF  = 
1  position.  Equating  equation  (36)  with  equation  (33),  we  find 


V  =  C'N“2/3  =  C"Nd2 


dB  =  C"' 


moot  WIDTH  w«l 


Figure  38.  Same  breakdown  data 
plotted  against  diode  width  and 
with  doping  density  (cm”3)  as 
parameter.  Solid  curves  are 
from  Sze  (ref  16)  and  points 
are  results  from  our  calcula¬ 
tions. 


Figure  39.  First  breakdown 
voltages  plotted  as  parameter 
for  variations  in  both  diode 
width  and  doping  density. 
Diagonal  line  is  locus  of 
points  for  PTF  »  1 . 


51 


This  shows  that  the  depletion  region  at  breakdown  varies  (inversely) 
slightly  less  than  linearly  with  doping.  Equation  (37)  is  again  the 
condition  for  PTF  =  1  and  is  also  shown  in  figure  39.  Since  for 
constant  doping, 

Em,x  =  2V/d, 
max 

we  find,  by  using  equation  (37)  for  the  breakdown  condition,  that 

Emax  =  C»»Nl/6  .  (38) 

This  shows  that  the  junction  field  at  breakdown  is  not  constant  as  often 
assumed,  but  varies  slightly  with  the  doping  density.  A  plot  of  calcu¬ 
lated  E  at  breakdown  as  a  function  of  N  is  shown  in  figure  41. 
max 


Figure  40.  Calculated 


B1 


function  of  doping  density  for 
various  diode  widths,  from  equa¬ 
tions  (33)  and  (34).  Dashed 
lines  give  depletion  width  for 
voltages  below  breakdown.  Each 
intersection  with  solid  diagonal 
lines  indicates  that  PTF  =  1  for 
that  doping. 


Figure  41.  Maximum  field  at  break¬ 
down  as  function  of  doping  dens¬ 
ity.  Curve  is  plot  of  equation 
(38);  points  are  for  diffused  one¬ 
sided  junctions. 


10.  DIFFERING  DOPING  PROFILES 

The  doping  profiles  considered  in  the  previous  sections  were  one¬ 
sided  diffuse  junction  diodes.  Most  approximations  of  semiconductor 
breakdown  have  assumed  abrupt  junction  one-sided  diodes.  This  approx¬ 
imation  is  often  valid  and  does  make  junction  positions  and,  therefore, 
depletion  widths,  PTF's,  etc,  more  definite.  It  also  facilitates  com¬ 
parison  with  other  calculations. 


10.1  One-Sided  Diodes 


Calculated  static  current-density /voltage  characteristics  for 
abrupt  junction  one-sided  diodes  at  300  K  are  shown  in  figure  42.  The 
doping  levels  are  1  x  i o* 7  cm”3  for  diode  widths  of  0.4,  0.5,  and  0.6  u m 
and  2  x  ip1 7  cm-3  for  a  0.3-pm  diode.  The  diodes  are  P-type  except  for 
one  N-type  diode.  Carrier  densities  across  the  N-type  diode  are  shown 
in  figure  43.  The  electric  field  distributions  across  the  same  diode  at 
the  same  current  densities  are  shown  in  figure  44.  Similar  field  dis¬ 
tributions  across  the  complementary  (same  width  and  doping  level)  P-type 
diode  are  shown  in  figure  45. 


For  the  lowest  voltages  shown  in  both  figures  44  and  45,  the 
field  is  linear  from  the  junction  to  the  edge  of  the  high  doping 
region.  Since  these  voltages  are  essentially  first  breakdown  voltages 
(VBl)  for  these  diodes,  the  PTF  =  1.  The  0. 6-vim  diode  does  not  punch 
through  at  breakdown  (PTF  =  0.5/0. 6  =  0.83),  so  its  Vg1  is  equal  to  that 
of  the  0.5-vim  diode,  as  may  be  seen  in  figure  42.  The  0.4-pm  diode 
punches  through  before  breakdown  (PTF  =  0.5/0. 4  =  1.25);  hence,  its  Vg^ 
is  less  than  that  for  the  wider  diode.  The  transition  from  the 
quadratic  current  dependence  upon  voltage  for  the  0.4-um  diode  to  the 
linear  dependence  for  the  0.6~im  diode  may  also  be  seen  in  figure  42. 


Figure  42.  Calculated  charac¬ 
teristics  for  abrupt  one-sided 
diodes.  One  curve  is  for  dop¬ 
ing  of  2  x  10 17  cm”3;  other 

curves  are  for  1  x  10 17  cm”3. 
Diode  width  is  given  for  each 
curve.  Dashed  curve  is  for  ll- 
type  device;  other  curves  are 
for  P-type  devices. 


Figure  43.  Carrier  densities 
for  N-type  diode  of  figure 
42.  Doping  profile  is  given 
by  short-dashed  curve; 
electron  density  is  given  by 
solid  curves;  hole  density  is 
given  by  long-dashed  curves. 
Parameter  is  current  density. 
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Figure  44.  Field  distri¬ 
butions  corresponding  to 
carrier  profiles  of  figure 
43.  Diode  voltage  is 
curve  parameter. 


DISTANCE  (fdn) 


Figure  45.  Field  distri¬ 
butions  for  P-type  diode 
which  is  complementary  to 
diode  of  figure  44. 
Again,  diode  voltages  are 
given . 


As  shown  in  figure  42,  the  P-type  diode  has  a  higher  VB1  than 

does  the  complementary  N-type  diode.  The  p-type  diode  also  has  a  higher 

space-charge  resistance  than  does  the  N-type  diode,  as  noted  by  the 
voltage  difference  increasing  with  current. 

At  the  intermediate  current  level  shown  in  figure  43,  the 
space  charge  of  the  mobile  carriers  becomes  an  appreciable  fraction  of 
the  doping  level,  and  the  field  profile  is  no  longer  linear  (fig.  44). 
At  the  highest  current  level  shown  in  figure  43,  the  electron  density 
exceeds  the  sum  of  the  doping  level  and  the  hole  density.  This  results 
in  a  field  minimum  (saddle)  as  shown  in  figure  44. 

The  peak  electric  field  at  the  junction  is  essentially 
constant,  independent  of  voltage  (above  VB1 )  for  the  N-diode,  as  shown 

in  figure  44.  The  field  maximum  does,  however,  move  slightly  from  the 

junction  toward  the  P+  doping  region.  The  peak  field  for  the  p-diode 
decreases  with  voltage  as  the  voltage  increases  above  VB1 ,  as  shown  in 
figure  45.  At  low  current,  the  maximum  field  for  the  P-diode  slightly 
exceeds  that  for  the  N-diode,  but  at  higher  currents  it  drops 
appreciably  below  the  constant  field  maximum  of  the  N-diode.  The  total 
voltage  difference  between  the  two  diodes  at  the  same  current  increases 
because  the  space-charge  field  in  the  P-diode  increases  with  current 
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faster  than  it  does  in  the  N-diode.  The  saddle  appears  at  about  a  20 
percent  larger  field  for  the  p-device.  These  differences  arise  from  the 
higher  ionization  rate  for  electrons  than  for  holes  in  silicon. 

The  results  reported  in  this  section  for  the  abrupt  one-sided 
diodes  differ  minimally  from  those  reported  in  the  previous  sections  for 
diffuse-junction  one-sided  diodes.  The  maximum  field  is  more  sharply 
pointed  for  the  abrupt  doping  profile. 

10.2  Two-Sided  Diodes 


Static  current-density/voltage  characteristics  for  abrupt 
two-sided  diodes  doped  to  1  x  1017  cm-3  are  shown  in  figure  46.  A 
logarithmic  current-density  scale  is  used  to  cover  a  wide  current 
range.  These  characteristics  were  calculated  point  by  point,  since 


errors  resulted  even  at  low  dV/dt 
rates.  The  PTF  equals  1  for  the  0.8-pm 
diode.  The  small  ellipses  at  the 
highest  current  densities  indicate 
IMPATT  oscillations25  for  three  of  the 
diodes.  Avalanche  oscillations  will  be 
discussed  in  section  14.  The  doping 
profile  and  the  mobile  carrier  dens¬ 
ities  are  shown  in  figure  47  for  the 
0.8-pm  diode.  The  corresponding  field 
distribution  at  each  current  density  is 
shown  in  figure  48.  For  the  symmetric 
two-sided  diode  doped  to  1  x  10 17  cm”3, 
PTF  =  1  for  a  0.8-pm  width.  All  diodes 
wider  than  0.8-pm  have  equal  Vg^,  but 
narrower  diodes  have  lower  Vg^.  The 
field  distribution  ?t  breakdown  has  the 
form  of  an  isosceles  triangle.  At 
higher  currents,  the  effect  of  the  mo¬ 
bile  space  charge  is  to  lower  the  peak 
field  at  the  junction  as  it  is  raised 
near  each  diode  edge  (fig.  48).  The 
effect  is  to  broaden  the  region  wherein 
current  multiplication  occurs,  as  may 
be  seen  in  figure  47.  The  higher 
fields  in  the  P  region  than  in  the  N 
region  again  result  from  the  difference 
in  the  ionization  rates  for  electrons 
and  holes. 


Figure  46.  Characteristics  of 
symmetric  two-sided  diodes 
doped  to  1  x  1017  cm”3  for  four 
diode  widths.  Ellipses  at 
highest  currents  are  equilib¬ 
rium  IMPATT  oscillations. 
Voltage  for  0.8-pm  diode  at 
J  ■  6  x  101*  A/ cm2  is  high 
because  too  low  doping  used  in 
N+  and  P+  regions  for  this 
current • 


2SA.  i,  ward.  Modes  of  Avalanche  Oscillations  in  Silicon  Diodes,  IEEE 
Trans.  Electron  Devices,  ED-25  (June  1978),  683-687. 
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U  0.2  0.4  0.0 

DISTANCE  (Mm) 

Figure  47.  Carrier  densi¬ 
ties  for  symmetric  two- 
sided  diode  doped  to  1  x 
Doping  profile 
as  dotted  curve, 
electron  densities  as  solid 
curves,  and  hole  densities 
as  dashed  curves.  Param¬ 
eter  is  current  density. 


1017  cm-3, 
is  shown 


10' l 
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Figure  48.  Field  distri¬ 
bution  corresponding  to 
densities  of  figure  47. 
Diode  voltage  is  given  for 
each  distribution. 


10.3  Linearly  Graded  Diodes 

Static  current-density/voltage 
characteristics  for  linearly  graded 
doping  profiles  are  shown  in  figure 
49.  Again,  the  current-density  scale  is 
logarithmic  to  cover  the  range  of  inter¬ 


est.  Three  doping  gradients  are  illustrated:  1  x  1020,  1  x  1021,  and  1 
x  1022  cm-4.  Diode  widths  are  4,  5,  and  8  vim;  1.6  and  2.0  pmj  and  0.5 
and  1.0  pm,  respectively.  Higher  current  densities  are  required  for 
space-charge  effects  to  appear  for  the  higher  doping  gradients.  Recall 
that  the  currents  increase  exponentially  in  the  SCF  current  regime,  but 
increase  less  rapidly  with  voltage  due  to  the  mobile  carrier's  space- 
charge  resistance.  The  doping  profile  and  the  carrier  densities  across 


the  5-pm  diode  with  a  doping  gradient  of  1  x  10 


20 


cm 


-4 


are  shown  in  fig¬ 


ure  50.  Both  N  and  P  doping  are  shown  as  positive  in  all  figures;  there 
is  no  change  in  gradient  at  the  metallurgical  junction.  The  correspond¬ 
ing  field  profiles  are  shown  in  figure  51.  Fbr  low  currents,  the  field 
is  parabolic;  with  increased  current  density,  the  field  maximum  de¬ 
creases  and  moves  into  the  P+  region.  The  PTF  equals  1  for  a  4-pm  diode 
with  a  gradient  of  1  x  102°  cm-4.  For  the  5-pm  diode  shown  in  figure 
51,  the  depletion  (high-field)  region  widens  at  higher  currents.  At  the 
highest  current  shown,  the  field  distribution  closely  resembles  that  for 
a  lightly  doped  (PTF  >  1),  N-type,  broadly  diffused  junction  diode.  In 
fact,  at  even  higher  current  densities,  field  distributions  tend  to 
become  independent  of  doping  profiles. 
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CURRENT  DENSITY  (A/em*) 


Figure  49.  Characteristics  of 
three  families  of  linearly 
graded  diodes,  with  gradients 
of  1  x  1020,  1  x  1021,  and  1  x 
1022  cm”1*.  Diode  width  is 
parameter  for  each  curve. 


Figure  51.  Field  distri¬ 
butions  corresponding  to 
current  densities  of  fig¬ 
ure  50.  Diode  voltage  is 
given  for  each  distribu¬ 
tion. 


Figure  50.  Carrier  densi¬ 
ties  for  5-pm  diode  with 
linear  doping  gradient  of 
1  x  1020  cm-4.  Doping 
densities — straight  dotted 
lines — are  all  plotted  as 
positive  numbers.  Elec¬ 
tron  densities  are  given 
by  solid  curves  and  hole 
densities  by  dashed 
curves.  Current  densities 
are  given. 


10.4  Other  Doping  Profiles  and 
Discussion 

With  the  advent  of  the 
use  of  ion  implantation  to  form  P-N 
junctions,  other  doping  profiles  are 
feasible.  In  particular,  high- low 
(Hi-Lo)  and  Lo-Hi-Lo  profiles,  both 
single-  and  two-sided,  have  been 
used  with  advantage  in  IMPATT  di¬ 
odes.  It  is  planned  to  make  calcu¬ 
lations  simulating  these  profiles, 
since  they  should  have  superior 
thermal  properties. 


* 
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It  was  shown  in  section  8  that  thermal  injection  currents 
play  a  major  role  in  second  breakdown  for  times  greater  than  the  thermal 
time  constant.  For  short  times,  essentially  adiabatic  conditions,  the 
temperature  distribution  is  qualitatively  the  same  as  the  field  distrib¬ 
ution.  Therefore,  the  field  distributions  calculated  in  this  section 
can  be  used  to  study  the  effect  of  doping  profiles  upon  temperature 
distributions.  Two-sided  diodes  have  their  highest  temperature  farther 
from  the  low-field  thermal  injection  regions  than  do  one-sided  diodes. 
Abrupt  two-sided  diodes  have  smaller  temperature  gradients  at  the  low- 
field  boundary  than  do  linearly  graded  diodes.  Thus,  the  former  should 
have  smaller  thermal  injection  currents  and  be  less  susceptible  to 
second  breakdown  than  the  latter.  The  two-sided  Hi-Lo  or  Lo-Hi-Lo 
diodes  have  even  smaller  temperature  gradients  at  the  low-field  boundary 
and  should  be  the  least  susceptible  to  second  breakdown.  Of  the  one¬ 
sided  diodes,  the  P-diode  should  be  superior  to  the  N-diode,  since  the 
high-temperature  junction  region  of  the  former  is  farther  from  the 
source  of  the  thermally  injected  electrons,  which  ionize  more  easily 
than  do  holes. 

In  conclusion,  we  predict  the  following  order,  from  the  most 
susceptible  to  thermal  damage  to  the  least,  of  doping  profiles: 
(1)  one-sided  N-diodes,  (2)  one-sided  P-diodes,  (3)  linearly  graded 
diodes,  (4)  abrupt  two-sided  diodes,  and  (5)  Lo-Hi-Lo  two-sided 
diodes.  Experimental  tests  will  be  required  to  verify  these  predic- 
ti ons. 


1 1 .  EXTENDING  THE  TIME  SCALE 

The  most  serious  limitation  of  the  computer  program  is  the  cost 
limit  of  simulating  problem  times  of  more  than  a  few  nanoseconds. 
Generally,  the  time  step  is  limited  by  the  dielectric  relaxation  time, 
but  in  all  cases  there  is  also  a  transit-time  limitation  Which  is  also 
important.  At  the  room-temperature  saturation  velocity,  the  transit 
time  is  approximately  10  ps  for  each  micrometer  of  diode  width.  One 
transit-time  effect  is  the  generation  of  avalanche  oscillations.  It 
will  be  shown  in  section  14  that  the  presence  of  oscillations  has  a 
marked  effect  on  the  current-voltage  characteristic  of  a  diode.  At 
times  when  the  current  is  quasi -stationary,  when  dl/dt  <  e  (where  e  may 
be  chosen  by  experience),  it  should  be  possible  to  use  large  values  of 
At  to  solve  the  thermal  diffusivity  equation  (eq  (6))  and  then  use  small 
values  of  At  to  solve  the  electronic  equations  with  the  new  temperature 
distribution. 

11.1  Procedure 


Since  the  present  computer  program  is  limited  to  the  nano¬ 
second  time  regime,  it  is  useful  to  explore  ways  of  obtaining  informa¬ 
tion  at  longer  times.  One  method  is  to  combine  the  calculated  static 
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characteristics  at  various  temperatures  with  the  knowledge  of  the 
temporal  rate  of  temperature  rise.  The  rate  of  temperature  rise  as  a 
function  of  current  density  and  field  for  silicon  was  given  in  figure 
20.  The  procedure  is  illustrated  for  a  4-pm  P-type  diode  doped  to  5  x 
1015  cm""3.  The  300  K,  dV/dt  =  1  x  1010  V/s  characteristic  of  this  diode 
was  shown  in  figure  13.  This  characteristic,  together  with  those  at 
selected  higher  temperatures,  is  shown  in  a  linear  plot  in  figure  52  (p 
60).  Load  lines  corresponding  to  Rg  =  50  fl  and  an  area  of  5  x  10~5  cm? 
are  also  shown  in  the  figure  for  applied  voltages  of  120,  130,  180,  and 
230  V.  The  calculated  temporal  rate  of  temperature  rise  as  a  function 
of  J  is  shown  in  figure  52(d)  for  temperatures  of  300,  400,  and  500  K. 
The  temperature  increase  rate  is  greater  at  higher  temperatures,  since 
higher  fields  are  required  to  obtain  the  same  current  density. 

To  proceed,  consider  the  180-V  load  line  (the  second  from  the 
top  in  fig.  52(a)).  The  current  density  at  300  K  is  3.1  x  10**  A/cm2. 
As  shown  in  figure  52(d),  the  temperature  increases  5  K  per  nanosecond 
at  300  K  for  J  =  3.1  x  104  A/cm-2.  Therefore,  20  ns  are  required  to 
reach  400  K  at  the  applied  voltage.  The  intersection  of  the  load  line 
with  the  400  K  characteristic  gives  new  J  and  V  values.  Thus,  J  versus 
time  curves  are  constructed  point  by  point,  as  shown  in  figure  51(b). 
Similarly,  V  versus  t  curves  are  constructed,  as  illustrated  in  figure 
52(c)  for  an  applied  voltage  of  120  V.  These  curves  show  an  initial 
increase  of  V  and  a  decrease  of  J  with  time,  followed  by  a  decrease  in  V 
and  increase  of  J,  indicating  the  onset  of  second  breakdown.  Other 
applied  voltages  may  be  used  to  extend  the  time-to-breakdown  scale  to 
higher  and  lower  values. 

The  curves  of  figure  52(b)  and  (c)  may  be  used  to  plot  power 
as  a  function  of  time.  A  log- log  plot  of  the  average  power  versus  the 
breakdown  time  is  shown  in  figure  53  (p  61).  Breakdown  time  in  this 
case  is  determined  by  the  time  to  reach  600  K.  The  slope  of  this  curve 
is  1.16.  Experimental  measurements  in  this  time  range  usually  have  a 
slope  of  about  one  on  the  log-log  plot. 

11.2  Simulation  of  1N4148 

Extensive  experimental  studies26  have  been  made  on  second 
breakdown  in  the  switching  diode,  1N4148.  This  diode  has  an  N-Epi 
region  of  11  ym,  doped  to  1.5  £)-cm.  The  diode  used  for  the  measurements 
of  interest26  has  a  3-mil  diffusion  window.  To  simulate  this  diode, 
calculations  were  made  for  a  12-ym  N-type  diode  doped  to  3  x  1015 
cm-3.  The  P  region  was  1  im.  The  diode  area  was  set  at  5  x  10~5  cm2. 

26D.  M.  Tasca,  J.  C.  Peden,  and  J.  Miletta,  Non-Destructive  Screening 
for  Second  Breakdown,  IEEE  Trans.  Nucl.  Sci.,  NS- 19  (December  1972),  57- 
67. 
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Figure  52.  Method  of  obtaining  current  and  volt¬ 
age  versus  time  curves  from  static  characteris¬ 
tics.  (a)  Characteristics  for  4-pm  P-type  diode 
doped  to  5  x  10 15  cm”3.  Parameter  is  temperature 
in  multiples  of  100  K.  Characteristic  for  550  K 
shown  by  unlabeled  dashed  curve.  Dashed  load 
lines  correspond  to  50  ft  and  diode  area  of  5  x 
10"5  cm”2.  (b)  Derived  current  density  versus 
time  curves  with  applied  voltage  as  parameter. 
Dashed  connections  between  (a)  and  (b)  help  show 
construction  process,  (c)  Derived  voltage  versus 
time  curves.  Temperature  at  various  times  can  be 
obtained  from  dashed  connection  to  (a).  Applied 
voltage  is  given  for  each  curve,  (d)  Rate  of  tem¬ 
perature  rise  as  function  of  current  density  that 
was  used  in  construction  of  (b)  and  (c).  Param¬ 
eter  is  temperature  (K). 
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Constant  temperature 
static  characteristics  were 
calculated,  point  by  point,  for 
this  diode  from  300  to  600  K. 
These  are  shown  in  figure  54. 
Calculations  were  also  made  for 
thermal  injection  currents 
corresponding  to  100  K  less 
than  the  bulk  temperature. 
These  curves,  for  temperatures 
from  400  to  650  K,  are  also 
shown  in  figure  54.  The  tem¬ 
perature  difference  of  100  K 
was  chosen  arbitrarily.  Dif¬ 
ferent  temperature  differences 
for  different  curves  might  be 
more  realistic.  It  may  be  seen 
that  the  curves  with  the  lower 
temperature  boundary  condition 
merge  with  the  corresponding 
isothermal  calculation  at  a 
current  density  that  increases 
with  temperature. 

Portions  of  the  two-t^mperature  characteristics  are  replotted 
in  figure  55.  A  linear  scale  is  used  and  the  ordinate  is  total 
current.  Load  lines  for  5000  fl  are  shown  in  figure  55(a).  This  value 
was  chosen  for  agreement  with  the  experimentally  observed26  changes  of 
voltage  and  current  upon  second  breakdown.  The  method  of  calculating 
current  versus  time  and  voltage  versus  time  was  described  in  section 
11.1.  The  results  for  applied  voltages  of  from  500  to  1200  V  are  shown 
in  figure  55(b)  and  (c).  The  area  of  5  x  1 0~6  cm  should  be  used  to 
convert  J  to  I  in  figure  55(d).  As  seen  in  figure  55(b),  the  current  is 
nearly  constant  in  time,  whereas  the  voltage  changes  considerably  with 
time  (and  temperature)  as  seen  in  figure  55(c). 

Prom  the  current  and  voltage  at  each  time  step,  the  average 
power  dissipated  may  be  calculated.  The  time  to  reach  the  temperature 
of  650  K  is  assumed  to  be  the  breakdown  time.  The  power  thus  calculated 
as  a  function  of  the  time  to  breakdown  was  found  to  agree  satisfactorily 
with  experimental  measurements  for  breakdown  times  of  less  than  ~1  ps. 
Power  to  breakdown  obtained  from  the  isothermal  calculations  fell 
considerably  below  the  measured  curve. 


Figure  53.  Power  to  breakdown  ver¬ 
sus  breakdown  time  for  diode  of 
figure  52.  Slope  of  straight  line 
is  1.16. 


26Z5.  M.  Tasca,  J.  C.  Peden,  and  J.  Miletta,  Non-Destructive  Screening 
for  Second  Breakdown,  IEEE  Trans .  Nucl.  Sci.,  NS-19  (December  1972 J,  57- 
67. 
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VOLTAGE  (VI 


Figure  54.  Calculated  static  character¬ 
istic  for  12-ym  N-type  diode  doped  to  3  x 
1015  cm-^,  simulating  1N4148.  Solid 
curves  are  for  isothermal  calculations, 
and  temperature  (K)  is  given  as  param¬ 
eter.  Ebr  dashed  curves,  thermal  injec¬ 
tion  current  corresponds  to  T  -  100,  where 
T  is  bulk  temperature,  again  given  as  pa¬ 
rameter. 


In  the  above  calculations  it  was  assumed  that  the  thermal 
time  constant  was  infinite,  i.e.,  the  temperature  continued  to  increase 
as  if  no  heat  flow  occurred.  The  thermal  time  constant,  rth,  of  the 
1N4148  is  not  available,  but  500  mW  is  given  as  a  safe  operating 
level.  Exploratory  measurements  in  our  laboratory  indicate  that  the 
diodes  fail  at  about  1  to  2  w. 

The  construction  diagram  of  the  1N4148  shows  two  heat  dumet 
heat  sinks;  on  the  small  area  P+  side,  the  sink  is  about  50  ym  from  the 
junction  and  on  the  large  area  (~2  x  10  ^  cm^)  N+  side,  the  sink  is 
about  135  ym  from  the  junction.  Although  the  equations  given  in  section 
4  apply  only  to  one  directional  heat  flow,  their  use  can  give  some  idea 
of  the  thermal  resistance  (0)  and  thermal  time  constants  (Tth) 
expected.  The  dependence  of  and  0  in  silicon  upon  the  distance  to 
the  heat  sink,  L,  is  shown  in  figure  56.  The  equilibrium  junction 
temperature  is  given  by 
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Figure  55.  Derivation  of  current  versus  time  and 
voltage  versus  time  for  1N4148s  (a)  Static  char¬ 
acteristics  of  figure  54,  replotted  on  linear 
scale.  Parameter  is  bulk  temperature  in  multiples 
of  100  K.  Dotted  and  dashed  characteristics  are  for 
550  and  625  K,  respectively.  Area  is  5  x  10-5  cm-3 
and  load  lines  shown  are  5000  ft.  (b)  Resultant  cur¬ 
rent  versus  time  curves.  See  figure  52  caption  for 
details.  (c)  Resultant  voltage  versus  time 
curves.  (d)  Time  derivative  of  temperature  as  func¬ 
tion  of  current  density.  Parameter  is  temperature 
and  unlabeled  curves  are  for  400  and  500  K  in  mono¬ 
tonic  order. 


Tj  =  Tq  +  LP/k.jA  ,  (39) 

as  seen  from  equations  (21)  and  (23).  For  a  power  of  1  W,  L  =»  135  pm, 
and  A  =  5  x  io-5  cm2,  T.  -  TQ  is  calculated  to  be  186  K.  This  is  also 
equal  to  0,  since  P  =  1  W.  The  temperature  dependence  of  kT  is 
neglected.  This  temperature  increase  per  watt  is  in  reasonable 
agreement  with  the  laboratory  measurement  of  1  to  2  W  for  failure.  The 
thermal  time  constant  would  be  dominated  by  the  closer  heat  sink.  From 
equation  (22)  we  find  Tth  -  2.8  x  io”5  s  for  t  «■  50  pm. 
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t/T 


DISTANCE  (mn) 


Figure  56.  Thermal  time  constants  and  thermal  resist¬ 
ance  of  silicon  as  function  of  distance  (left  and 
bottom  scales).  See  equations  (22)  and  (23).  Area 
assumed  is  5  x  10"5  cm2.  Shown  also,  right  and  upper 
scales,  is  the  plot  of  exp(-t/i)  versus  t/x. 


For  a  further  check,  we  can  use 

P0/Tth  =  JE/pc  ,  (40) 


which  is  required  at  t  =  0  from  equations  (25)  and  (26).  From  computer 
calculations,  JE/pc  =  1.5  x  io7  K/s  for  a  power  of  1  W.  Since  P0  is  the 
steady-state  temperature  difference  (eq  (21)),  0  is  better  known  than 
Tth*  To  a9ree  with  the  computer  value  for  dT/dt  at  t  =  0,  one  finds 
from  equation  (40)  that  Tth  =  1.2  x  10-5  s  for  0  =  186  K/W. 


64 


Selecting  1  x  10“5  s  for  xth,  the  computed  power  to  break¬ 
down,  P0,  versus  time  for  breakdown,  t0,  for  the  1N4148  was  recalculated 
using  equation  (26),  The  results  showed  an  appreciable  increase  in  P0 
for  tB  >  Tth/2.  The  experimental  data26  showed  an  appreciable  deviation 
from  the  Pg  «  T-1/2  curve  at  tB  =  1  x  10-6  s.  Therefore,  the  calcu¬ 
lations  were  repeated  with  Tth  =  2  x  io~6  s.  The  results,  shown  in 
figure  57,  show  good  agreement  with  the  experimental  data. 

Considering  the  large  uncertainty  in  r^,  the  factor  of  5 
needed  for  agreement  does  not  seem  disturbing.  Moreover,  it  might  be 
expected  that  would  not  be  constant,  but  vary  with  time,  especially 

since  the  heat  sink  is  dual  and  the  geometry  not  planar. 


Figure  57.  Comparison  of  calcu¬ 
lated  and  experimental  breakdown 
power  for  1N4148  diode.  Range  of 
experimental  data  of  Tasca  et  al 
(ref  26)  is  given  by  error  bars. 
Calculated  values  are  given  by 
crosses  for  Tth  =  “,  by  pluses 
for  xth  =  1  x  io“5  s,  and  by 
triangles  for  tth  =  2  x  10“6  s. 
Straight  line  is  drawn  with  slope 
of  -1  and  indicates  constant 
energy. 


While  the  calculated  PB 
agrees  with  the  measured  values 
for  tg  <  1  x  10” 6  without  adjust¬ 
ment,  the  calculated  breakdown 
currents,  Ig,  are  lower  than  those 
measured  by  about  40  percent. 
This  is  shown  in  figure  58.  There 
are  a  number  of  plausible  expla¬ 
nations  for  this  discrepancy. 
First,  the  doping  level,  3  x  10*  5 
cm-3,  chosen  for  the  calculations 
applies  to  ideal  1.5  (2-cm  silicon; 
a  higher  doping  level  correspond¬ 
ing  to  lower  mobility  material 
would  increase  the  current  at  a 
given  voltage.  Second,  the 
spreading  resistance  of  this  diode 
is  less  than  the  calculated  space- 
charge  resistance  in  the  SCL 
current  regime.  Since  the  area  of 
the  N-type  Epi  region  is  40  times 
that  of  the  P-dif fusion  window, 
current  spreading  should  be  ex¬ 
pected  to  occur.  Recall  that  the 
maximum  field  at  the  junction  is 
nearly  independent  of  current 
density,  but  the  field  near  the 
edge  of  the  depletion  region 
increases  sharply  with  current 


26D.  M.  Tasca,  J.  C.  Peden,  and  J.  Miletta,  Non-Destructive  Screening 
for  Second  Breakdown,  IEEE  Trans,  Nu cl,  Sci,,  NS- 19  (December  1972),  57- 
67. 
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density,  both  in  the  SCL  regime. 
Again,  current  spreading  would 
increase  the  current  at  a  given 
voltage.  Third,  the  calculations 
of  the  static  characteristics  were 
made  with  high  diffusion  coeffi¬ 
cients  to  deliberately  damp  out 
IMPATT  oscillations  which  would 
otherwise  appear.  The  presence  of 
avalanche  oscillations  results  in 
higher  average  currents  at  a  given 
voltage.  Avalanche  oscillations 
will  be  discussed  further  in  sec¬ 
tion  14.  Finally,  some  of  the 
material  parameters  chosen  may  be 
erroneous.  In  particular,  the  use 
of  the  ionization  coefficients 
measured  by  Grant27  reduces  the 
breakdown  voltages  by  about  5  to 
10  percent,  depending  on  the  diode 
doping.  Again,  higher  currents 
would  be  expected  at  a  given  volt¬ 
age. 


TIME  bd) 

Figure  58.  Comparison  of  calcu¬ 
lated  and  measured  breakdown  cur¬ 
rents  for  1N4148  diode.  Open 
circles  with  error  bars  give  mea¬ 
surements  of  Tasca  et  al  (ref 
26).  Crosses  indicate  values  cal¬ 
culated  with  =  “>  and  pluses 

indicate  Tth  =  2  x  10“^  s 

calculations. 


12.  IONIZING  RADIATION 

In  section  7  it  was  mentioned  that  the  boundary  condition  on  the 
current  density,  JQ,  could  be  used  to  simulate  the  effect  of  ionizing 
radiation  upon  the  diode.  Only  the  injected  photoionization  current  has 
been  simulated.  Bulk  photoionization  in  the  high-field  region  is  not 
simulated,  but  since  JQ  is  multiplied  by  a  multiplication  factor  of  from 
102  to  10i(J  after  breakdown,  the  bulk  photoionization  may  be  neglected. 

12.1  Two  Examples 


The  calculated  characteristics  of  section  8  may  be  used  to 
study  the  effect  of  ionizing  radiation  upon  second  breakdown.  Portions 
of  the  characteristics  of  figure  25  have  been  redrawn  on  a  linear  scale 
in  figure  59.  An  area  of  1  x  i  O-4  cm2  has  been  assumed  and  50-ft  load 
lines  are  also  shown.  Arrows  indicate  currents  for  the  voltage 
maxima.  For  illustrative  purposes,  let  us  consider  a  half-cycle  sinu¬ 
soidal  voltage  with  an  amplitude  of  160  V  and  a  half-period  of  180  ns. 
With  no  ionizing  radiation  it  is  assumed  that  JQ  =  1  x  i 0~7  A/cm2,  or  IQ 
=  1  x  io"n  A.  From  figure  59,  no  appreciable  current  flows  until  the 

2 7W.  nT  Grant,  Electron  and  Hole  Ionization  Rates  in  Epitaxial 
Silicon  at  High  Electric  Fields,  Solid-State  Electron, ,  (1973),  1189- 

1203, 
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Figure  59.  Calculated  character¬ 
istics  of  4-pm  N-type  diode  doped 
to  5  x  10 crtf^  with  ionizing 
radiation  as  parameter.  Same 
characteristics  shown  in  figure 
25(b),  but  plotted  on  linear 
scale.  Area  of  1  x  IO-1*  cm2  is 
assumed.  50-il  load  lines  are 
shown  for  applied  voltages  of  120 
and  160  V. 


VOLTAGE  <V| 

applied  voltaqe,  V  ,  exceeds  110  V.  Then  the  diode  current  and  voltage 
are  determined  by  the  intersection  of  the  static  characteristic  and  the 
moving  load  line.  This  may  be  done  graphically  at  desired  time 
intervals.  The  resulting  current-versus-ti me  and  diode-voltage-versus- 
time  curves  are  plotted  in  figure  60(a)  together  with  V  as  applied. 
Similar  curves  simulating  ionizinq  radiation,  with  Io  =  1  x  10  c  A,  are 
plotted  in  figure  60(b).  The  presence  of  the  ionizinq  radiation  results 
in  an  increase  in  the  peak  current  of  about  20  percent  and  of  peak  power 
by  18  percent.  At  first  thought,  one  might  easily  conclude  that  the 
presence  of  ionizinq  radiation  increases  the  susceptibility  to  damage. 
However,  note  that  for  IQ  =  1  x  io-11  A,  a  negative  differential  resis¬ 
tance  (NOR)  reqime  is  reached  at  I  H  0.5  A,  whereas  the  NDR  regime  is 
only  reached  at  I  =  1.7  A  for  the  I  =  1  x  io"2  A  curve.  Thus,  for  ID  = 
1  x  io-11  A,  the  current  exceeds  the  threshold  for  the  NDR  regime  for  a 
period  of  60  ns,  whereas  for  I  =  1  x  io  2  A,  the  maximum  current  is  well 
below  the  NDR  threshold.  Since  current  contraction  into  a  channel  is 
likely  in  the  NDR  reqime,  the  presence  of  ionizing  radiation  may 
decrease  the  susceptibility  of  some  diodes  to  second  breakdown. 


Figure  60.  Current  and  voltage  as  function  of  time  for  half-cycle 
sinusoidal  applied  voltage  (Va)  pulse  shown.  Characteristics  of 
figure  59  were  used  to  deduce  these  curves.  (a)  I  =  1  x  10”11 
A.  Negative  resistance  is  encountered  in  "unsafe”  area.  (b)  I  = 
10'2  A.  Negative  resistance  area  is  off  scale. 


In  the  above  calculation  we  assumed  a  constant  room  temper¬ 
ature  characteristic  and  a  moving  load  line  to  calculate  I  and  V  versus 
t.  In  the  last  section  we  used  a  stationary  load  line  and  a  moving 
characteristic  (as  a  function  of  temperature)  to  calculate  dynamic 
curves.  The  two  techniques  may  be  readily  combined  to  allow  both  the 
load  line  and  the  static  characteristic  to  move  in  time  in  order  to 
simulate  a  varying  applied  voltage  and  a  heating  of  the  diode,  respec¬ 
tively. 


A  second  illustrative 
calculation  assumes  a  constant 
applied  voltage  of  120  V,  but  IQ 
is  varied  with  time.  The  pulse 
of  ionizing  radiation  is  assumed 
to  cause  I  to  rise  exponential¬ 
ly  from  1  10'"  to  1  x  10~2  A 
in  18  ns,  to  reach  a  peak  of  1  x 
10"1  A  at  24  ns,  decay  to  1  x 
10“2  A  at  60  ns,  and  thereafter 
decay  by  a  factor  of  10  each  10 
ns.  The  variation  of  I  with 
time  is  shown  in  figure  61.  The 
calculation  is  again  made  graph¬ 
ically  with  the  use  of  figure 
59,  and  the  results  are  also 
plotted  in  figure  61.  The  peak 
power  during  the  irradiation 
pulse  is  greater  by  a  factor  of 
just  over  2.5  than  the  non- 
irradiated  power.  In  this  case, 
the  susceptibility  to  second 
breakdown  is  increased  by  the 
presence  of  the  ionizing  radi¬ 
ation. 

In  the  second  illustrative  calculation  in  the  paragraph  just 
above,  we  assumed  a  constant  temperature,  a  stationary  load  line,  and  a 
characteristic  moving  as  a  function  of  IQ.  In  order  to  further  include 
heating  during  the  ionizing  pulse,  a  full  set  of  curves  such  as  figure 
57  would  be  required  at  each  temperature  chosen.  It  would  seem  that  a 
requirement  for  information  on  a  particular  diode  and  a  particular 
ionizing  pulse  would  be  needed  in  order  to  justify  such  a  complex  calcu¬ 
lation. 

12.2  Discussion 


TIME  (ns) 

Figure  61.  Calculated  current, 
voltage,  and  power  for  same  diode 
(fig.  59)  with  V,  =  120  V  and 

u 

exposed  to  ionizing  radiation 
pulse  described  in  text.  Unmulti¬ 
plied  injection  current  is  labeled 


There  are  very  limited  experimental  reports  upon  the  combined 
effects  of  an  ionizing  radiation  environment  and  high  electrical  stress 
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on  semiconductors.  Vault28  found  that  the  effect  of  ionizing  radiation 
was  to  reduce  the  power  level  for  second  breakdown,  with  the  reduction 
greater  for  reverse  bias  than  for  forward  bias.  However,  Habing29  had 
earlier  found  that  ionizing  radiation  lowered  the  pulse  current  failure 
threshold  of  a  forward-biased  transistor  base-emitter  junction  by  a 
factor  of  five,  but  made  no  change  in  the  reverse-biased  threshold. 
Hartman  and  Evans88  found  that  the  same  level  of  ionizing  radiation  had 
minimal  effect  upon  the  burnout  thresholds  of  several  transistor  base- 
emitter  junctions.  There  have  also  been  oral  reports  of  ionizing  radi¬ 
ation  increasing  the  level  of  stress  required  for  junction  failure. 

In  order  to  compare  the  above  calculations  with  experimental 
measurements,  it  is  necessary  to  determine  the  level  of  photocurrents 
that  result  from  given  dose  rates.  The  chief  unknown  in  this  calcu¬ 
lation  is  the  effective  volume  in  which  ionization  occurs.  As  argued 
above,  ionization  that  occurs  outside  the  depletion  region  is  more 
effective  through  avalanche  multiplication  than  is  the  ionization  that 
occurs  in  the  high-field  region.  However,  Habing29  calculates  a  photo¬ 
ionization  current  of  100  mA  for  his  device  with  a  dose  rate  of  5  x  1011 
rad(Si)/s.  This  is  just  the  maximum  initiating  photocurrent  that  was 
assumed  in  figure  61. 

In  experimental  measurements  of  power-to-fail,  one  usually 
adjusts  the  power  measured  to  a  selected  time-to-fail  by  means  of  a 
power  law  relation  between  the  two.  This  is  necessary  since  neither 
power-  nor  time-to-fail  can  be  preselected.  However,  this  procedure 
could  erroneously  eliminate  any  effect  of  ionizing  radiation  upon 
second-breakdown  susceptibility.  For  example,  in  the  illustrative 
calculation  above  of  the  effect  of  a  pulse  of  ionizing  radiation  (fig. 
60),  an  increased  power  was  calculated.  However,  the  time-to-fail  would 
also  be  decreased.  The  fact  that  the  data  fit  the  same  power-versus- 
time  curve  masks  the  fact  that  the  ionizing  radiation  did  modify  the 
power  dissipated.  This  consideration  makes  the  experimental  determi¬ 
nation  of  the  effect  of  ionizing  radiation  upon  stress  measurements 
extremely  difficult. 

A  second  indication  of  possible  decreased  susceptibility  to 
second  breakdown  due  to  ionizing  radiation  is  indicated  from  calcula¬ 
tions  made  of  switching  from  forward  to  reverse  bias.  These  calcula¬ 
tions,  which  will  be  discussed  in  section  16,  have  shown  that  the  stored 

2fy.  L.  Vault,  Response  of  Integrated  Circuits  to  a  Simulated  IEMP 
Environment,  Harry  Diamond  Laboratories,  HDL-TR-1649  (October  1973), 

2%>.  H •  Habing,  Response  of  Bipolar  Transistors  to  Combined  EMP  and 
Ionization  Environment,  IEEE  Trans,  Nucl ,  Sci,,  NS-1 7  (December  1970), 
360-363. 

30b.  f.  Hartman  and  D.  C.  Evans,  Electrical  Pulse  Burnout  of 
Transistors  in  Intense  Ionizing  Radiation ,  IEEE  Trans.  Nucl.  Sci.,  NS-22 
(December  1975),  2528-2532. 
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charges  remaining  from  the  forward  conduction  prevent  the  maximum  field 
at  the  junction  from  attaining  as  high  a  value  as  was  found  in  calcula¬ 
tions  for  a  fully  depleted  diode.  Similarly,  the  mobile  space  charge 
due  to  ionizing  radiation  could  lower  breakdown  voltages.  In  the  case 
of  extremely  high  levels  of  ionizing  radiation,  one  would  expect  little 
difference  between  forward  and  reverse  bias  susceptibility. 

As  far  as  is  known,  this  is  the  first  theoretical  attempt  to 
quantitatively  calculate  the  combined  effect  of  ionizing  radiation  and 
high  electrical  stress  upon  reverse-biased  silicon  diodes.  It  has  been 
shown  that,  depending  upon  the  devices  and  experimental  conditions 
chosen,  susceptibility  to  high  electrical  stress  may  be  increased  or 
decreased  by  the  presence  of  ionizing  radiation.  There  is  some  exper¬ 
imental  support  of  this  conclusion. 


13.  DISCUSSION  OF  NEGATIVE  RESISTANCE 

Negative  static  and  dynamic  characteristics  have  been  calculated  in 
virtually  all  diodes  simulated.  The  threshold  for  the  NDR  regime  varies 
from  around  1  04  A/cm2  for  near  intrinsic  material,  to  a  few  times  10s 
A/cm2  for  material  doped  to  1017  cm-3.  By  definition,  the  static  char¬ 
acteristic  separates  a  lower  voltage  region,  where  currents  decrease  in 
time,  from  a  higher  voltage  region  where  currents  increase  in  time. 
Therefore,  dynamic  characteristics  cross  the  static  characteristic  at 
maximum  or  minimum  points  in  curren-,  except  for  the  unusual  case  when 
the  mobile  space-charge  distribution  is  drastically  different  from  the 
steady  state.  In  contrast,  dynamic  characteristics  attain  their  maximum 
and  minimum  voltages  when  they  cross  the  load  line.  These  rules  will  be 
illustrated  in  the  later  section  on  avalanche  oscillations. 

The  stability  of  a  diode  for  currents  on  the  positive  resistance 
region  is  much  greater  than  for  currents  in  the  NDR  regime.  Assume 
that,  for  reason  of  material  inhomogeneity  or  temperature  nonuniformity, 
for  example,  the  current  density  of  a  local  area  is  greater  than  that  of 
the  surrounding  region.  If  the  current  is  in  the  positive  resistance 
regime,  like  point  A  in  figure  62,  the  local  area  is  undervoltaged,  and 
its  current  will  decrease.  This  stabilizing  effect  helps  keep  the 
current  density  uniform  despite  material  defects  and  temperature  gra¬ 
dients.  However,  if  the  current  is  in  the  NDR  regime,  like  point  B  in 
figure  62,  the  local  area  is  overvoltaged  and  the  current  will  increase 
further.  Thus,  the  small  area  will  rob  current  from  surrounding  areas 
and  form  a  channel  or  filament.  One-dimensional  calculations  are  no 
longer  valid  when  f ilamentation  occurs.  However,  once  f ilamentation  is 
initiated,  damage  will  occur  if  the  current  is  not  limited  sufficiently 
by  the  external  circuit. 
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The  size  of  the  external  re¬ 
sistance  is  of  utmost  importance 
when  operating  in  the  negative 
resistance  region.  Two  load  lines 
are  shown  in  figure  62.  If  the 
external  resistance  is  less  than  the 
absolute  magnitude  of  the  negative 
resistance,  as  for  Rg  =  5  ft  in  this 
illustration,  the  current  will  jump 
from  the  point  of  tangency  (point  C) 
to  intersect  the  static 
characteristic  at  a  much  higher 
current  (point  D) .  If  the  external 
resistance  is  greater  than  the 
magnitude  of  the  negative  resistance 
and  the  load  line  intersects  the 
static  characteristic  in  the 
negative  resistance  region  (point  E 
for  Rs  =  50),  oscillations  will 
result.  These  will  be  discussed  in 
the  next  section. 


Figure  62.  Illustration  of  di¬ 
ode  stability  in  negative  dif¬ 
ferential  resistance  region. 
Characteristic  was  calculated 


14.  AVALANCHE  OSCILLATIONS 


14.1  Background 
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Several  times  in  the  pre¬ 
ceding  sections,  avalanche  oscilla¬ 
tions  have  been  mentioned.  While  it 
is  not  appropriate  to  discuss  these 
oscillations  fully  in  this  report, 
there  are  important  consequences  of 
these  oscillations  for  second  break¬ 
down.  Read31  originally  proposed  a 
high-frequency  negative  resistance  diode  in  1958.  The  IMPact-ionization 
Avalanche  Transit  Time,  or  IMP ATT,  mode  was  first  observed  by  Johnston 
et  al32  in  1964.  The  anomalous,  or  high-efficiency,  mode  was  first 
observed  by  Prager  et  al33  in  1967,  and  was  shown  to  entail  a  trapped 
plasma  and  given  the  acronym  TRAPATT  by  Johnston  et  al.34  The  relaxing 


3 V.  T.  Read,  A  Proposed  High-Frequency  Negative  Resistance  Diode, 
Bell  Syst .  Tech.  J.,  37_  (March  1958),  401-446. 

32r.  L.  Johnston,  B.  C.  DeLoach,  Jr.,  and  B.  G.  Cohen,  A  Silicon 
Diode  Microwave  Oscillator,  Bell  Syst.  Tech.  J.,  47  (1964),  369-372. 

3%.  J.  Prager,  K.  K.  N.  Chang,  and  S.  Weisbrod,  High-Power,  High- 
Efficiency  Silicon  Avalanche  Diodes  at  Ultra-High  Frequencies,  Proc. 
IEEE  (Lett.),  55_  (April  1967),  586-587 . 

3<4r.  l.  Johnston,  D.  L.  Scharfetter,  and  D.  J.  Bartelink,  High- 
Efficiency  Oscillations  in  Germanium  Avalanche  Diodes  Below  the  Transit- 
Time  Frequency ,  Proc.  IEEE  (Lett.),  56  (September  1968),  1611-1613 . 
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avalanche  mode  (RAM)  was  calculated  by  Ward  and  Udelson^  in  1968*  A 
variation  of  the  latter  mode,  named  the  MULTIPATT  mode ,  was  reported  by 
Ward36  in  1973.  These  modes  will  be  described  in  the  following  para¬ 
graphs  . 

1 4. 2  Calculations 

In  section  11,  in  discussing  the  calculations  simulating  the 
1N4148,  it  was  mentioned  that  avalanche  oscillations  were  present  unless 
prevented  by  increasing  the  diffusion  coefficients  and/or  the  external 
RgC  time  constant.  Examples  of  these  calculated  oscillations  are  shown 
in  figure  63;  the  series  resistance,  Rs,  was  chosen  to  be  50  Si,  and 
various  shunt  capacitances,  C,  of  from  one  to  several  picofarads  were 
selected.  Conventional  IMPATT  oscillations  are  calculated  with  applied 

voltaqes  below  210  V  at  300  K.  The 
lower  current  limit  has  not  been  inves¬ 
tigated.  These  oscillations  are  shown 
in  figure  63  as  the  small  ellipses 
joined  by  a  dashed  line.  In  this 
applied  voltage  range,  IMPATT  oscilla¬ 
tions  appear  centered  about  the  inter¬ 
section  of  the  static  characteristic 
and  the  load  line,  but  leave  the  static 
characteristic  as  they  grow  in  ampli¬ 
tude  and  average  current.  Equilibrium 
oscillations  are  attained;  these  are 
those  which  are  shown  in  figure  63.  At 
higher  applied  voltages,  the  IMPATT 
oscillations  continue  to  grow  in  ampli- 
Figure  63.  Various  avalanche  tude;  equilibrium  is  not  attained.  The 
oscillations  calculated  for  current  increases,  and  the  voltage  de- 
simulated  1N4148  diode.  creases  in  a  sequence  of  loops  along 
Static  characteristics  are  and  above  the  load  line.  When  the 
shown  for  300  and  400  K.  voltage  across  the  diode  can  no  longer 
Small  ellipses  indicate  IMPATT  sustain  the  current,  the  current  decays 
oscillations.  A  MULTIPATT  rapidly,  looping  below  the  load  line, 
oscillation  is  shown  for  300  The  voltage  recovery  is  determined  by 
K,  V  =  220  V,  and  Rg  =  50  n  the  external  RgC  time  constant.  When 
(load  line  is  shown).  A  RAM  the  current  and  voltage  return  to  the 
oscillation  is  shown  for  400  load  line  near  the  static  character- 
K,  V  =  260  V,  and  Rg  =  50  SI •  istic,  the  oscillation  is  complete. 

35A.  L.  Ward  and  B.  J.  Udelson,  Computer  Calculations  of  Avalanche- 
Induced  Relaxation  Oscillations  in  Silicon  Diodes,  IEEE  Trans.  Electron 
Devices,  ED-15  (November  1968),  847-851. 

36A.  L.  Ward,  Further  Calculations  of  the  Relaxation  Avalanche  Mode, 
Proc.  4th  Cornell  Electrical  Engineering  Conference,  Ithaca,  NY  (1973), 
391-400. 


This  mode  of  oscillation  has  been  named  the  MULTIPATT  mode,  since  mul¬ 
tiple  transit  time  oscillations  are  superimposed  on  each  RAM  oscilla¬ 
tion.  Near  its  threshold,  the  MULTIPATT  oscillation  is  extremely  noisy; 
an  example  is  shown  in  figure  64.  She  variation  of  the  oscillation 
period  is  mainly  due  to  the  residual  amplitude  of  the  IMPATT  oscillation 
as  it  returns  to  the  load  line.  In  calculations  where  the  temperature 
is  increased  as  power  is  dissipated,  the  oscillations  may  abruptly 
cease.  The  similarity  to  microplasma  behavior  is  obvious. 


Figure  64.  MULTIPATT  current  oscillations  calculated 
for  300  K,  V  =  232  V,  R  =  50  S  ,  and  C  =  1  x  10" 12  F. 


At  voltages  well  above  the  threshold,  the  MULTIPATT  period 
becomes  much  more  constant.  The  number  of  transit  time  oscillations'  per 
MULTIPATT  period  averages  about  20  for  this  diode  and  circuit  combina¬ 
tion.  Other  combinations  have  produced  from  2  to  over  30  smaller 
current  pulses  superimposed  on  the  large  current  pulse.  For  the  ratio 
of  2  and  3,  the  corresponding  number  of  bunches  of  electrons  and  holes 
are  found  to  be  in  transit  across  the  diode  at  a  given  time.  For  higher 
ratios,  the  number  of  bunches  present  at  any  time  is  usually  about  3  or 
4.  An  example  of  the  electron  distributions  at  selected  times  for  the 
simulated  1N4148  is  shown  in  figure  65. 


As  the  applied  voltage  is  increased,  the  maximum  current  and 
voltage  excursion  increase.  In  some  cases  no  equilibrium  is  reached 
before  the  computations  become  unstable  and  are  halted.  Ttie  highest 
current  peak  in  equilibrium  calculations  was  1.5  x  105  A/cm2  and  the 
minimum  voltage  was  94  V.  The  latter  is  well  below  the  V_.  of  125  V. 
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Figure  65.  Distribution  of 
electrons  across  diode  during 
current  growth  stage.  Inset 
shows  portion  of  character¬ 
istic,  whose  numbered  times 
correspond  to  distributions 
shown.  Calculations  are  for 
300  K,  V  =  300  V,  Rs  =  50  (1, 
and  C  =  3  x  10“12  F. 


At  still  higher  applied  voltages,  the  MULTIPATT  mode  oscillations  decay 
into  a  new  mode  of  small-ampli tude  transit  time  oscillations,  as  shown 
in  figure  63.  This  mode  is  termed  the  traveling-wave  IMPATT  mode, 
because  multiple  bunches  of  carriers  are  traversing  the  diode  at  any 
time.  At  still  higher  currents,  the  carrier  bunches  are  confined  to  two 
high-field  regions  near  each  (P+N  and  NN+)  junction.  Ibis  mode,  called 
the  trapped-plasma  IMPATT  mode,  is  also  shown  in  figure  63.  It  requires 
a  much  higher  voltage  than  does  the  traveling -wave  IMPATT  mode,  but  less 
than  the  nonoscillating  state.  These  modes  have  been  described  more 
fully. 25 


IMPATT  oscillations  were  alsb  calculated  for  this  diode  at 
400  K  at  both  low  and  high  currents.  No  oscillations  occurred  at  inter¬ 
mediate  currents.  At  the  higher  temperature,  the  amplitudes  of  the 
oscillations  are  smaller.  In  the  NDR  regime,  suppressing  the  oscil¬ 
lations  resulted  in  higher  voltages  required  to  obtain  the  same  current, 
as  can  be  seen  in  figure  63. 

The  variation  of  frequency  of  the  IMPATT  mode  with  current 
density  is  shown  in  figure  66  for  several  temperatures.  The  decrease  in 
frequency  with  temperature  is  noted  experimentally  in  pulsed  IMPATTs. 
It  results  from  the  decrease  in  saturation  velocity  with  temperature. 
The  straight  lines  in  figure  66  are  drawn  with  a  slope  of  1/2.  There  is 
an  approximate  fit  over  two  orders  of  magnitude  of  J.  The  half-power 


25A.  L.  Ward,  Modes  of  Avalanche  Oscillations  in  Silicon  Diodes,  IEEE 
Trans,  Electron  Devices,  ED-25  (June  1978),  683-687. 


Figure  66.  Variation  of  IMPATT 
frequency  with  current  density 
at  indicated  temperatures. 
Straight  lines  are  drawn  with 
slope  of  1/2. 


law  has  been  observed  experimentally  and  attributed16  to  the  constancy 
of  the  inductance-current  (LI)  product,  since  the  frequency  is  given  by 

f  cc  (LC)-1/2  , 

and  C  is  a  constant.  The  constancy  of  the  LI  product  was  verified  by 
calculations  for  the  RAM  oscillations.36 


In  the  absence  of  IMPATT  oscillations,  steady-state 
RAM  oscillations  can  only  occur  in  the  negative  resistance  regime.  A 
RAM  oscillation  is  shown  for  400  K  in  figure  63.  That  these  are  relax¬ 
ation  oscillations  is  more  evident  from  plots  of  I  and  V  as  a  function 
of  time.  These  are  shown  in  figure  67  for  a  near-intrinsic  2.5-pm 
diode.  It  is  seen  that  the  period  is  essentially  determined  by  the  time 
required  for  the  diode  to  reach  the  breakdown  voltage,  Vg,  from  the 

minimum  voltage,  V  . 

m 

The  time  for  a  capacitance,  C,  to  charge  up  to  Vg  from  Vm 
through  Rg  with  an  applied  voltage,  Vg,  is 

T  =  RgC  *n[(va  -  Vj/(va  -  Vg)]  .  (41) 


1%.  M.  Sze,  Physics  of  Semiconductor  Devices,  Wiley-Interscience , 
New  York  (1969),  39-41,  57-59 ,  121-126. 

36A.  L.  Ward,  Further  Calculations  of  the  Relaxation  Avalanche  Mode, 
Proc .  4th  Cornell  Electrical  Engineering  Conference,  Ithaca,  NY  (1973), 
391-400. 
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Figure  67.  Typical  relaxation  oscillations.  Calculation  is  for 
2. 5-vim  P-type  diode  with  area  of  4  x  10_l4  cm2  and  doped  to  2  x 
1014  cm-^ .  (a)  Temporal  variation  of  voltage  across  and  current 

through  diode.  Shown  on  enlarged  scale  is  current  through  series 
resistor.  (b)  Dynamic  current-voltage  trace  for  data  shown  in 
(a).  Data  points  are  for  time  intervals  of  5  *  10_12  s. 

V. 

However,  VB  and  V  are  not  constant,  but  depend  on  dV/dt,  for  example, 
which  is  a  function  of  RgC.  An  extensive  series  of  calculations  was 
made  to  test  equation  (41).  The  diode  used  was  a  nearly  intrinsic  5-um 
diode.  For  large  amplitude  oscillations,  equation  (41)  was  satisfied 
when  R  and  C  were  varied  only  when  V„  and  Vm  were  adjusted  to  their 

S  dm! 

calculated  values. 

By  reducing  the  external  RgC  time  constant,  oscillations  were 
found  to  decay  into  nearly  sinusoidal  oscillations.  To  analyze  these 
decaying  (or  growing)  sinusoidal  oscillations,  the  diode  may  be 
considered  to  be  an  inductance,  L,  representing  the  avalanche  delay,  and 
a  differential  resistance,  R^,  representing  the  dV/dl  of  the  diode. 
Then  it  can  be  shown37  that  current  excursions  from  the  steady-state 
value  are  given  by 

I(t)  =  I(o)[sin  (wt)]  exp [At]  ,  (42) 


where 


-i(3t 

2  \L  R  c/ 


37A.  L,  Ward  and  L.  G.  Schneekloth,  Calculations  of  Relaxation 
Oscillations  in  Gas  Tube  Circuits,  Harry  Diamond  Laboratories,  HDL-TR- 
1166  (August  1963). 
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and 
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1 

LC 


(44) 


From  equations  (42)  and  (43),  it  is  seen  that  oscillations  will  decay 
(A  <  0)  if  R.  is  positive,  or,  for  R.  negative,  if  RgC  is  less  than 
| L/Rd | ,  the  diode  time  constant.  Oscillations  will  grow  (A  >  0)  if  R^ 
is  negative  and  RgC  >  jL/R^I.  Equations  (43)  and  (44)  may  be  solved  for 
Rd  and  L  to  obtain 


and 


-R  C 
s 

1  +  2R  CA 
s 


(45) 


-R3(l  +  2RSCA) 

1  +  2RgCA  +  R|c2(w2  +  A2) 


(46) 


These  two  equations  have  been  used  to  analyze  the  near- 
sinusoidal  oscillations  as  computed  with  varying  Rg,  C,  and  V  .  The 
calculated  period  for  the  near-intrinsic  5-ym  diode  is  shown  in  figure 
68(a)  as  a  function  of  the  total  capacitance,  the  external  C  plus  the 
diode  self-capacitance.  The  range  of  C  was  from  1  x  10~13  to  1  x  )0~11 
F,  the  range  of  Vg  was  from  130  to  170  V,  and  Rg  was  kept  constant  at 
500  0.  The  straight  lines  through  the  data  at  each  value  of  V  have  a 
slope  of  1/2.  The  dashed  line  separating  growing  oscillations  (• 
points)  and  decaying  oscillations  (o  points)  is  drawn  at  a  slope  of 
0.625. 


A  plot  of  the  average  current  against  Vg  showed  that  the 
current  extrapolated  to  zero  at  127.5  V,  the  breakdown  voltage.  Figure 
68(b)  shows  a  plot  of  the  calculated  frequency,  f,  against  V  -  127.5 
for  C  =  2  x  10  12  F  and  Rg  =  500  0.  The  straight  line  has  a3  slope  of 

V2.  Since  1  is  proportional  to  Vg  -  127.5  this  again  shows  that 

f  *  1  /  •  Also  shown  in  figure  68(b)  is  a  plot  of  f  against  R  for  v 

=  145  V  and  C  =  1  x  10~12  F.  The  straight  line  is  drawn  with  aSslope  of 

-1/2.  The  information  in  these  two  figures  may  be  included  in  the 

single  expression 


f  -  9.5  x  103(va  -  127.5)l/2/(R8C)1/2  .  (47) 
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Figure  68.  Plots  of  oscillation  period  and  fre¬ 
quency.  (a)  Oscillation  period  as  function  of  total 
capacitance  for  Rs  *  500  ft.  Log  scales  are  used. 
Open  circles  denote  decaying  and  solid  circles  denote 
growing  oscillations.  Parameter  of  each  cur.ve  is 
Vft.  (b)  Oscillation  frequency  against  -  127.5  for 
Rs  =  500  ft  and  C  =  2  x  10~12  F,  and  against  Rfi  for  Va 
=145  V  and  C  =  1  x  10"1Z  F. 


Comparison  of  equations  (41)  and  (47)  does  not  give  much 
promise  of  compatibility.  However,  with  the  end  point  known,  equation 
(41)  can  be  manipulated  to  give  the  form  of  equation  (47).  Consider  the 
log  term  of  equation  (41): 


V 

a 


V 

a 


V 

a 


+  v. 


V 

m 


(48) 


Consider  the  term  Vg  -  Vm,  the  maximum  voltage  excursion  of  the  oscil¬ 
lation.  This  should  be  equal  to  AQ/C,  where  AQ  is  the  charge  flow  for 
the  voltage  difference.  This  charge  flows  through  Rg,  so  we  may  expect 
A Q  to  be  inversely  proportional  to  Rg.  Then, 

VB  -  Vm  =  const/RgC  .  (49) 
Substituting  equations  (48)  and  (49)  into  equation  (41),  one  finds 


T  =  R  C  Jin 
s 


const j 

+  RsC(va  -  VB)_ 


(50) 


Now  we  introduce  an  approximation  that 

£n(1  +  x)  ®  0 . q/x  .  (51) 


Plots  of  fn(1  +  x)  and  0.8/x  are  shown  in  figure  69.  It  is  seen  that 
the  approximation  is  a  good  one  for  0  <  x  <  ~8.  Identifying  1  +  x  with 
(va  ~  vmV(va  “  vb ) *  we  find  x  >  0,  since  Vg  >  vm,  and  x  of  more  than  8 
would  be  unlikely.  For  the  illustration  in  figure  67,  Va  =  100,  VB  a 
90,  and  Vm  =  30.  This  gives  x  =  7,  one  of  the  highest  values  yet  calcu¬ 
lated.  Applying  the  approximation  of  equation  (51)  to  equation  (50),  we 
find 


or 


T 


const2 


f 


const3 


fva  -  vb)l/2 

(RgC)l/2 


(52) 


which  is  the  expression  sought,  equation  (47).  Therefore,  we  may 
identify  relaxation  oscillations  in  semiconductor  diodes  by  shunting  the 
diode  with  various  capacitances  and  measuring  the  frequency  variation. 
IMPATT  frequencies  are  unaffected  by  the  external  circuit. 
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Figure  69.  Plots  of  £n(1  +  x) 

and  of  0.8v^x  as  a  function  of  x. 


The  growing  or  decaying 
oscillations  for  Rg  =  500  ft  and  C 
=  2  x  10“ 12  F  were  analyzed  by 

determining  the  growth  rate,  X, 
and  a)  as  a  function  of  Vg.  Then 
R^  and  L  were  calculated  from 
equations  (45)  and  (46)  and  are 
shown  in  figure  70.  Also  shown  in 
the  same  plot  is  the  product  of  L 
and  the  average  current,  I.  This 
product  is  seen  to  be  constant,  as 
stated  previously.  Figure  71 

shows  plots  of  the  diode  time 
constant,  L/|Rdl,  and  X  as  a 
function  of  Vg.  This  figure  shows 
that  X  is  negative  (decaying 
oscillations)  for  L|Rdl  >  RgC  =  1 
x  10“9  s. 
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Figure  70.  Calculated  values  of 
negative  resistance,  inductance 
(L),  and  LI  as  function  of  Va  for 
diode  of  figure  68.  Rg  =  500  ohms 
and  C  =  2  x  10“ 12  F. 
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Figure  71.  Calculated  time  con¬ 
stant  ( L/ | R | )  and  growth  constant 
(X)  as  function  of  V#.  Circuit 
parameters  are  given  in  figure  70 
caption. 
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14.3  Experimental  Measurements 

During  the  course  of  EMP  vulnerability  studies,38 
oscillations  were  observed  in  the  reverse-biased  T1DF606A  diode 
(equivalent  to  1N645  and  so  designated  hereafter)  at  stresses  just  below 
second  breakdown.  Oscillations  were  noted  in  the  vhf/uhf  range;  figure 
72  shows  current  amplitudes  of  greater  than  1  A  and  voltage  excursions 
of  over  50  V.  The  onset  of  the  large  amplitude  oscillations  is  delayed 
about  400  ns  (in  this  example)  from  the  time  the  voltage  reaches  its 
nominal  maximum.  This  indicates  that  these  are  device  oscillations,  not 
merely  circuit  oscillations.  Various  capacitors  were  shunted  across  the 
diode  to  test  for  relaxation-type  oscillations.  The  dependence  of  the 
frequency  upon  the  shunt  capacitance  is  shown  in  figure  73.  The 
straight  line  drawn  through  the  data  has  a  slope  of  -1/2,  in  agreement 
with  equation  (52)  for  the  relaxing  mode. 


VCLTBGE 


Figure  72.  Typical  pulsed  avalanche  oscillation  observed  in  1N645. 
Time  scale,  200  ns/div.  Top  trace;  voltage,  200  V/div.  Bottom 
trace:  current,  1.0  A/div. 


38R.  (/,  garver,  R.  /(.  Fenenias,  Vasco  C.  Martins,  .7.  L.  Washington, 
Gary  L.  Roffman,  R.  J.  Royzer ,  Samuel  A.  Clark,  Jr.,  and  Morris  Campi , 
Systems  Performance  Analysis  of  the  Electromagnetic  Pulse  Effects  on 
Defense  Communications  System  Transmission  Systems,  Final  Report  (U ) , 
Harry  Diamond  Laboratories,  HDL-TR-19S1  (February  1982).  (SECRET- 
RESTRICTED  DATA/NOFORN ) 


Figure  73.  Frequency  of  oscillations  in 
1N645  as  function  of  C.  Straight  lines 
are  drawn  with  slope  of  -1/2.  Multiple 
calculated  frequencies  arise  from 
dependence  upon  amplitude  of  oscillations. 

The  doping  profile  of  the  1N645  is  not  known.  The  devices 
were  measured  in  our  laboratory  to  have  a  VR1  of  from  about  500  to  600 
V.  Since  the  oscillations  occur  at  about  BOO  V  we  know  that  the  device 
is  not  punched  through  at  breakdown.  From  fiqure  40  it  is  seen  that  a 
dopinq  level  of  about  4  x  I0il*  cm-3  will  give  Vgl  -  500  V.  Figure  38 
shows  that  a  diode  width  of  at  least  50  pm  would  be  needed.  Calcula¬ 
tions  were  made  for  an  N-type  50-gm  diode  doped  to  4  x  i01Lt  cm’^.  The 
calculations  for  this  diode  showed  stronq  IMPATT  and  MULTIPATT  oscilla¬ 
tions.  The  current-density  threshold  for  MULTIPATT  oscillations  was 
just  under  lx  1  o3  A/cm2.  Comparing  this  threshold  with  the  exper¬ 
imental  threshold  of  just  under  1  A,  as  seen  in  figure  72,  we  find  an 
area  of  about  1  x  1 0-3  cm2.  Capacitance  versus  reverse  voltage  measure¬ 
ments  of  the  1N645  were  consistent  with  the  deduced  doping  and  area.  As 
noted  in  previous  sections,  oscillations  usually  may  be  eliminated  in 
calculations  if  desired.  Static  characteristics  thus  calculated  at 
three  temperi Lures  are  shown  in  figure  74.  These  characteristics  are 
composites  of  both  dynamic  and  steady-state  calculations,  chosen  for 
best  accuracy  in  several  current  ranges.  The  voltage  increases  about  10 
V  per  decade  increase  in  current  in  the  SCF  current  regime,  which  agrees 
well  with  laboratory  measurements. 

Both  IMPATT  and  relaxation  oscillations  were  calculated  for 
the  1N645.  A  sequence  of  oscillations  is  shown  in  figure  75,  together 
with  the  static  characteristic  and  the  appropriate  50-fi  load  lines. 
IMPATT  oscillations  (small  ellipses)  are  noted  at  both  lower  and  higher 
currents  than  the  MULTIPATT  relaxation  osci  nations .  The  calculated 
current  versus  time  and  voltage  versus  time  for  an  applied  voltage,  Va, 
of  738.4  V  are  shown  in  figure  76.  The  2.2-GHz  IMPATT  oscillation  is 
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Figure  74.  Calculated  static  char¬ 
acteristics  for  simulation  of 
1N645.  Parameter  is  temperature  (K). 


VOLTAGE  (V) 

Figure  75.  Calculated  oscilla¬ 
tions  for  1N645  simulation.  Tem¬ 
perature  is  300  K,  A  =  1  x  10"3 
cm2,  Rs  =  50  Q,  C  =  1  x  10“12 

F.  MULTIPATT  types  of  relaxation 
oscillations  occur  at  inter¬ 
mediate  applied  voltages;  IMPATT 
oscillations  at  both  higher  and 
lower  voltages. 


Figure  76.  Voltage  and  current  oscillations  calculated  for 
1N645.  Change  of  frequency  in  time  is  due  to  change  in  am¬ 
plitude  of  oscillations. 


seen  to  grow  slowly  to  the  current  peak.  Reducing  Vg  a  few  tenths  of  a 
volt  extends  this  time  appreciably.  It  is  believed  that  growing  IMPATT 
oscillations  occur  during  the  200-ns  delay  before  large-amplitude  oscil¬ 
lations  are  seen  experimentally  in  figure  72.  The  frequency  of  the 
relaxation  oscillations  in  figure  76  increases  from  145  to  182  MHz  as 
time  increases.  The  last  current  pulse  continued  to  grow  until  the 
calculations  were  halted  due  to  instability. 

One  major  factor  in  the  calculated  frequency  change  of  the 
oscillations  is  the  amplitude  of  the  oscillations.  This  is  complicated 
by  mode  changes  which  are  not  fully  understood.  Unexplained  mode 

changes  are  prevalent  among  IMPATT  diode  measurements.  The  effect  of  a 
temperature  increase  in  the  diode  is  also  quite  noticeable.  Figure  77 
shows  that  a  rise  of  only  6  K  of  the  average  temperature  across  the 
diode  reduces  the  frequency  from  ~317  to  ~310  MHz.  During  oscillations, 
heating  is  uniform  neither  in  time  nor  in  space.  The  maximum  and 

average  temperatures  across  a  diode  are  shown  in  figure  78(a)  as  a 

function  of  time.  The  temperature  distribution  across  the  diode  at 
selected  times  is  shown  in  figure  78(b).  The  heating  is  quite 

nonuniform. 


The  dependence  of  the  RAM  frequency  upon  the  current  level 
and  changes  of  modes  complicates  the  study  of  the  dependence  o£ 
frequency  upon  shunt  capacitance.  The  calculated  frequency  of  the  1N645 
as  a  function  of  C  is  given  in  figure  73.  There  is  about  a  factor  of 
five  range  of  calculated  frequency  for  each  capacitance.  Generally,  the 
higher  the  average  current,  the  higher  the  resulting  frequency,  but  mode 
changes  produce  exceptions.  The  highest  calculated  frequencies  obey  the 
inverse  square  root  of  C  quite  well,  the  lowest  frequencies  less  well. 
It  should  be  mentioned  that  the  experimentally  measured  oscillations 
were  also  subject  to  mode  changes.  In  particular,  there  were  observably 
different  oscillations  of  a  lower,  but  irregular,  frequency.  Those 
plotted  in  figure  73  were  of  a  uniform  amplitude. 


Figure  77.  change  of  frequency  with  temperature.  (a)  Current  as 
function  of  time.  Frequency  drops  from  317  MHz  at  300  K  to  310  MHz 
at  averag>  temperature  of  306  K.  (b)  Voltage  as  function  of  time. 
Dashs  -3  «,_.rve  was  calculated  with  constant  temperature  of  300  K; 
so? id  curve  was  calculated  with  increasing  temperature. 
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Figure  78.  Diode  temperature  in  time  and  space  for  run  of  figure  77. 
(a)  Average  and  maximum  temperature  as  function  of  time.  Arrows  indi¬ 
cate  times  of  current  peaks.  (b)  Distribution  of  temperature  across 
diode  at  selected  times.  Arrow  at  left  top  indicates  position  of  metal¬ 
lurgical  junction. 


Restricting  the  calculations  to  a  given  mode  and  external 
circuit,  the  RAM  frequency  was  found  to  increase  with  the  square  root  of 
the  average  current  density.  This  is  shown  in  figure  79.  However,  the 
frequency  of  the  IMPATT  mode  of  this  50-gm  diode  increased  with  only 
about  the  0.4  power  of  J,  as  seen  in  figure  79.  Also  shown  in  figure  79 
are  the  IMPATT  oscillations  of  a  12-um  diode  and  the  RAM  oscillations  of 
a  5-pm  diode,  illustrated  earlier.  Further  shown  in  figure  79  are  two 
points  for  higher  temperature  oscillations  in  the  1N645,  the  only  cases 
where  oscillations  have  been  observed  at  these  or  higher  temperatures. 
The  modes  and  circuit  are  not  comparable  to  the  300  K  data. 

It  is  well  known  that  the  frequency  of  IMPATT  oscillations  is 
strongly  temperature  dependent,  but  the  frequency  range  of  5  GHz  and 
above  precludes  using  this  variation  as  a  practical  dynamic  temperature 
probe  of  the  junction.  Our  measurements  and  calculations  show  that 
relaxation  oscillations  of  less  than  100  MHz  may  be  used  as  a 
temperature  probe.  Figure  80  shows  an  example  of  the  decrease  in 
frequency  observed  in  a  long  (4-ys)  pulse.  A  change  of  oscillation  mode 
is  also  evident. 

Low-amplitude  oscillations  were  also  observed  in  the  1N4148 
during  pulse  testing  in  the  laboratory. 
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Figure  79.  Variation  of  frequency 
with  current  density  for  several 
diodes.  IMPATT  frequencies  at  300 
and  400  K  are  shown  for  12-um  diode 
(1N4148)  and  at  300,  400,  and  500  K 
for  50-ym  diode  (1N645).  Relaxation 
oscillations  are  shown  for  5-um  and 
50-pm  diodes.  Solid  lines  are  drawn 
with  slope  of  1/2. 
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Figure  80.  Observed  change  of  frequency  of  re¬ 
laxation  oscillation  with  temperature  in 
1N645.  Diode  was  shunted  by  940  pF.  Upper 
traces  200  V/cm;  lower  traces  2  A/cm.  Time 
scales  0.5  ys/cnt. 


15.  FORWARD  BIAS  CHARACTERISTICS 
15.1  Zero  Current 

Before  discussing  calculations  of  forward  bias  character¬ 
istics,  it  is  useful  to  discuss  the  condition  of  zero  current  flow, 
i.e.,  no  applied  bias. 

For  J  =  0,  equation  (3)  becomes 

ne^E  =  eD(dn/dx)  ,  (53) 


where  the  subscripts  have  been  dropped.  Solving  for  E  we  have 

E  =  (d/jj)  (dn/dx)/n  =  (kT/e )  (dn/dx)/n  ,  (54) 

where  the  last  form  follows  from  the  Einstein  relation.  For  room 
temperature,  kT/e  =  0.026  V.  Neglecting  p,  equation  (5)  becomes 

dE/dx  =  e(N  -  n)/<EQ  .  (55) 

The  sign  convention  has  been  changed  back  to  the  conventional. 


The  last  two  equations  can  be  solved  for  the  two  unknowns,  E 
and  n,  by  a  direct  difference  substitution.  That  is. 

An  =  (e/kT)nE(Ax)  ,  (56) 

AE  =  e(n  -  N)Ax/xeQ  .  (57) 

The  r roblem  is  determining  the  boundary  conditions.  First,  consider  a 
step  junction  at  x  =  0.  That  is,  N  =  NQ  for  x  <  0,  and  N  =  0  for  x  > 
0.  We  know  E  =  0  at  x  =  ±®  and  n  =  NQ  at  x  =  -»  and  n  =  0  at  x  =  +», 
but  these  conditions  cannot  be  used  to  solve  equations  (56)  and  (57). 
However,  we  can  start  computations  from  x  =  0  and  proceed  in  each 
direction  in  turn.  The  problem  is  that  E(0)  and  n(0)  are  both 
unknown.  One  could  guess  that  n{0)  =  NQ/2,  choose  Ax,  a  trial  value  of 
E(  0) ,  and  proceed  to  use  each  equation  alternately.  Even  a  pocket 
calculator  would  quickly  convince  one  that  n(0)  <  NQ/2.  In  fact,  it  was 
soon  found  that  n{0)  =  NQ/e.  This  should  have  been  anticipated,  since 

n(x)  =  No  exp[-ev(x)/kT]  (58) 


87 


and  V( 0)  =  kT/e  for  a  step  junction.  The  value  of  E(0)  and  V(0)  for 
various  PN  junctions  may  be  found  in  Kennedy,^  where  it  is  shown  that 
V(0)  “  e/kT  for  P  >>  N  or  N  >>  P.  For  Nq  =  1  x  10*9  cm-3  and  Ax  =  1  x 
10-7,  the  value  of  E(0)  was  found  to  be  1.67  x  io5  V/cm,  again  using  a 
pocket  calculator.  This  field  maximum  agrees  with  Kennedy's  value  for  N 
>>  P.  The  field  distribution  for  the  Nq  =  1  x  10* 9  cm-3  step  junction 
is  shown  in  figure  81  .  The  total  high-field  space-charge  region  is 
about  3  x  io-6  cm  wide.  This  reverse-bias  field  everywhere  exactly 
counterbalances  the  diffusion  current  that  otherwise  would  flow  (eq 
(53)).  The  junction  potential  is  known  as  the  built-in,  diffusion,  or 
contact  voltage.  It  cannot  be  measured  directly  because  of  other 
contact  potentials  to  the  measuring  instrument. 


Figure  81.  Field  distributions  for 
step  junctions  with  zero  current 
flow.  Junction  is  at  x  =  0.  Doping 
on  x  <  0  side  is  labeled  for  each 
curve.  Doping  is  zero  for  x  >  0. 


Field  distributions  for 
other  step  junctions  are  shown  in 
figure  81.  It  may  be  noted  that 
for  N  =  N^IOO,  E  is  reduced  by  a 
factor  of  10  and  x  is  increased  by 
a  factor  of  10.  This  is  easily 
verified  by  substituting  in  equa¬ 
tions  (56)  and  (57);  the  equations 
are  unchanged.  The  variation  of 
E(  0 )  with  for  a  step  junction 

is  shown  in  figure  82. 

A  comment  should  be 

made  on  guessing  the  boundary 
field,  E(0).  Consider  the  x  >  0 
direction  where  Nq  =  0.  If  E  is 
too  high,  then  An  is  too  high  by 
equation  (56).  This  means  that  n 
is  too  low  in  equation  (57); 

therefore,  AE  is  too  low.  The 

eventual  result  is  that  n  becomes 
negative.  Conversely,  for  E  too 
small.  An  is  too  low  or  n  too 
high.  It  follows  that  AE  is  too 
large  and  E  becomes  negative. 
Similar  criteria  apply  to  the  x  < 
0  direction;  N  either  exceeds  N 

o 

or  decreases,  or  E  becomes  nega¬ 

tive  or  changes  direction.  For 
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programming  a  computer,  equation 
(58)  could  be  used  as  a  test  for 
the  initial  boundary  condition. 

A  second  type  of  junc¬ 
tion  studied  is  the  linear  grad¬ 
ient  of  width  d.  This  is  easily 
programmed,  since  NQ  changes  a 
constant  amount  each  Ax  step. 
However,  the  boundary  conditions 
are  now  more  tricky.  For  very 
steep  gradients  the  midpoint  of 
the  gradient  may  be  used,  and 
values  of  n(0)  and  E(0)  may  be 
estimated  from  the  step  junc¬ 
tion.  FOr  less  steep  gradients, 
the  edge  of  the  gradient  where  N  = 
Nq  should  be  chosen  for  x  =  0. 

Then  a  guess  of  n(0)  fixes  E(0)  in 
the  x  <  0  direction.  Calculations 
converge  when  E(d)  and  n(d)  agree 
with  a  possible  solution  to  the 
step  junction  N  =  0  region. 


The  field  and  carrier 
distribution  for  a  graded  junction 
with  Nq  =  1  x  10 153  and  d  =  2  x 

10-f>  cm  is  shown  in  figure  83. 
The  maximum  field  for  the  graded 
junction  is  only  6.1  x  101*  V/cm  as 
compared  to  1.67  x  105  V/cm  for 
the  step  junction. 

The  two  equations  were 
programmed  for  a  desk  calculator 
and  the  accuracy  of  n(0)  =  NQ/e 

and  the  fit  to  equation  (58)  veri¬ 
fied  to  several  more  significant 
figures.  Even  for  these  more 
accurate  calculations,  E  becomes 
sufficiently  small  while  n  is 
still  so  much  greater  than  n^  that 
the  neglect  of  p  is  justified. 


Figure  82.  Junction  field  for 
step  junctions  with  zero  current 
flow  as  function  of  doping 
level.  Equation  of  straight 
line  is  E  =  5.4  x  1 0-5  . 


Figure  83.  Field  distribution  for 
linearly  graded  junction  with  zero 
current  flow.  Doping  is  shown  as 
dashed  curve;  equilibrium  carrier 
density  as  dotted  curve.  Field 
maximum,  solid  curve,  is  near  zero 
doping  intercept. 
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This  discussion  of  the  field  distribution  for  J  =  0  has  been 
presented  to  introduce  the  fact  that,  for  reasonable  currents  the  field 
in  regions  of  doping  gradients  will  be  negative  with  zero  or  small 
forward  bias.  The  current  density  needed  to  assure  E  =  0  when  the 
gradient  is  dn/dx  is  given,  from  equation  (3),  by 

J  =  eD(dn/dx)  .  (59) 

Since  D  »  40  for  electrons,  for  AN  =  1  x  1016  in  1  pm  (a  very  small 

gradient  of  1  x  102°  cm-*4),  J  »  640  A/cm2.  Thus,  for  smaller  current 
densities,  the  field  will  remain  negative  in  the  junction  region,  but 
will  be  positive  where  there  is  a  smaller  carrier  gradient.  The 
presence  of  fields  of  both  polarities  causes  problems  in  the  present 
computer  program. 

When  Nq  is  only  a  few  orders  of  magnitude  greater  than  n^ ,  it 
will  be  necessary  to  include  p  in  equation  (55).  Then  p  can  be 
calculated  from  the  relation 


p  =  n?/n  , 

and  p  and  n  will  each  approach  n^  as  E  0. 

For  very  small  applied  voltages,  equation  (56)  can  be 
modified  to  include  J,  as  follows. 

An  =  (e/kT)nE( Ax)  ±  (J/eD)(Ax)  ,  (60) 

where  the  positive  si<m  pertains  to  forward  bias  and  the  negative  sign 
to  reverse  bias.  The  ^ign  is  verified  by  setting  E  =  0  (valid  only  for 
a  forward  bias)  and  regaining  equation  (59). 

No  exploratory  calculations  have  been  made  with  p  or  J  not 
equal  to  zero. 

15.2  Static  Characteristics 

The  initial  conditions  chosen  to  calculate  forward-bias 
static  characteristics  are  those  that  are  exactly  SCF;  i.e.,  the  mobile 
carriers  are  set  equal  to  the  doping  density.  The  applied  voltage  is  an 
important  parameter  in  the  initial  transient.  For  low  voltages, 
carriers  cannot  flow  in  unlimited  densities  from  the  highly  doped 
region,  but  are  limited  by  space  charge.  After  the  carriers  cross  the 
diode,  the  charges  neutralize  and  the  densities  continue  to  grow.  The 
formation  of  a  negative  field  at  one  or  both  high  doping  gradients 
eventually  causes  instability  in  the  present  computer  program  (which  was 
not  designed  for  forward-bias  calculations),  and  the  calculations  are 
halted.  By  reducing  the  doping  levels  and/or  gradients  in  the  highly 
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doped  regions,  the  problem  is  reduced.  Further,  at  high  enough 
currents,  the  negative  field  disappears,  as  discussed  in  the  first  part 
of  this  section.  Therefore,  second-breakdown  studies  are  not  impeded. 

When  higher  applied  voltages  are  used  for  forward-bias  calcu¬ 
lations,  initial  space-charge  limiting  is  reduced,  and  the  current  grows 
to  an  equilibrium  value  quite  rapidly.  A  large  series  resistance 
permits  a  high  initial  voltage  and  a  sharp  voltage  drop  as  the  current 
increases.  The  current  overshoots  and  approaches  equilibrium  from 
above.  It  has  been  verified  in  a  few  cases  that  the  same  equilibrium  is 
reached  whether  the  equilibrium  current  is  approached  from  above  or 
below.  The  negative  field  problem  is  reduced  by  the  approach  from 
above . 


The  carrier  distribution  during  the  turn-on  transient  for  a 
12-im  PIN  diode  is  shown  in  figure  84.  The  maximum  doping  was  held  to  1 
x  10  36  cm-3  in  the  P+  and  N+  regions.  The  calculations  were  made  for  Va 
=  20  V,  Rg  =  100  ft,  and  an  area  of  1  x  lO-*4  cm2.  The  corresponding 
field  distributions  are  shown  in  figure  85.  The  voltage  and  current 


transients  are  shown  in  figure  86. 
attained. 
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Figure  84.  Distribution  of  car¬ 
riers  during  turn-on  tran¬ 
sient.  Doping  profile  is  shown 
as  solid  curves,  electrons  as 
dashed  curves,  and  holes  as 
dotted  curves.  Parameter  is 
time  in  ps. 


Reasonable  equilibrium  has  been 
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Figure  85,  Distribution  of 
field  during  turn-on  tran¬ 
sient.  Parameter  is  time  in 
ps. 
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Figure  86.  Voltage  and  current  turn-on  transients, 
(a)  Diode  voltage  and  (b)  diode  current  density  as 
function  of  time.  (c)  Current  density  as  function  of 
voltage;  100-S2  load  line  is  shown.  (d)  Portion  of 
forward  characteristic;  V  =  20  and  Rg  =  100,  200,  and 
500  «. 


Whereas  doping  levels  are  the  dominant  factor  in  reverse 
breakdown,  they  are  a  minor  factor  in  forward-bias  characteristics.  The 
two  dominant  factors  are  diode  width  and  the  recombination  lifetime. 
Calculated  forward  characteristics  for  different  diode  widths  are  shown 
in  figure  87.  The  accuracy  of  these  calculations  at  the  lowest  currents 
is  questionable,  but  the  strong  dependence  on  diode  width  is  observed 
experimentally. 


A  much  more  extensive  study  has  been  made  of  the  dependence 
of  forward  characteristics  on  recombination  lifetimes,  t.  Static  char¬ 
acteristics  calculated  for  a  12-ym  PIN  diode  for  various  recombination 
lifetimes  are  shown  in  figure  88.  For  most  of  the  data,  the  conduction 
current  nearly  equals  the  total  current,  but  at  the  lowest  currents 
shown,  diffusion  currents  are  an  appreciable  part  of  the  total  current. 


VOLTAGE  (V) 


Figure  87.  Forward  characteris¬ 
tics  for  several  diode  widths. 
Diode  width  in  pm  is  parameter. 


1.0  10  100  1000 
VOLTAGE  (V) 


Figure  88.  Forward  character¬ 
istics  for  12-pm  PIN  diode  for 
recombination  lifetimes  of  1,  10, 
and  100  ps.  Straight  lines  have 
slopes  of  1  or  2.  l 


For  still  lower  currents,  not  calculable  with  the  present 
program,  the  diffusion  component  is  nearly  the  total  current,  and  the 
current  increases  exponentially  with  voltage.  A  recent  paper40  presents 
calculations  of  forward  characteristics  up  to  J  =  10  A/cm2  and  V  S  1  V, 
showing  the  deviation  from  exponential  growth  at  the  highest  currents. 


It  is  noted  in  figure  88  th^t  the  current  increases  as  the 
square  of  the  voltage  for  over  an  order  of  magnitude  (in  the  low  current 
regime)  for  t  =  100  ps ,  and  over  smaller  ranges  for  shorter  lifetimes. 
At  higher  currents,  there  is  a  range  of  ohmic  conductivity,  i.e., 
I  «  V.  For  the  higher  values  of  r,  there  follows  a  region  of  sublinear 
current  increase  with  voltage.  At  a  diode  voltage  of  about  100  V,  the 
field  attains  the  avalanche  threshold,  and  the  current  again  increases 
with  the  square  of  the  voltage.  There  is  no  sign  of  a  negative  resis¬ 
tance  regime,  even  though  the  current  is  more  than  a  factor  of  10  larger 
than  the  threshold  of  the  NDR  regime  with  reverse  biasing. 


The  field  and  carrier  distributions  across  the  diode  vary 
markedly  as  the  current  increases  and  are  quite  different  for  different 
lifetimes.  The  carrier  densities  drop  monotonically  across  the 
intrinsic  region  at  all  currents  for  t  =  1  ps,  as  shown  in  figure 


40P.  C.  H.  Chan  and  C.  T.  Sah,  Experimental  and  Theoretical  I-V 
Characteristics  of  Zinc-Doped  Silicon  p-n  Junctions  Using  Exact  DC 
Circuit  Model,  IEEE  Trans.  Electron  Devices,  ED-26  (June  1979),  937-941. 


93 


89(a).  Recombination  prevents  nearly  all  carriers  from  completely 
traversing  the  intrinsic  region.  The  resulting  space  charge  causes 
large  fields  to  develop,  as  shown  in  figure  89(b).  For  t  =  10  us,  many 
more  carriers  completely  traverse  the  intrinsic  region  and  accumulate  in 
the  opposite  high  doping  region  where  they  recombine;  this  is  shown  in 
figure  90(a).  The  net  space  charge  is  now  smaller,  resulting  in  smaller 
fields  for  the  same  current,  as  shown  in  figure  90(b).  For  x  =  100  us, 
the  total  number  of  minority  carriers  in  the  highly  doped  region  is 
greater  than  those  in  the  much  wider  intrinsic  region.  Figure  91(a) 


shows  that  the  carrier  densities  of 
equal  across  the  intrinsic  region, 
at  the  same  current  density,  as  seen 


<b)  DISTANCE  (urn) 


Figure  89.  Field  and  carrier  distri¬ 
butions  for  t  =  1  us.  (a)  Solid 

curves--electrons;  dashed  curves — 
holes;  parameter — current  (mA).  Area 
is  1  *  10-4  cm2.  (b)  Field 

distributions:  parameter-diode  volt¬ 

age  . 
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the  electrons  and  holes  are  near 
Again,  the  field  is  greatly  reduced 
in  figure  91(b). 


Figure  90.  Field  and  car¬ 
rier  distributions  for  x  = 
10  ys.  (a)  Solid  curves — 
electrons;  dashed  curves — holes; 
parameter — current  (mA).  Area 
is  1  x  10~4  cm2.  (b)  Field 
distributions:  parameter-diode 

voltage . 


In  published  theoretical 
studies  of  PIN  diode  switching,  it 
has  been  generally  assumed  that  p  = 
n  in  the  intrinsic  region  and  that 
there  are  no  accumulated  minority 
charges  in  the  junction  region. 
These  calculations  indicate  that 
minority  charges  accumulate  in  the 
junction  region  in  any  condition 
where  n  and  p  could  be  considered 
equal . 

The  maximum  field  in  the 
forward-biased  diode  increases  with 
current  as  seen  in  figures  89  to 
91.  Plots  of  the  maximum  field  as  a 
function  of  J  for  each  t  are  shown 
in  figure  92.  A  power  law  region  is 
evident  for  each  tj  the  slope  varies 
from  about  0.63  for  t  =  100  ps  to 

I. 0  for  t  =  1  us.  The  falloff  of 

EMAX  at  J  is  due  to  avalanch¬ 

ing.  Note  that  there  is  a  higher 
threshold  for  avalanching  for  a 
smaller  T.  An  example  of  calculated 
fields  with  and  without  avalanching 
is  shown  in  figure  93(a).  The 
corresponding  carrier  densities  are 
shown  in  figure  93(b),  When  E  = 
eMAX'  p  =  n  in  the  intrinsic 
region.  The  value  of  ne<^  r  p  =  n  is 
shown  in  figure  94  as  a  function  of 

J.  Again,  power  laws  are  noted  at 
small  currents.  The  slopes  vary 
from  0.77  to  0.92.  At  higher  cur¬ 
rents,  the  curves  merge  into  a 
single  curve  with  a  slope  of  1.07. 


0 1 - ‘ -  j - 1 _ I 

0  3  8  8  12 

DISTANCE  (pm) 


It  is  also  of  interest  to  know  how  the  maximum  field  varies 
with  n  •  This  is  shown  in  figure  95.  The  j/lopes  of  the  straight  lines 
as  drawn  vary  from  0.84  to  1.21.  There  is  ^ome  supporting  evidence  that 
for  t  =  100  us,  the  sharp  dropoff  of  EMflX  f't  the  lowest  ngC!  densities  is 
real.  Calculations  can  be  made  to  higher  t  values  for  smaller  diode 

as  a 
to 

the  best  silicon.  Also  shown  are 
there  is  a  strong  indication  that 


widths.  Figure  96  shows  a  plot  of  Em, 
diode,  with  x  varying  from  2.5  x  10' 
value  of  x  is  that  given  by  Sze1*’  fi 
a  few  calculations  for  a  3-um  diodi 


function  of  neg  for  a  4-pm 
2.5  x  10~3  s.  The  highest 


M.  S ze.  Physics  of  Semiconductor  Devices,  W iley-Interscience , 
New  York  (1969),  39-41,  57-59,  1$1-126. 
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eMAX  droP8  sharply  at  the  lowest  currents  calculable,  i.e.,  when 
diffusion  currents  begin  to  dominate.  This  is  to  be  expected,  since 
diffusion  is  independent  of  field.  In  all  cases,  in  both  figures  95  and 
96,  the  dropoff  of  Ejj^  occurs  at  approximately  E  *  5  *  103  V/cm.  It 
may  or  may  not  be  significant  that  this  is  about  the  field  where  the 
electron  velocity  begins  to  saturate;  see  figure  1. 


All  the  data  thus  far  have  been  for  diodes  at  room 
temperature.  Temperature  effects  are  due  to  the  variation  of  mobility 
and  of  the  intrinsic  density,  n^,  with  temperature.  The  latter  effect 
is  evidenced  in  the  recombination  rate;  see  equation  (10).  A  plot  of  J 
versus  V  for  a  12~Mm  PIN  diode  is 
shown  in  figure  97  for  temperatures 
up  to  400  K,  t  *  1  ys.  Comparing 
figure  88  with  figure  97,  we  see  that 
the  400  K,  t  =  1  s,  curve  follows 
very  closely  with  the  300  K,  t  = 

100  ps,  curve.  The  decreased  V 
required  to  obtain  a  given  J  at 
higher  temperatures  is  due  to  the 
reduced  recombination.  Calculations 
were  also  made  for  450  K,  but  the 
points  showed  considerable  scatter. 

Generally,  the  450  K  curve  followed 
the  400  K  curve,  showing  that  the 
changes  due  to  mobility  and  lifetime 
temperature  dependencies  tended  to 
counterbalance  one  another.  For  long 


Figure  92.  Maximum  field  as  func¬ 
tion  of  current  density  for  three 
recombination  lifetimes •  Slopes 
of  straight  lines  are  given. 


Figure  93.  Comparison  of  field 
and  carrier  distributions  cal¬ 
culated  with  and  without  avalanche 
multiplication.  (a)  Electron  and 
hole  distributions  with  avalanche 
multiplication  are  shown  as  Bolid 
curves  and  as  dashed  curves  with¬ 
out  avalanche  multiplication.  Pa¬ 
rameter  gives  current  in  A.  (b) 
Solid  curve  gives  field  with  ava¬ 
lanche  multiplication  and  dashed 
curve  without.  Parameter  is  diode 
voltage. 
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lifetime  material,  the  temperature  effect  on  recombination  is  small,  so 
an  increase  in  V  would  be  required  to  obtain  a  given  J ,  due  to  the 
decrease  in  mobility  with  temperature.  The  variation  of  resistivity 


with  temperature,  shown  in  figure  3, 
biased  diode. 


CUMEMT  DENSITY  (A/M*) 

Figure  94.  Carrier  density  at  po¬ 
sition  where  electron  and  hole 
densities  are  equal ,  as  function 
of  current  density  for  three  re¬ 
combination  lifetimes.  Parameters 
are  slopes  of  straight  lines. 


DENSITY  (cm-*) 


Figure  96.  Maximum  field  as  func¬ 
tion  of  equal  carrier  density  for 
3-  and  4-iim  diodes.  Curve  labeled 
3  ym  is  for  recombination  lifetime 
of  2.S  ms.  All  other  curves  are 
for  4-ym  diode  and  recombination 
lifetime  is  given  in  ys. 


is  thus  correlated  in  the  forward- 


owsity  (cm-)) 


Figure  95.  Maximum  field  as  func¬ 
tion  of  equal  carrier  density  at 
that  position.  Diode  width  is  12  yn 
and  three  recombination  wlifetimes 
are  used  in  calculations.  Slope  of 
each  straight  line  is  indicated. 
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Figure  97.  Forward  characteris¬ 
tics  for  12-ym  diode  at  different 
temperatures.  Recombination  life¬ 
time  is  1  ys. 
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16.  REVERSE  SWITCHING  TRANSIENTS 


16.1  Low- Voltage  Switching 

The  steady-state  distributions  of  electrons  and  holes  across 
the  forward-biased  diode  are  used  as  initial  conditions  for  simulation 
of  reverse  switching  transients.  For  applied  reverse  voltages  less  than 
VB1'  the  current  will  decay  monotonically  with  time  to  the  reverse 
saturation  current.  For  the  PIN  diode,  the  turn-off  occurs  in  two 
stages  and  is  known  as  a  step  recovery.  During  the  charge  storage  time, 
the  reverse  current  is  essentially  constant,  controlled  by  the  series 
resistance  and  the  applied  voltage.  When  the  carriers  begin  to  deplete 
the  diode,  the  current  decays,  and  the  field  concurrently  builds  up 
across  the  depletion  region. 

An  example  of  the  calculated  decay  of  current  in  a  12-pm  PIN 
diode  is  shown  in  figure  98.  The  calculations  are  for  a  temperature  of 
400  K.  Two  applied  voltages  are  used,  as  well  as  two  initial  values  of 
neq  (i.e.,  n  =  p  in  the  intrinsic  region).  The  recombination  lifetime 
(t)  is  2.5  ms  and  the  series  resistance  is  100  Q.  The  distribution  of 
charges  across  the  diode  is  shown  at  selected  times  in  figure  99. 
Comparing  figures  98  and  99  we  see  that  the  initial  constant  current 
period  corresponds  to  the  recovery  of  the  charges  stored  in  the  highly 
doped  regions.  Published  theories  of  the  step-recovery  diode  do  not 
include  these  stored  charges,  but  identify  the  end  of  the  constant 
current  period  with  the  drop  to  zero  of  the  charges  at  the  infinitely 
narrow  NI  junction.  Since  they  assume  that  all  the  stored  charge  is  in 
the  intrinsic  region,  they  deduce  a  charge  level  in  the  I  region  several 
times  too  high.  The  field  deduced  thereupon  is  similarly  too  low. 
Invariably,  much  more  charge  is  in  the  constant  current  period  than  in 
the  current  decay  period.  Many  calculations,  not  presented  here,  verify 
that  there  is  a  one-to-one  relation  between  the  initial  charge  in  the 
high  doping  region  and  the  charge  flowing  in  the  constant  current 
stage.  Calculations  made  with  differing  reverse  voltages  (see  fig.  98) 
and  different  resistances  (not  illustrated)  verify  that  the  amount  of 
charge  recovered  is  essentially  constant,  independent  of  the  current 
level.  The  amount  of  charge  recovered  during  the  current  decay  is  found 
to  agree  with  the  total  initial  charge  in  the  intrinsic  region.  This  is 
also  illustrated  in  figure  98.  If  T  is  decreased  to  2.5  x  i o“4  from  2.5 
x  io“ 3  s,  some  measurable  charge  is  lost  by  recombination,  and  less  is 
collected  than  was  stored.  However,  identifying  the  measured  reverse 
recovery  time  of  a  diode  with  t  is  a  too  prevalent  erroneous  practice. 

The  field  distributions  corresponding  to  the  carrier 
densities  shown  in  figure  99  are  presented  in  figure  100.  The  second 
sharp  decay  in  current  observed  in  figure  98  is  identified  in  figures  99 
and  100  as  the  transition  from  SCL  to  8CF  current  flow.  This  transition 
has  not  been  calculated  previously.  For  a  3-j a  Pill  diode,  the  current 
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reverse  switching  of  12-jim  PIN 
diode.  Solid  curves  are  for  ini¬ 
tial  equal  density  in  the  intrin¬ 
sic  region  of  2.4  x  10 15  cm2, 
whereas  dashed  curve  is  for  ini¬ 
tial  density  of  5.4  x  lO1^  cm2. 
Parameter  for  each  curve  is 
applied  voltage. 


tions  during  reverse  re¬ 
covery  of  diode  of  figure  98 
with  15  V  applied.  Solid 
curves  are  for  electrons  and 
dashed  curves  are  for 
holes.  Parameter  is  time  in 
ns.  Distributions  of  major¬ 
ity  carriers  in  high  doping 
regions  are  shown  only  for  t 
-  0. 


level  at  which  the  transition  from  SCL  to  SCF  conditions  occurred  was 
calculated  to  be  just  one-half  of  the  constant  current  level.  Exper¬ 
imental  measurements  by  others  in  our  laboratory  for  such  a  diode  showed 
such  a  current  break  at  that  level. 


16.2  Switching  into  Avalanche 

When  the  reverse  voltage  applied  for  switching  approaches 
V  ,  the  field  becomes  sufficient  to  cause  ionisation,  and  the  current 
decay  rate  is  reduced  during  its  later  stages.  Figure  100  shows  that 
the  maximum  fields  occur  at  about  one-half  of  the  time  that  it  takes  for 
the  current  to  decay  to  its  SCR  value.  Figure  101  shows  the  current 
decay  for  V  below  V-.  and  the  decay  followed  by  growth  for  V  above 
VB1.  Wiese  calculations  are  for  the  12-ym  PIN  diode  and  R#  »  50  w.  the 
shunt  capacitance,  C,  is  either  lor}*  10“12  F.  Wia  effect  6f  the  RBC 
time  constant  is  dramatic,  since  it  determines  how  fast  the  voltage  can 
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build  up  across  the  diode.  If  the  TIME(m) 

voltage  increases  faster,  the  current 

does  not  drop  as  low  before  the  ava-  Figure  101.  Current  densities 

lanche  threshold  is  reached.  The  during  reverse  switching  into 

threshold  applied  voltage  to  prevent  avalanche  for  12-gm  PIN 

current  decay  is  about  232  V;  above  diode.  Load  resistance  is  50  ft 

that  voltage  an  equilibrium  current  and  parameter  is  applied  volt- 

is  attained.  Calculations  made  with  *  shunt  capacitance  is  1  x 

depletion  (rather  than  charge  stor-  .  £  ^or  solid  curves  and  2 

age)  initial  conditions  have  shown  x  10  r  ^or  dashed  curves, 

either  no  current  growth  or  current 

increasing  indefinitely  until  calculation  instability  halted  the  pro¬ 


gram;  no  equilibrium  could  be  attained.  There  is  considerable  experi¬ 
mental  support  for  the  stabilizing  effect  of  the  role  of  an  "approach" 
current  (of  various  kinds)  to  the  breakdown  transition. 

For  the  highest  applied  voltages  shown  in  figure  101,  the 
current  overshoots  its  equilibrium  value.  At  still  higher  voltages,  the 
oscillation  does  not  decay,  but  builds  up  into  a  MULTIPATT  oscil¬ 
lation.  The  initial  buildup  of  a  MULTIPATT  oscillation  is  shown  in 
figure  102,  where  J  is  plotted  against  V  for  various  applied  voltages. 
The  initial  stored  charge  for  the  calculations  shown  in  figure  102  was 
less  than  that  used  in  those  shown  in  figure  101,  but  the  final  steady- 
state  currents  (when  attained)  are  the  same. 


Figure  102,  t*  'amic  current  den¬ 
sity  versus  voltage  curves  for 
12-yra  PIN  diode.  Applied  voltage 
is  given  as  parameter.  R  ■  50  (2 
and  C  *  1  x  10-12  F.  stored 
charge  (t  =  0)  was  less  for  these 
calculations  than  for  those  of 
figure  101.  Points  along  solid 
lines  indicate  time  intervals  of 
0.1  ns.  Dashed  curve  is  static 
characteristic. 


The  dynamic  impedance 
during  the  switching  transient 
may  be  calculated  by  superimpos¬ 
ing  a  small  ac  signal  (1  or  2  V) 
on  the  applied  voltage.  Since 
the  various  switching  stages  are 
of  nanosecond  or  sub-nanosecond 
duration,  frequencies  of  from  5 
to  20  GHz  have  been  used.  For 
reverse  switching,  the  impedance 
is  low  during  the  charge  storage 
time,  and  the  phase  angle  be¬ 
tween  the  current  and  voltage  is 
found  to  be  about  5  to  10  deg 
capacitive.  The  impedance  and 
phase  angle  increase  during  the 
carrier  sweep  out.  The  dynamic 
diode  capacitance  may  be  deter¬ 
mined  from  the  impedance  and 
phase  angle.  This  capacitance 
is  found  to  agree  with  that  cal¬ 
culated  from  the  width  of  the 
depletion  regions  during  ‘•the  SCL 
sweep  out.  During  the  final  SCF 
sweep  out  of  carriers,  the  phase 
angle  is  about  80  to  90  deg 
capacitance .  The  calculated 
capacitance  is  then  smaller  than 
that  calculated  for  a  fully 
depleted  diode.  This  may  be 
because  of  the  inductive  effect 
of  the  residual  current. 


When  the  diode  is  reverse  switched  at  a  high  enough  voltage 
to  induce  avalanche,  the  phase  angle  becomes  more  inductive.  Hie  calcu¬ 
lated  phase  angle  for  the  12-ym  PIN  diode  is  shown  as  a  function  of 
current  density  in  figure  103.  At  each  frequency  there  is  a  transition 
from  a  capacitive  reactance  (@  =  270  deg)  to  an  inductive  reactance  (0  « 
90  to  110  deg)  as  J  increases.  Hie  transition  occurs  at  higher  J  for 
higher  frequencies.  Hie  calculated  impedance  as  a  function  of  J  is 
shown  in  figure  104,  where  an  area  of  1  x  1 0“"**  cm2  has  been  assumed. 
Hie  impedance  increases  slightly  as  J  increases  and  reaches  a  maximum 
when  the  diode  is  a  pure  negative  resistance  (0  -  180  deg).  Hie  calcu¬ 
lated  inductance  as  a  function  of  I  is  shown  in  figure  105.  Hie  induc¬ 
tance  decreases  approximately  as  I-*'2,  as  found  more  directly  in  other 
calculations  discussed  earlier.  Hiese  calculations  are  for  lower 
current  densities  than  previous  ones;  a  possible  maximum  value  for  L  was 
anticipated,  but  not  found. 
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Figure  103*  Calculated  phase  angle  as 
function  of  current  density  deter¬ 
mined  from  superimposed  ac  signal  of 
frequency  in  G.  labeled  for  each 
curve.  Capacitive,  inductive,  and 
positive  and  negative  resistance 
regions  are  marked. 


Figure  104.  Impedance  as  function 
of  current  calculated  at  frequency 
(GHz)  labeled. 


17.  FURTHER  DISCUSSION 

I  feel  that  it  is  unfortunate 
that  the  term  "breakdown"  is  used  so 
extensively  in  the  semiconductor 
field  to  designate  the  first  appear¬ 
ance  of  current  in  a  reverse-biased 
junction.  Rather  than  the  breakdown 
voltage,  I  would  prefer  the  term 
"avalanche  voltage"  as  more  accu¬ 
rate.  Then  the  awkward  term  "second 
cmmmt  wmitt  (»/■■*)  breakdown*  would  be  avoided.  Break- 

Figure  105.  Inductance  as  func-  down  would  then  be  defined  as  the 
tion  of  current  density  deter-  maximum  voltage  attainable  under  the 
mined  from  figures  103  and  '104,  cited  condition.  Further,  it  would 
assuming  resonance  with  device  be  logical  to  define  the  avalanche 
capacitance.  voltage  in  terms  of  current  density, 

'  rather  than  total  current.  When 

/  device  areas  range  over  more  than  a 

'actor  of  10s,  meaningful  conparisons  can  be  made  only  in  terms  of 
current  density.  It  should  be  mandatory  for  diode  manufacturers  to 
publish  the  area  of  each  device. 
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It  would  also  be  useful  in  predicting  second  breakdown  if  the 
nominal  PTF  of  each  diode  were  available.  Alternatively,  the  doping 
level  and  the  diode  width  could  be  given.  These  data  would  be  useful  in 
predicting  the  ratio  of  second  to  first  breakdown  voltages  and  the 
magnitude  of  the  expected  space-charge  resistance.  Further,  the  doping 
level,  if  available  directly,  or  indirectly  from  the  PTF,  gives  the  best 
criterion  for  the  appearance  of  double  avalanching.  Specifically,  the 
criterion  is  that  the  diode  current  density  should  not  exceed  the 
product  of  the  doping  density,  the  saturation  carrier  velocity  (~1  x  io7 
cm/s  at  300  K) ,  and  the  electronic  charge.  To  use  this  criterion  it  is 
necessary  to  know  both  the  diode  area  and  the  doping  level.  The  appli¬ 
cation  of  the  criterion  to  diodes  of  different  doping  levels  is  shown  in 
figure  106.  It  is  seen  that  the  critical  current  increases  with  the 
doping  level  and  indicates  a  current  level  on  the  static  characteristic 
just  below  the  point  of  negative  resistance.  Since  double  avalanching 
occurs  at  roughly  one-half  the  current  for  a  negative  resistance  and  its 
resultant  current  channeling,  the  criterion  offers  a  safety  margin  for 
the  prediction  of  second  breakdown.  This  criterion  is  not  adequate 
where  the  applied  voltage  pulse  length  is  the  same  or  greater  than  the 
device's  thermal  time  constant.  Then,  lower  currents  over  long  times 
can  cause  damaging  heating.  However,  it  is  the  short-duration  high- 
power  pulse  regime  where  the  experimental  scatter  is  largest  and  the 
thermal  models  are  most  inadequate.  The  proposed  criterion  would  be  of 
both  general  and  particular  use  here.  In  general,  for  low  doping 
levels,  even  low  currents  are  unsafe;  damage  can  occur  when  the  first 
breakdown  voltage  is  barely  exceeded.  A  large  series  resistance  may  be 
called  for.  In  particular,  the  criterion  is  most  useful  if  impedance 
measurements  can  be  made  above  first  breakdown.  Then  the  voltage 
required  to  exceed  the  current  criterion  becomes  known  and  that  can  be 
used  as  an  alternative  criterion  on  voltage  or  power.  Since  the  satur¬ 
ation  velocity  decreases  with  temperature,  the  critical  current  will 
also  decrease  with  temperature,  whereas  the  critical  voltage  will 
increase  with  temperature. 


Figure  106.  Static  characteris¬ 
tics  of  12-pa  N-type  diodes. 
Doping  levels  are  given  as 
parameter.  Horizontal  lines 
indicate  equality  of  mobile 
electron  and  of  doping  dens¬ 
ities,  a  predictor  of  double 
avalanching. 
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Hie  criterion  that  the  mobile  charge  density  not  exceed  the  doping 
level  is  one  that  has  been  introduced  previously.  However,  the  physical 
justifications  of  the  criterion  have  been  faulty.  It  has  often  been 
correlated  with  "avalanche  quenching,"  elimination  of  the  "junction,"  or 
similar  concepts.  Supposedly,  the  removal  of  the  high  field  at  the 
junction  allows  the  thermal  current  to  "run  away."  While  the  present 
calculations  show  that  the  junction  avalanche  multiplication  decreases 
with  temperature,  the  junction  field  actually  increases.  In  any  case, 
avalanche  quenching  would  reduce,  not  increase,  the  total  current.  The 
separation  of  the  high  resistivity  region  of  a  diode  into  junction  and 
bulk  regions  is  also  misleading.  Bulk  regions  should  be  limited  to 
doping  densities  of  ~1 01 9  to  1021  cm'3.  The  high  resistivity  region  of 
the  diode  may  be  correctly  separated  into  depleted  and  undepleted 
regions.  Then  the  undepleted  region  can  be  characterized  by  its  resis¬ 
tivity.  One  must  be  aware  that  the  edge  of  the  depletion  region  moves 
as  a  function  of  V  until  the  edge  reaches  the  high  doping  region. 
Finally,  if  the  "junction"  were  eliminated  at  high  currents,  the 
reverse-biased  diode  would  be  equivalent  to  the  forward-biased  diode. 
Since  the  forward-biased  diode  will  sustain  higher  currents  than  the 
reverse-biased  diode,  a  contradiction  is  evident.  The  presence  of  a 
negative  resistance  current  regime  in  the  reverse-biased  diode  and  its 
absence  in  the  forward-biased  diode  make  the  critical  difference. 

The  other  common  criterion  for  second  breakdown  is  that  the  temper¬ 
ature  increase  causes  the  diode  to  become  "intrinsic."  To  quote 
Kalab:5  "Hie  intrinsic  state  within  the  junction  is  achieved  when 
minority  carrier  density  becomes  equal  to  the  majority  carrier  density 
due  to  a  temperature  rise  within  the  junction.  At  the  critical 
intrinsic  temperature,  it  was  postulated  that  the  device  could  no  longer 
sustain  the  electric  field  stress  necessary  to  maintain  avalanche 
behavior,  and  thus  would  enter  a  second  breakdown  mode  leading  to 
eventual  damage."  Usually,  the  intrinsic  temperature  is  identified  from 
resistivity  versus  temperature  curves  (fig.  3).  As  discussed  earlier 
(sect.  4),  the  low-field  resistivity  curves  are  not  valid  for  the  high- 
field  junction  region  of  the  diode.  Hie  question  has  been  asked,  what 
replaces  the  resistivity  versus  temperature  curves  at  high  fields? 
Unfortunately  there  is  no  replacement.  it  is  necessary  to  use  the 
static  characteristics  at  various  temperatures  for  each  device  doping 
and  width.  Further,  the  intrinsic  temperature  defined  by  the  equality 
of  minority  and  majority  carriers  is  a  confused  picture.  As  discussed 
in  section  10,  thermal  generation  is  insignificant  in  the  high-field 
region,  since  the  recombination-thermal  generation  lifetime  is  three  or 
more  orders  of  magnitude  greater  than  the  carrier  transit  time.  At  very 
high  currents,  the  minority  and  majority  carriers  do  become  nearly 
equal,  forming  a  plasma. 

*B.  Kalab,  Analysis  of  Failure  of  Electronic  Circuits  for  ENP 
Induced  Signals— Review  and  Contribution,  Harry  Diamond  Laboratories, 
HDIs-TR-1615  (August  1973). 
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The  analogy  between  the  glow  and  arc  transitions  in  gases  with  the 
first  and  second  breakdowns  in  semiconductors  has  often  been  made.  I 
believe  this  analogy  is  more  complete  than  has  been  postulated.  First, 
all  the  basic  equations  used  in  the  computer  program  are  essentially  the 
same  for  semiconductors  and  gases,  except  for  the  pressure  variation  in 
gases,  the  doping  level  in  semiconductors,  and  the  form  of  the  recombin¬ 
ation  rate  expressions.  The  doping  level  in  the  semiconductor  is  in 
many  ways  an  inverse  analog  of  pressure,  p  ,  in  the  gas.  Gas  breakdown 
is  a  function  only  of  the  p  d  product,  whereas  in  the  semiconductor  it 
is  a  function  of  1/N  (eq  (26))  for  d  large  and  of  d  (eq  (27))  for 
intrinsic  material.  In  particular,  the  field  dependence  of  the  ioni¬ 
zation  coefficients  (eq  (12),  (13))  and  of  the  mobilities  (eq  (14), 

(15))  is  the  same  in  gases  and  semiconductors.  In  contrast,  a  major 
difference  is  that  ions  in  a  gas  are  orders  of  magnitude  less  mobile 
than  are  electrons,  and  consequently  the  same  difference  arises  in  the 
ionization  coefficients,  whereas  each  is  of  the  same  magnitude  in  semi¬ 
conductors  . 


As  is  the  case  in  semiconductors,  breakdown  in  gases  has  been 
defined  in  many  ways:  appearance  of  a  small  fixed  current,  a  certain 
large  multiplication  factor,  attainment  of  a  maximum  voltage,  and 
others.  Studies  of  breakdown  and  of  the  glow  and  arc  transition  in 
gases  have  major  advantages  over  those  in  semiconductors:  visibility, 
larger  physical  dimensions,  and  a  slower  time  scale.  Although  details 
are  still  controversial,  the  basic  glow  and  arc  transitions  are  well 
understood.  The  contraction  into  the  glow  is  the  consequence  of  an 
electrical  negative  resistance  and  the  transition  into  the  arc  is  a 
thermal  electrode  process.  How  close  is  the  semiconductor  analogy? 


One  test  of  the  analogy  is  to  compare  field  and  charge  distribu¬ 
tions  for  the  gaseous  glow  discharge  and  those  of  the  semiconductor 
analog.  Hie  400  K  characteristics  of  the  simulated  1N4148  diode,  shown 
in  a  log  plot  in  figure  62,  are  replotted  on  a  linear  scale  in  figure 
107.  As  previously  indicated,  IMPATT  oscillations  occur  at  both  low  and 
high  current  densities,  but  may  be  damped  out  with  the  use  of  higher 
diffusion  coefficients.  Field  distributions  corresponding  to  the 
numbered  points  in  figure  107  are  shown  in  figure  108,  and  the  corre¬ 
sponding  charge  distributions  are  shown  in  figure  109.  In  the  negative 
slope  region,  the  field  and  charge  distributions  are  remarkably  like 
those  observed  in  gaseous  glow  discharges.  In  particular,  cathode  fall, 
plasma,  and  anode  fall  regions  may  be  distinguished.  The  function  of 
the  cathode  fall  (high  field  at  the  PN  junction)  is  to  inject  electrons 
(through  avalanching)  into  the  plasma.  A  smaller  anode  fall  (high  field 
at  the  NN+  junction)  is  required  to  inject  a  smaller  number  of  holes 
into  the  plasma. 
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Figure  107.  Characteristics  of  12-pm 
N-type  diode  doped  to  3  x  1015  cm”3 
at  400  K.  These  same  data  are 
plotted  on  a  log  scale  in  figure 
63.  Ellipses  indicate  presence  of 
IMP ATT  oscillations;  50-fl  load  line 
for  a  diode  area  of  5  x  10”5  cm2  and 
Va  =  350  V  is  shown. 

Dumin14  *  has  observed  visible 
radiation  from  plasmas  in  silicon 
diodes  during  second  breakdown. 

He  found  that  the  junction  was 
necessary  to  ignite  the  plasma. 

At  the  highest  currents  the  aver¬ 
age  field  was  quite  low.  Visible 
radiation  during  second  breakdown 
has  also  been  observed  at  Harry 
Diamond  Laboratories*  Portnoy  and 
Gamble1*2  found  two  or  more  stages 
(voltage  steps)  in  second  break¬ 
down.  It  is  quite  likely  that 
they  observed  the  effect  of  one  or 
more  avalanche  oscillations.  As 
voltage  changes  when  oscillations  a 
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Figure  108.  Field  distributions  for 
current  density/voltage  points  as 
numbered  in  figure  107. 


Figure  109.  Distribution  of  elec¬ 
tron  and  hole  densities  for 
numbered  points  of  figure  107. 
Solid  curve  gives  doping  profile; 
dashed  lines  are  for  electrons 
and  dotted  lines  for  holes. 

shown  in  figures  63  and  107,  the 
pear  or  change  modes. 


I*1D.  J.  Dumin,  Emission  of  Visible  Radiation  from  Extended  Plasmas  in 
Silicon  Diodes  During  Second  Breakdown,  IEEE  Trans .  Electron  Devices, 
ED-16  (May  1969),  479-485. 

**2IF.  M.  Portnoy  and  F .  R.  Gamble,  Fine  Structure  and  Electromagnetic 
Radiation  in  Second  Breakdown,  IEEE  Trans.  Electron  Devices,  ED-11 
(October  1964),  470-478. 
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Gaseous  glow  discharges  may  be  classified,  in  the  order  of 
increasing  current,  as  subnormal,  normal,  or  abnormal.  The  subnormal 
discharge  is  in  the  negative  resistance  region;  it  is  usually 
oscillating.  The  normal  discharge  has  a  constant  J;  the  area  changes  as 
the  applied  voltage  is  increased.  When  the  cathode  is  fully  covered,  a 
higher  voltage  is  required  to  increase  the  current  density,  resulting  in 
the  abnormal  glow.  There  are  strong  indications  of  a  minimum  voltage 
near  the  highest  currents  calculable  in  some  of  the  present 
semiconductor  computer  calculations.  Hie  decreasing  slope  (from  the 
vertical)  of  the  characteristic  is  noted  at  the  highest  currents  shown 
in  figure  107.  During  oscillations,  pulse  currents  can  reach  a  much 
higher  magnitude  than  can  be  calculated  in  static  conditions.  When  the 
current  peaks  of  the  dynamic  characteristics  of  differing  applied 
voltages  are  joined,  the  voltage  shows  a  minimum  at  a  current  slightly 
above  that  calculable  for  the  static  case  and  increases  considerably  for 
higher  currents.  As  shown  in  figure  63,  the  current  maxima  and  minima 
during  oscillations  are  on  the  static  characteristic. 

The  glow  to  arc  transition  in  a  gaseous  discharge  tube  occurs  as  a 
result  of  cathode  heating.  Thermal  emission  replaces  the  cathode  fall 
as  an  electron  source.  In  the  same  way,  the  thermal  injection  current 
(sect.  10)  in  the  semiconductor  can  negate  the  need  for  a  high  junction 
field.  In  the  gaseous  discharge,  there  is  an  exact  correlation  with  the 
appearance  of  a  negative  resistance  and  the  contraction  of  the  visible 
discharge.  Current  densities  can  be  calculated  from  the  observed  dis¬ 
charge  area.  The  measured  J-V  characteristic  then  agrees  with  the 
calculated  characteristic  from  breakdown  to  the  abnormal  glow.  in  the 
s ami conductor,  the  channeling  cannot  be  observed  directly  in  conven¬ 
tional  devices,  but  is  evidenced  by  melt  channels  in  damaged  devices. 

Relaxation  oscillations  in  the  gaseous  subnormal  glow  discharge  are 
well  known.  Section  14  shows  that  they  are  observed  in  silicon 
diodes.  Moving  striations  are  generally  present  in  the  positive  column 
(plasma)  of  a  gaseous  glow  discharge.  While  these  striations  are  not 
completely  understood,  they  have  been  proven  to  arise  in  the  cathode 
and/or  the  anode  fall  regions.  There  is  an  almost  complete  analogy  tc 
what  has  been  named  the  traveling-wave  IMPATT  mode  (Ward25).  One  must 
conclude  that  the  analogy  between  the  gaseous  glow  and  arc  discharges 
and  the  steps  of  the  second  breakdown  process  is  a  most  useful  one. 


18.  SUMMARY  AND  RECOMMENDATIONS 


18.1  Summary 


After  a  brief  historical  introduction  to  second  breakdown,  an 
electrothermal  model  computer  program  is  described.  The  formulation  and 


25A.  L.  Ward,  Modes  of  Avalanche  Oscillations  in  Silicon  Diodes,  IEEE 
Trans •  Electron  Devices,  ED-25  (June  1978),  683-687 . 
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procedures  are  listed  and  discussed.  The  material  parameters  of  silicon 
are  described.  “Hie  various  current  regimes  of  the  reverse-biased  static 
characteristics  are  discussed.  It  is  shown  that  space  charge  causes 
double  avalanching,  which  in  turn  results  in  a  negative  differential 
resistance  and  current  channeling.  Dynamic  characteristics  are 
presented  for  constant  overvoltages  and  for  constant  temporal  rates  of 
voltage  increase.  Both  isothermal  and  adiabatic  thermal  calculations 
are  made.  It  is  shown  that  thermal  injection  currents  are  the  critical 
thermal  process.  One-sided,  two-sided,  and  linearly  graded  doping 
profiles  are  compared.  Although  direct  calculations  are  limited  to  the 
nanosecond  time  scale,  it  is  shown  that  longer  times  may  be  simulated  by 
combining  knowledge  of  the  temperature  variation  of  static  character¬ 
istics,  the  rate  of  temperature  rise  as  a  function  of  current,  and  the 
thermal  time  constant  of  the  diode.  Good  agreement  has  been  obtained 
with  published  measurements  of  power  as  a  function  of  breakdown  time  up 
to  100  ps  for  the  1N4148  diode.  The  combined  effect  of  ionizing  radi¬ 
ation  and  high  electrical  stress  is  calculated.  It  is  shown  that  the 
frequency  of  avalanche  oscillations  can  be  used  as  a  dynamic  temperature 
monitor  of  the  junction  region.  High-current  forward-bias  static  char¬ 
acteristics  are  calculated.  It  is  shown  that  the  higher  power 
capability  of  the  forward-biased  diode  results  from  the  absence  of  a 
negative  resistance.  Reverse  switching  transients  into  reverse  break¬ 
down  are  calculated.  The  analogy  of  second  breakdown  processes  in 
semiconductors  with  the  glow  and  arc  discharges  in  gases  is  discussed 
and  found  to  be  reasonably  accurate. 

18.2  Recommendations 


Until  recently,  the  experimental  work  on  second  breakdown  far 
outpaced  the  theoretical  studies.  Now,  the  limitations  on  the  available 
measurements  are  handicapping  theoretical  progress.  This  has  resulted 
from  the  lack  of  pertinent  information  in  the  experimental  publi¬ 
cations.  Very  few  papers  furnish  the  diode  area,  doping  profile  (or 
even  the  doping  level  or  diode  width),  and  the  thermal  time  constant  (or 
the  thermal  resistance)  needed  for  calculations.  Close  cooperation  is 
needed  between  the  theoretical  and  experimental  work  in  second  breakdown 
in  order  to  expedite  solutions  to  diode  damage  susceptibility  problems. 

It  is  planned  to  continue  the  exploitation  of  the  elec¬ 
trothermal  model  computer  program  in  simulating  as  wide  a  variety  of 
problems  as  possible.  Feasibility  studies  should  be  made  in  further 
improvements  to  the  program.  Of  highest  priority  is  the  investigation 
of  using  larger  time  steps  for  thermal  processes  whenever  the  electrical 
processes  are  in  quasi-equilibrium.  Second  priority  should  be  given  to 
extending  the  calculations  to  two  dimensions.  This  would  allow  quanti¬ 
tative  calculations  of  current  channeling  in  the  negative  resistance 
regime  and  would  allow  nonplanar  geometries  to  be  studied. 
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At  present,  the  time  scale  is  extended  by  hand  calcu¬ 
lations.  This  could  be  programmed  for  computer  calculation,  either 
directly  or  by  circuit  simulation  codes.  Oily  lack  of  money  and  person¬ 
nel  has  delayed  this  important  step  in  developing  an  engineering  model; 
it  should  be  developed  as  soon  as  possible. 
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SYMBOLS  AND  ACRONYMS 


A 


Ai 


Bi 


L 

M 


ND  *  NA 


P 


R,  R' 


's'  Rd 


T 


Diode  area  (cm) 

Constants,  electron  ionization 
Constants,  hole  ionization 
Capacitance  ( P) 

Constants,  mobility 
Diffusion  coefficient  (cm2/s) 

Thermal  diffusivity  (cm2/s) 

Electric  field  (V/cm) 

Band  gap  of  semiconductor  (ev) 

Fraction,  defined  by  equation  (8) 

Total  current  (A) 

Total  current  density  (A/cm2) 

Boundary  current  density  (A/cm2) 
Electron,  hole  current  density  (A/cm2) 
Distance  to  heat  sink  (cm) 

Inductance  (H) 

d/Ax,  multiplication  factor 

Net  doping  level  (cm-3) 

Power  (W) 

Recombination  rate 

Series  resistance,  diode  resistance  (Q) 
Temperature  (K) 

Voltage  (V) 

Applied  voltage  (V) 
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t.  At 


x.  Ax 


V+.V_ 


SYMBOLS  AND  ACRONYMS  (Cont'd) 

First  breakdown  voltage  (V) 

Second  breakdown  voltage  (V) 

Maximum  depletable  width  (cm) 

Heat  capacity  (J/g  C) 

Diode  width  (cm) 

Depletion  width  at  breakdown  (cm) 
Electron  charge  (Coulombs) 

Boltzmann  constant  (J/C) 

Thermal  conductivity  (W/cm  C) 

Index 

Electron  density  (cm) 

Intrinsic  density  (cm) 

Hole  density  (cm) 

Time,  time  step  (s) 

Carrier  velocity  (cm/s) 

Distance,  distance  step  (cm) 

Ionization  coefficient,  electrons  (cm) 
Ionization  coefficient,  holes  (cm) 
Permittivity  of  free  space  (P/cm) 
Thermal  resistance  (C/M) 

Phase  lag  (deg) 

Dielectric  constant 
Growth  constant 
Low-field  mobility  (cm /Vs) 

Mobility  of  holes,  electrons  (em/ve) 
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p 

o 

T 

V 

Tth 

BMP 

NDR 

PIN 

PTF 

SCF 

SCL 

RAM 


SYMBOLS  AND  ACRONYMS  (Cont'd) 

Density  (g/cm) 

Electrical  conductivity  (A/V  cm) 

Dielectric  relaxation  time  (s) 

Recombination  lifetime  (s) 

Thermal  time  constant 
Electromagnetic  pulse 

Negative  differential  resistance  (resistivity) 

Diode  with  intrinsic  region 

Punch-through  factor 

Space  charge  free 

Space  charge  limited 

Relaxing  (or  relaxation)  avalanche  mode 
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